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I. INTRODUCTION 


A. BACKGROUND 


The Joint Army/Navy Rotorcraft Analysis and Design (JANRAD) computer code 
was originally developed in response to the 1993 American Helicopter Society (AHS) 
Design Competition. The emphasis on the use of locally developed analysis tools as 
opposed to packaged commercial software was an integral part of the Naval Postgraduate 
School design course philosophy. This code has been instrumental in the NPS Helicopter 


Design Teams winning consecutive 1* or 2" place awards in the last 5 annual contests. 


Originally developed as a preliminary design tool, the main performance function 
calculated 25 different helicopter performance parameters given 35 input parameters. 
These parameters are described in more detail in Nicholson (1993). The program was 
written using MATLAB” PC version 3.5 using a combination of script and function M- 
files. Inputs were made at the command line or loaded as a filename.mat file from the 
current working directory. The output was displayed in the main workspace and an 


option to save the output was also built into the program. 


Since JANRAD version 1.0 was initially developed, many features have been 
added. A routine to analyze Stability and Control was added by Wirth (1993) and the 
latest, a section for Rotor Dynamics, was developed by Hiatt (1995). Minor updates have 
included time varying tip loss analysis, the ability to calculate performance characteristics 


over a range of selected input parameters and the addition of airfoil choices. 


Eccles (1995) was the lastest thesis written concerning JANRAD. His thesis 
validated the code by a detailed comparison with Sikorsky UH-60A and H-34 measured 
flight test data with calculated results from JANRAD version 3.1. It was determined that 
JANRAD results predicted power required within 2% for altitudes below 6000 feet MSL. 


These results were considered very good in light of the fact NACA 0012 airfoil data was 
used instead of actual Sikorsky UH-60A airfoil data. 


B. JANRAD 98 


During the last two AHS design competitions and numerous projects within the 
helicopter curriculum at NPS, students initially unfamiliar with the program described the 
interface as cumbersome and non-intuitive. Once familiar with the program interface 
however, the tool quickly demonstrated its usefulness. Since the program was difficult to 
run on the outset, during the design project, students would designate a specific team 
member to learn and run the JANRAD program throughout the design process. Because 
of the compressed nature of the curriculum, this was not the most efficient use of 
manpower and highlighted the need for an update of JANRAD. A concensus of former 
team members agreed that a graphical user interface (GUI) similar to Windows” would 
improve the overall look, feel and speed of the program. It would also allow any student 
to quickly learn and run models without having to refer to detailed instructions or an 


experienced user. 


Additionally, there have been major changes to MATLAB”, the language 
originally used to develop JANRAD. Although graphical user interface features existed 
in prior versions, these features had to be coded line by line. Each user interface control 
such as a push button on a figure window could easily require 10 to 20 lines of 
MATLAB*® code depending on the complexity of the interface function. MATLAB® 
version 5.0 included a new function that would automatically generate the required code 
while the programmer "Dragged and Dropped" controls or objects from a pallet to a 
blank MATLAB® figure window. The function called GUIDE” (Graphical User 
Interface Design Environment) maked the process of creating a GUI quicker and easier. 
The programmer now spent the major portion of his time with the task of tying the GUI 
to the analysis part of the program. 


The MATLAB” Editor/Debugger was also instrumental in creating the GUI. This 
editor program automatically indented and highlighted specific MATLAB” commands 
with different colors. These features allowed easier and quicker identification of code 


and simplified the extensive modifications required in developing JANRAD 98. 


With the necessary tools made available, JANRAD transitioned from a command 
line program, reminiscent of MS DOS, to the look and feel of some of the latest 
professional engineering software packages. It also gave the user flexibility and quick 
response to changing design specifications. The new JANRAD version 4.0 was named 


JANRAD 98. 
OF NEW FEATURES 


A basic Users Guide, included in Appendix A, describes an example run using 


UH-60 input data and incorporates important new features over version 3.1: 


e Any or all input parameters can be changed and run without saving or altering 
an existing file. 

e Input and output parameters are all displayed on one window. 

e User can retrace steps back through program to verify selections or options. 

e User can exit or quit program at any point of analysis. 

e Actions such as saving or printing input/output files are executed by a single 
mouse click. 

e Creating a new file does not require prompting for each input. Values are all 
typed directly into edit boxes on a single input window. 

e A total of 5 airfoil choices are now available. 


e User interface minimizes keyboard entries. 








ll. DESIGNING THE FRONT END 


A. WHY CREATE A GRAPHICAL USER INTERFACE? 


"Why would anyone want to create a GUI in MATLAB®?” The short answer 
according to Hanselman/Littlefield (1996) is probably no one. However, an engineer or 


engineering student might want to develop a GUI when: 


e Writing code that will be use often or repetitively ; or 


e When writing code that will be used mostly by other engineers or students. 


Both of these reasons applied to JANRAD. It is also important to note that students have 
developed JANRAD primarily with other students in mind. This concept was also true 
for the JANRAD 98 GUI development. 


B. PAPER PROTOTYPE 


The first and probably most important step in developing the JANRAD 98 
interface was construction of a paper prototype. Dean (1997) describes and summarizes 
the prototyping procedure as a low tech design strategy leading to quick iterative GUI 
designs. The design is done on paper using pens, tape, transparencies etc. Usability tests 
are run on the paper prototype by simulating the actions of the GUI using these props. 
By shuffling these props in response to user actions, the prototype logic is tested before 


any code is written. 


The paper prototype used for JANRAD 98 was drawn on engineering paper using 
colored pencils and straight edge. The prototype showed the general layout of each 


proposed window and the function of the controls. At this point, no MATLAB” code had 


actually been written for the GUI. Because very little was known about the capabilities 


or limitations of MATLAB® GUI's, the prototype started with the ideal design. 


Usability tests were conducted by allowing a user to physically touch buttons on 
the paper prototype as if their finger was the mouse pointer. The individual 
windows/controls were then moved to indicate the next action of the GUI. Figure 1 from 


The MathWorks” shows this process with a more detailed prototype. 
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Figure 1. Paper Prototype Usability Testing at the MathWorks°. 


The overall effect of paper prototyping allowed the programmer and testers to 
suggest and make changes quickly and without hesitation. The result was a better initial 
design and more confidence in the initial phase of development. It was noted that 
JANRAD 98 implemented 100% of the features contained in the original paper prototype 
plus several additional capabilities. Building GUI’s with MATLAB (1997) detailed 
much of the GUI design principles used throughout this project. 


C. USING MATLAB” GUIDE” 


MATLAB®’s Graphical User Interface Design Environment (GUIDE), included 
in version 5.0, was the primary tool used to develop JANRAD 98. The GUIDE® function 
is basically a GUI used to create GUI’s. It is started from the command line by typing 
guide. Figure 2 depicts the GUIDE® Control Panel. In addition to the control panel, a 


blank MATLAB® figure window will also appear on start up. As shown in Figure 2, the 
ten different user interface controls are displayed along the bottom of the panel. These 
controls are generically referred to in MATLAB” code as uwicontrols, meaning user input 


controls. 
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Figure 2. The GUIDE” Control Panel from MATLAB? 5. 


By clicking on one of the uicontrols and minimizing the control panel, the 
programmer can then place the object on the blank figure window. After the desired 
uicontrols have been placed on the figure, they can be sized, moved and aligned using the 
mouse or other GUIDE” tools. Figure 3 shows a JANRAD 98 GUI window under 
construction using GUIDE®. 
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Figure 3. JANRAD 98 GUI Window Development Using GUIDE” . 





GUIDE® will prompt the programmer to save the figure window when finished. 
At this point GUIDE® has machine generated all of the code necessary to redraw the 
figure when the file is executed. GUIDE” saves two files for each GUI figure. The first 
is a MATLAB?” function with extension .m, and the second is a data file with the same 
filename and .mat extension. To recreate the figure, the file name is typed at the 
command line. The first line of the function will load the filename.mat. Then the 
remainder of the code will be executed and the figure is placed above the base 


workspace. 


GUIDE® was also used to customize the figure by selecting specific properties for 


each uicontrol. These uicontrol properties determine the size, shape, color, position, 


label, and even font size and weight among other things. Properties not only determine 
the look of the control, but also its functions. The most critical of these is the Callback 
property. The Callback property contains the MATLAB” commands that will be 
executed each time the uicontrol is activated. Any code that can be executed from the 
command line can be executed from the Callback. Two other related properties are the 
CreateFcn and DeleteFcn. These properties are similar to Callback except they will 
execute the commands they contain only once, when the uicontrol is either first created or 
deleted respectively. It is important to note here that although each uicontrol is created in 
a function workspace, the Callback, CreateFcn and DeleteFcn are executed in the base 
workspace. This property and its ramifications will be discussed in more detail in the 


next chapter. 


GUIDE® generated the code for individual uicontrols in the order they were 
created. This often made finding a particular line within the file difficult. Particulary 
when doing detailed editing. In addition to the difficulty in finding specific code, while 
operating the GUI, the Tab key followed the order that the controls appear in the M-file. 
The decision was made to cut and paste the completed GUI uicontrols for all JANRAD 
98 files to a more logical order. This allowed the user to use the tab key to navigate 
through the GUI quicker and easier. It also allowed the programmer to identify specific 


code wihin M-files when the GUI program was connected to the computational end. 


Using GUIDE® to edit figure windows had some limitations worth noting. 
Specifically, when editing an existing figure window, subsequent saving deleted any 
customized Graphic Handles. Hanselman/Littlefield (1996) described Graphic Handles 
as variables attached to particular uicontrols. These handles were used extensively to 
change uicontrol properties during the operation of the GUI. Because of this limitation, 
large-scale changes to the figure window were very difficult. Small changes also 


necessitated a strong working knowledge of uicontrol code and Handle Graphics. 


Hanselman/Littlefield (1996) and Marchand (1996) contain the most in depth discussions 
on the subjects of Handle Graphics. 


Il. CONNECTING THE BACKEND 


A. FILES AND VARIABLE CONVENTIONS 


To keep track of the extensive number of M-files and variables it was important to 
establish a convention for naming these files and variables. Many of the new files 
created for the interface were given longer and more descriptive file names to minimize 


confusion. The following lists the file name conventions: 


e File names beginning with an Upper case first letter indicates a computational 
file. Example: Perfim, Trim.m, Thrcalc.m, Tmcalc.m and Dmcalc.m 

e File names in all lower case letters indicates a GUI file. Example: janrad98.m 

e File names with lower case and _fcn.m endings indicate Switchyard Callback 


Functions. Example: janrad98_fcn.m 


Because of the increased size and complexity of JANRAD 98, some conventions 
had to be developed to manage the additional variable names and Graphic Handles. 


These conventions are as follows: 


e Variables in all UPPER CASE letters indicate global. Example: ASPECT 

e Variables beginning with an H_ indicate a Graphics Handle. Example: 
H_CHECK1 indicates a global graphics handle for a check box. 

e Variables beginning with an S_ indicate a structure variable. Example: 


S_USER_INPUT indicates a global structure. 


Although not a strict convention, Graphic Handle variables were sometimes 


selected as the initials of the label or string that appeared next to the uicontrol in the GUI. 
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Particularly, this method was used on most radio buttons. For example, the handle for the 


Create New File radio button was H_CNF, indicating a global graphics handle. 
B. THE SWITCHYARD CALLBACK 


Dean (1997) describes in detail the concept of the Switchyard Callback. This is a 
MATLAB® programming technique that increases speed and avoids the proliferation of 
function files. As discussed earlier, the Callback property lists MATLAB® code that will 
be evaluated when the uicontrol is pushed or otherwise activated. Shorter Callback 
function names run faster than long lists of commands or functions. Additionally, 
complex Callback strings will be harder to de-bug and modify if they are nested within 
the GUI code itself. A Switchyard function is a single function that 1s called by most or 
all of the controls within a single figure window GUI. In other words, every control on a 
single window will call the same function — the Switchyard function. The Switchyard 
function contains a MATLAB? switch to determine which control was activated and then 


execute the applicable code. 
CG: STRUCTURES 


MATLAB® version 5 incorporates new variable types, one of which is called a 
Structure. Simply stated, a structure is an array of data containers referred to as fields. 
The struct command is used to define the structure. The result is the ability to group a 
large number of variables under one name to pass to a function, assign to a workspace or 
declare global. This feature was use in several areas of JANRAD 98. First, the 36 
separate user input variables were stored in a single structure called S USER_INPUT. 
Then, after the computational routines were completed, the 25 output values defined in a 
Structure called S PERF OUTPUT. Both of these structures were declared global to be 
available to multiple functions and Callbacks. The variables within a structure can be 


accessed individually or as a group. They can be displayed by typing the structure name 


at the command line. In JANRAD 98, the input structure is saved in a filename. mat file, 
which can be loaded and edited for subsequent JANRAD evaluations. When creating a 
new input file, JANRAD creates the structure after the user presses the continue button in 


the input window. 
D. MATLAB” WORK SPACES 


Because functions work in their own workspace, writing complex programs with 
function M-files is more efficient. The single most difficult part of the interface 
development was ensuring all the necessary variables and graphics handles were 


available to the GUI Callbacks and Switchyard functions. 


Initially, the plan to manage the architecture of the program was to write the 
entire code as functions. This included the computational files originally written as script 
files. As the program took shape however, this proved difficult. Because the Callback 
property works in the base workspace of MATLAB”, it became apparent that simply 
passing groups of variables from one function to the next did not work. Additionally, the 
Callbacks needed access to many of the graphics handles defined within the function 
workspace. The fix as described in Marchand (1996) required that all of the graphics 
handles be declared global. Many of the graphics handles, specifically the radio button 
handles, also needed to be assigned to the base workspace via the assignin command. 
The input and output structures were designated global to enable the user maximum 
flexibility to navigate freely within the GUI while being able to modify and store the 


input values throughout the process. 


These global variables were stored in their own separate MATLAB” workspace. 
Once designated a global, a variable was called at the beginning of the script or function 
file in order to be defined. However, if a global was initially declared within a function, 


the value of that global became empty once that function completed its execution. In 
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other words, it was still a global variable, but its value was erased after the function 
completed its execution. This property required careful attention to when and where 
global variables were declared. To avoid losing the value of a global variable, some 
cases required the global variable to be declared within the Callback before the 


Switchyard function was called. 


The computational files Perf'm, Trim.m, Thrcalc.m, Tmcalc.m and Dmcalc.m 
were called from within a Switchyard Callback Function. Because they worked within a 
single function workspace and shared many variables, they remained script files without 
significant performance degradation. The airfoil functions however, were passed a single 
angle of attack and a single Mach number to return a lift coefficient and drag coefficient. 


Therefore, they remained functions. 


E. ARCHITECTURE 


In general, the basic layout required the user to first start the GUI from the 
command line of MATLAB® by typing janrad98. The first GUI function executed 
created the first window. After making the appropriate selections, the Switchyard 
function changed or modified variables, closed the existing window if necessary, called 
the next GUI function and executed any other necessary script or function files. As the 
user progressed through the program, input structures were built, options saved and 
performance computation files were executed, all within the switchyard callback 


functions. 


The large number of files and functions (40+) required careful management of the 
GUI procedures. A simple flow chart was constructed to track the files and Callback 
operations and is included in Figure 4. This chart describes the basic process of starting, 
selecting a file and running the performance routine through one complete cycle. To 
simplify the figure, it does not depict all possible progressions through the file structure. 
Appendix D through AP include all of the M-files developed for JANRAD 98. 


The under lying principle in helicopter rotor design is that the main rotor forces 
and moments must be adjusted so the rotor provides the lift and propulsive thrust required 
to meet the conditions of flight. The process of adjusting the rotor’s forces is referred to 
as “trimming” the rotor. The performance module uses a routine which employs blade 
element theory and a harmonic balance method for rotor tim. This method was 
described in detail by Nicholson (1993). The flow chart in Figure 5 graphically shows 


the analysis method used to complete the performance analysis. 
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MATLAB COMMAND LINE 


>> janrad98 











stability_control.m 


performance_output_fcn.m 


iteration method.m 


Airfoil Airfoil Airfoil 
Functions Functions Functions 


: 
3 


Figure 4. JANRAD 98 File Structure Flow Chart. 
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Figure 5. JANRAD Computation Flow Chart . 
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IV. VALIDATION 


A. JANRAD 3.1/JANRAD 98 OUTPUT 


Identical input parameters were run in versions 3.1 and 4.0 to validate the output 
of JANRAD 98. Because the format of the input was changed, a single input file to be 
run in either program was not possible. Therefore, two separate files were built with 
identical values for the input variables. As expected, the two programs calculated 
identical outputs. These results are contained as actual JANRAD 3.1 and JANRAD 98 
printed output in Appendix B and C. 


There were however, suspicions about the origin of the JANRAD 3.1 files that 
were modified for use in JANRAD 98. Because a master copy of JANRAD 3.1 was not 
available, a previously modified but undocumented version was wail It was then 
decided to plot JANRAD 98 power required verses airspeed and compare with flight test 
results as done in Eccles (1995). Figure 3 showed close correlation between JANRAD 
98 and flight test data from a National Aeronautics and Space Administration UH-60A 
Test Flight. Again, as in Eccles (1995), VR-12 airfoil data was used in order to check the 
performance of the computational routines with respect to JANRAD 3.1. Figure 6 was 
then compared to the JANRAD 3.1 results in Figure 46 of Eccles (1995). They showed 
essentially identical plots. However, the exact JANRAD 3.1 inputs used to create Figure 
46 of Eccles (1995) were not specifically documented. Therefore, this check served only 
as a general indication of version 4.0 reliability relative to version 3.1 and was not 


intended to substitute for a detailed code validation. 
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UH-60A Power Required vs. Airspeed, Flight #85 










2800 J ____ a ain Seed aie 
| [ [ | 

2600} - - --- rene ere is ee 
[ | [ [ 

2400 a eee ee ne 5 

2200+ ----- | _|G@——O_  JANRAD 98 Results 

00 “-|+ + ‘Test Flight #85 Results 
2000) ie ee a i 


> ak 
a oO a 
oO Oo oO 
oO AF @®S 


Main Rotor Power Required (HP) 


0 20 40 60 80 100 120 140 160 
True Airspeed (KTS) 


Figure 6. UH-60A Power Required vs. Airspeed, Flight #85. 


B. USABILITY TESTS 


Although a formal usability test was not performed, it was noted during the 
validation procedure that JANRAD 3.1 lacked robustness. During several JANRAD 3.1 
sessions, the program locked up in what appeared to be an infinite loop. The program 
seemed prone to this condition particularly when editing multiple input parameters and 
subsequently saving these files. This was probably due to the large number of variables 
and evaluations taking place in the base workspace. Clearing the base workspace and 
restarting JANRAD usually corrected these problems. Other students have recently 
experienced similar problems. As discussed earlier, however, the JANRAD 3 program 


used may have undocumented modifications that caused these complications. 
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JANRAD 98 on the other hand, demonstrated superb robustness. During a single 
continuous session, the program ran several new and edited files as well as several input 
iteration runs without any apparent problems. Saving and printing operations performed 
well. Additionally, running several consecutive inputs proved much quicker and easier. 
The program could be toggled from output back to the input with a click of the mouse 


button. Subsequent computations were started with two additional mouse clicks. 
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V. RECOMMENDATIONS 


Having demonstrated the usefulness of the graphical user interface, there are 
unlimited possibilities to improving the look, feel and features of JANRAD 98. 
Recommendations are subdevided into three areas: (1) GUI recommendations; (2) 


analysis recommendations; and (3) general recommendations. 
A. GUI RECOMMENDATIONS 


First, the GUI should continue to be modified according to student/user 
preferences. If JANRAD 98 is being used differently than the author’s original intent, 
these changes should be reported to the Helicopter Design Team instructor, who inturn 
will document the necessary modifications. Enabling and disabling uicontrols within the 
program will need careful attention. Although every attempt has been made to prevent 
the user from making illogical selections within the GUI, it is impossible to eliminate all 
possible errors. If any of these errors are identified they should be corrected by disabling 
the appropriate uicontrols or the development of warning windows to notify the user of 


the limitations of the program. 


Presently, only the airspeed iteration method will plot a series of performance 
calculations with respect to the airspeed variations. The addition of plotting routines for 
the other iteration methods would allow for a more compete analysis, specifically when 
trade-off studies are desired. An additional window allowing the user to select 4 or 5 
specific pre-determined plots or subplots via check boxes would improve the program 
greatly. An ‘Other’ check box could refer the user to an empty script M-file to modify 


for custom plots of special interest without having to modify the performance files. 
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Because most student engineers would not be familiar with shapes of particular c; 
and cq curves for various airfoils, saved plots of these curves, available for viewing, 
would enhance the users ability to pick an appropriate airfoil. These curves could be 
plotted and their resulting figures saved to avoid having to execute a script file each time 
an airfoil selection is examined. The selection could be made though a separate pull 


down menu titled Airfoil Plots appearing above the input window. 


B. ANALYSIS RECOMMENDATIONS 


The analysis routines developed for JANRAD 3.1 and used in JANRAD 98 were 
the same ones validated in Eccles (1995). Although the results were generally good, 
there were some unexplained variations at altitudes above 6000 feet. With the role of 
helicopters expanding to include high density altitude operations, this limitation is 
significant. A similar analysis using the Sikorsky $C1094r8 and SC1095r8 airfoil data 
currently being added to JANRAD 98 should be completed to gather more insight to this 


documented discrepancy. 


Additional airfoil data from a variety of manufactures and type/class helicopters 
would allow users to better match the mission and type helicopter with appropriate 
airfoils. As a minimum, data should be included for a large transport, military attack, and 


a light civil helicopter. 


The performance routines currently assume the center of gravity of the aircraft to 
be located directly below the rotor mast. The capability for the user to offset the center of 
gravity from the rotor hub would greatly enhance the capabilities of the program. This 
capability should be included as a single user input as well as an iteration parameter to 


study the effects of CG travel on performance and stability and control. 
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A 3-D plot of airloads or induced velocity at selected azimuth intervals or at least 
at y = 0°, 90°, 120° and 180° should be produced as part of the output. These plots would 
improve the validation of each successive run. Should a trim routine not properly 


converge or have other errors, it would be immediately apparent by the shape of this plot. 
c. GENERAL RECOMMENDATIONS 


First, master copies of JANRAD 98 should be maintained and clearly marked. 
Any subsequent modifications should be saved on separate disk and clearly marked as an 
updated or modified version. JANRAD 98 as described by this thesis is version 4.0. Any 
modification hereafter should be clearly labeled version 4.1 or 5.0. A general explanation 
should be included as comments preceding the code in the janrad98.m file. Additionally, 
the author’s name should be added to the About JANRAD window found in the 
about_janrad.m file. If the modifications or additions were not part of a thesis, a detailed 
description of the theory and code should be included as an attachment to this thesis. Also 


a binder should be maintained with all attachments and theses related to JANRAD. 


To increase exposure to JANRAD 98 and to maximize user feedback, the 
program should be available for download on the MATLAB*® ftp site. Include Appendix 
A, JANRAD 98 Users Guide, as a Readme.txt file for those unfamiliar with its use. An 
email to select industry and university rotorcraft engineers announcing JANRAD 98’s 


posting on the site should also be considered. 
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APPENDIX A. JANRAD 98 USER’S GUIDE 


The JANRAD 98 Users Guide is written as a brief introduction to the Joint 
Army/Navy Rotorcraft Analysis and Design computer program. It is intended to explain 
the basic features and operation of the program and assumes a basic knowledge of 
helicopter mechanics and the use of the MATLAB® programming language by The 
MathWorks” Inc. 


A. SYSTEM REQUIREMENTS 


JANRAD 98 requires MATLAB® version 5.0 or MATLAB Student Edition 
version 5.0 or higher. It will not run on any previous versions. JANRAD 98 version 4.0 
will fit on a single 1.44 MB floppy disk and will need that much memory available for 
installation. JANRAD 98 requires only the hardware to support MATLAB? 5S. 


B. INSTALLATION 


The recommended installation of JANRAD 98 is accomplished by first creating a 
subdirectory of MATLAB called Janrad98. The entire contents of the JANRAD 98 


floppy disk should be copied into this directory. Include all M-files and .mat files. 
JANRAD 98 will not run without all of the .mat files. 


1 & MATLAB Path 


ee ee 
tie 
i 7 tovgt 


. i 4 5 g . La t, 4 c 

- > do ee c a r ~ pS AIS Sol 
Be OR nc Nits 4 a. Lin ra é g , q 
i f te 

M a] a * a 
atlab\janra ba: Ue a Ea mecomencentl uae 
bed s ~ a to Wy FR ‘ 
be “26 “a igs we a + 4 
eae te =. - 5: 
hog ie ee i J 7 D 


eaelee 






ee ee eae eae Ls Siar LM 

Pec \matlab\isnrad38 ees *! 

Z c:\matlab\gu examples ed 
4. c:\matlab\toolbox\matlab\datafun 
‘BJ c:\matlab\toolbox\matlab\datatypes 
a c:\matlab\toolbox\matlab\dde 
‘YA c:\matlab\toalbox\matlab\demos 
3 c:\matlab\toolbox\matlab \elfun 
#4 c:\matlab\toalbox\matlab\elmat 
4 c.\matlab\toolbox\matlabfunfun 
EB c:\matlab\toolbox\matlab\general 
ac \matlab\toolbox\matlab\graph2d 
@ c:\matlab\toolbox\matlab\graph3d 
#4 c: \matlab\toolbox\matlab\ graphics 
‘qc: \matlab\toolbox\matlab\iofun 

ZA c:\matlab\toolbox\matlab\lang 


%4 c:\matlab\toolbox\matlab\matfun Po ee Eee 
oe wd woe Habe the mile. mark an al abate pia + iacsdseiees ee EA 


GDOLDELEL_LOELLELOOODEEE LEE EOVELIE 








eates 


a: Ree 
Pe ea er, 


sue See ea ak 
SEAN NEEL ek 











fed) r60.mat 


ay 


ate 
Z 3 Nhe Sats 
SSS SEN NN AR 


2, 


o' ese 


im - wrastd 





Figure A.1. MATLAB 5 Path Window. 
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It is recommended to add this new subdirectory to the MATLAB 5 search path. 
This procedure will eliminate the need to change the working directory from the 
command line each time JANRAD 98 is run and allows you to work from a floppy disk if 
desired. Adding the subdirectory to the search path is accomplished by selecting File, Set 
Path... from the File menu. Change the current directory to the new Janrad98 
subdirectory by using the Browse button. Then press the Add to Path push button. You 
will then be given the option to save the new path or just use the new path for the active 
session. It is recommended to save this path. Figure A.1 shows the MATLAB® Path 
window. 


Cc. STARTING JANRAD 98 


Typing janrad98 (lowercase, one word) at the command line prompt of a current 
MATLAB session starts JANRAD 98. This action will launch the JANRAD 98 welcome 
window shown in Figure A.2. 
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Figure A.2. JANRAD 98 Start Up Window. 
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D. USING JANRAD 98 


As an example, the use of this program will be demonstrated by selecting a 
previously saved input data file and changing the weight, airspeed and pressure altitude. 
The input and output files will be saved and printed. After the performance analysis is 
complete, we will then iterate on airspeed from 80 to 100 knots in increments of 5 knots. 


First, from Figure A.2, select the Run/Edit Existing File radio button. It is usually 
easier to edit an existing file because Create New File will not give you the chance to 
change the working directory if desired. Once the selection has been made, press the 
Continue >> button. 
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Figure A.3. Selecting a File to Edit. 


The next figure window to appear is shown in Figure A.3. This window allows 
you to change the working directory and select an input data file. To change the 
directory, type or edit the desired path in the edit box. All of the .mat files listed in the 
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working directory are displayed in the list box. Input data files are saved as filename.mat. 
A note of caution here, each GUI window also has an associated guifilename.mat file. 
The user should name input/output data files using helicopter aircraft designations such 
as UH-60A, h-99 or h-design1 to differentiate from JANRAD 98 GUI files. Next, select 
an input file to edit by clicking on the file name. For this example, we will edit the UH- 
60A.mat file. Then press Continue >>. 


The Performance Input window will be displayed as in Figure A.4. The input 
data will be displayed within the appropnate edit boxes. Any or all of the parameters can 
be changed at this point without altering the original data file. You will have the 
opportunity to save the new data if you chose after the analysis has been completed. The 
airspeed, weight and pressure altitude can be edited by highlighting and typing 100, 
17,000 and 2,000 in the respective edit boxes. Pressing the enter key is not necessary to 
enter the new value. Using the Tab key or clicking on another edit box or control will 
enter any changes. The Print Screen button will print a draft copy of the GUI window 
with the displayed values. It however, will not record the file name for which the values 
are stored. Now press the Continue >> button. 
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Figure A.4. Input Parameters Window. 
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From the Input Parameters window, JANRAD 98 will call the Iteration Method 
window as shown in Figure A.5. JANRAD 98 will call its computational routines from 
this window. By choosing No Iteration and Analyze, JANRAD 98 will run the 
parameters selected from the previous window. The Analyze pushbutton initiates the 
computational routines. All controls on the GUI will be disabled except the Interrupt 
pushbutton. The Analysis Status Box will display the performance routine status and 
clock as JANRAD trims the rotor and adjusts the collective and cyclic mathematically. 
The Interrupt button will halt the routine and enable the Resume control and JANRAD 
Options menu on the GUI. This will allow the user to change parameters, quit or return 
to beginning. The Resume button will continue with the performance routine where it 
originally interrupted. It is worth noting that the Interrupt button will not always respond 
immediately. However, once MATLAB finishes its current line evaluation, the 
calculation will pause. 
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Figure A.5. Iteration Method / Analysis Window. 


The Status Box will inform the user when calculations are complete. The 
Performance Output window will be displayed automatically. The Performance Output 
window shown in Figure A.6 displays the performance results. These results can be 
saved and the screen printed from this window. However, it is recommended to print the 
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saved input and output files through the next window. By pressing the Options >> push 
button, the saved input and output files can be printed simultaneously and in a more 
usable format. The input/output files are saved after activating the checkboxes, typing a 
file name and pressing the Options >> push button. 
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Figure A.6. Performance Output Window. 


The Options window in Figure A.7 will provide the capability to print the latest 
files and select additional analysis routines. At this time however, the Stability and 
Control and Rotor Dynamics routines have not been completed. If selected, you will be 
reminded of this limitation. 


From here, we will go back to calculate performance parameters by varying 
airspeed. This is done by pressing the Change Iteration Method radio button and the 
Continue >> push button. This will go back to Figure A.5. To vary airspeed, press the 
Airspeed radio button and then Analyze. 
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Figure A.8 shows the Iteration Parameters window. By typing 80, 100 and 5 in 
the appropriate edit boxes and selecting Analyze, the performance routine will calculate 
and display plots of various performance results with respect to airspeed. Note: the 
Aspect Ratio edit box and HIGE check box are only enabled when Altitude Iteration or 
Wing Span Area is selected. Figure A.9 shows an example an airspeed iteration subplot. 
These plots are primarily used to examine trade off studies during the design process. 
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Figure A.7. Options Window. 
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Figure A.9. Example of Airspeed Iteration Output Plot. 
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The M-file janrad98_plots.m contains the details for the airspeed iteration plots. 
You can add, change or otherwise personalize the plots within this file as necessary. It is, 
however, recommended that any plots not desired be eliminated by commenting the 
appropriate lines and adding any additional code at the end of the file to preserve the 
original file structure. Always document changes made to the code for future users. 


E. HINTS FOR JANRAD 98 OPERATIONS 


The following paragraphs list some recommendations for the most efficient use 
of JANRAD 98. They are a guide based on observation, experience and knowledge of 
the code. Any other recommendations should be addressed to the Helicopter Design 
Instructor at the Naval. Postgraduate School for implementation into the next version of 
JANRAD 98. 


JANRAD 98 was designed for robust operations. However, because this is the 
first version to utilize a Graphical User Interface, not all cases of user inputs have been 
exercised. If the program appears not to be working properly, quit JANRAD using the 
options menu available on all but the first window. This action will quit JANRAD, close 
all MATLAB figure windows and clear the base workspace. Restart JANRAD by typing 
janrad98 at the command line. 


Using azimuth sectors greater than 24 and blade elements greater than 20 will 
dramatically increase computation time. The accuracy of the results does not improve a 
significant amount for these larger values. 


For more accurate results, use 0012, sc1094r8 or sc1095r8 airfoils when possible. 
These airfoils include Mach number inputs for ¢ and ca calculations. The HH-02 and 
VR-12 data do not depend on Mach number. — 


Using the print screen buttons on the performance input or performance output 


windows will take up to several minutes on older processors. Use the Send to Printer 
button on the JANRAD 98 Options page for faster and more compact printer output. 


A. 








APPENDIX B. JANRAD 3.1 INPUT/OUTPUT 


The following is the actual h60.prf? file from JANRAD 3.1 used to compare with 
JANRAD 98 output. 


toe .ReoU iS sa55 


h60 
Forward velocity = Li Ets 
Temperature = 66 degs F 
Pressure altitude = UEGOOEE 
Gross weight = 16500 lbs 
Number of blades = 4 
Rotor radius = -26.33-£t 
Blade mean chord = ee eek 
Blade twist = -18.00 degs 
Blade lift curve slope = Das 
Blade weight = 175.00 lbs 
Rotational velocity = 26.47 rads/sec 
Blade grip length = 300 LC 
Hinge offset = ao Le 
Equivalent flat plate area = 45.00 ft%*2 
Vertical projected area = 227.50 ft*2 
Wing area = O20 Sets 2 
Wing span = O02 00 "£E 
Wing CL = 0.00 
Wing CDo = 0.0000 
Wing CEELciency factor = O00 
Horizontal tail area = 45.00 ft%*2 
Horizontal tail span = 14.38 ft 
Horizontal tail CL = 0.80 
Horizontal tail CDo = 0.0100 
Vertical tail area = 32.30 ft%*2 
Vertical tail span = So 17 ie 
Vertical tail CL = 0.00 
Vertical tail CDo = 0.09100 
Fuselage drag = 1927 lbs 
Rotor drag = 336 lbs 
Wing lift = 0 lbs 
Wing drag = 0 lbs 
Horizontal tail lift = 1541 lbs 
Horizontal tail drag = 126 lbs 
Vertical tail side force = 0 lbs 
Vertical tail drag = 14 lbs 
Auxiliary thrust = 0 lbs 
Tip path angle = 9.06 degs 
Rotor coning angle = 4.16 degs 
Location of mean thrust (r/R) = 0.60 
Collective pitch at .7 r/R = 8.83 degs 
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lst lat cyclic term-Al (deg) 
ists rong cyclic term-Bl (deg) 
Solidity 

Disk loading 

Figure of Merit 

CT/sigma 

CQ/sigma 

CH/sigma 

Tip mach of the adv. blade 
Advance ratio 

Rotor thrust required (TPP) 
Rotor power required 


Rotor torque = 


1.46 
=4.25 
0.082 
0.00 lbs/f£t*2 
0.00 
0.071 
0.0049 
0.0016 
0.796 
0.270 
15147 lbs 
16) Tae 
27802 ft-lbs 
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APPENDIX C. JANRAD 98 INPUT/OUTPUT 


The following is the actual h60.prf file from JANRAD 98 used to compare with 
JANRAD 3.1 output. 


bt “RESULTS 2 a 


Forward velocity = iia ktes 
Temperature = 66 degs F 
Pressure altitude = LOO0> £t 
Gross weight = 16500 lbs 
Number of blades = 4 
Rotor radius = 26.83 ft 
Blade mean chord = eg Soe 
Blade twist = -18.00 degs 
Blade lift curve slope = S273 
Blade weight = 175.00 lbs 
Rotational velocity = 26.47 rads/sec 
Blade grip length = S200 Lt 
Hinge offset = Lazo Gt 
EQuUIVvalent: Elatiplate area = 45 005ft 2 
Vertical projected area = 227.50 ft*2 
Wing area = O200 fe 2 
Wing span = OF 00k Ec 
Wing CL = 0.00 
Wing CDo = 0.0000 
Wing efficiency factor = 0.00 
Horizontal tail area = 45.00 ft*2 
Horizontal tail span = 14.38 ft 
Horizontal tail CL = 0.80 
Horizontal tail CDo = 0.0100 
Vertical tail area = 32.30 ft%*2 
Vertical tail span = Sree en te 
Vertical tail CL = 0.00 
Vertical tail CDo = 0.0100 
Fuselage drag = W277 ibs 
Rotor drag = 336 ibs 
Wing lift = 0 lbs 
Wing drag = 0 lbs 
Horizontal tail lift = 1541 lbs 
Horizontal tail drag = 2G els 
Vertical tail side force = 0 lbs 
Vertical tail drag = 14 lbs 
Auxiliary thrust = 0 lbs 
Tip path angle = 9.06 degs 
Rotor coning angle = 4.16 degs 
Location of mean thrust (r/R) = 0.60 
Collective pitch at .7 r/R = 8.83 degs 


Shy) 


lst lat cyclic term-Al (deg) 
lst long cyclic term-Bl (deg) 
solidity 

Disk loading 

Figure of Merit 

CT/sigma 

CQ/sigma 

CH/sigma 

Tip mach of the adv. blade 
Advance ratio 

Rotor thrust required (TPP) 
Rotor power required 

Rotor torque 


1.46 
5 
0.082 
0.00 lbs/ft*2 
0.00 
Gncnt 
0.0049 
0.0016 
0.796 
0.270 
15147 lbs 
1338 h.p. 
27802 ft-lbs 
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APPENDIX D. JANRAD98.M 


This file launches JANRAD 98 and its Graphical User Interface. It is called by 
typing janrad98 at the MATLAB command line. 


function janrad98() 


% 

3 JANRAD98 .M 

3 Joint Army Navy Rotorcraft Analysis and Design 
3 (JANRAD) 

3 Version 4.0 

3 March 1998 

% Version 1.0 Designers 

3 MAJ Bob Nicholson 

% MAJ Walter Wirth 

% Version 2.0 Update Designer 
% LT Dale Feddersen 

3 Version 3.0 update Designer 
% LT Dave Eccles 

% Version 3.1 update Designer 
% LT Dan Hiatt 

% Version 4.0 update Designer 
% LCDR. Chris F. Lapacik 


% Version 4.0 added the Graphical User Interface. The GUI allows the 
user 

$ to create, run, save and print files with less effort and greater 
speed. 

% The basic performance calculation routines remain essentially the 
same as 

peEvctcolon 3. seuOwever, the Input, Output and file Struceure were 
modified 

% extensively. Sikorsky UH-60A airfoil data was also added. 


% Version 3.1 adds time varying tip loss and the corrected dynamics 
module. 

%$ The dynamics module provides Southwell plots and rotor blade 
response in 

% flap and lag motion. Also included is rotor flapping stability 
determination 

$ by Floquet analysis. 


4} 


Version 2.0 corrected minor bugs in ver 1.0 and incorporated 
Wheatley's Eqn. Additionally the user may now input a 
tapered rotor blade. Finally a menu was created for 

ease of performing various iterations and then saving that 
calculated data for later use/manipulation. 


0 &® ® © OP 


This program is an interactive preliminary design tool 
developed to aid the design student in determination of 
initial rotorcraft configurations and in the calculation 

of performance, stability and control, and other parameters. 
The program will work for conventional or compound rotorcraft. 
It will provide accurate data for airspeeds less than 10 
knots and greater than or equal to 50 knots. 


cP WH PP ®P OP CP Oo 


load janrad98 
global H JAN H_EREF H_CNF 


H JAN = figure('Units', 'normalized', 
eeolLor’, 10-0 0.8 0.8], 
"Colormap',matoO, 
"MenuBar','none', 
"Name', 'JANRAD 98', 
“NumberTitle*,*off",: « 
'PointerShapeCData',matl, ‘. 
Posi cion’, [-0.003125 0.0625 0.954683 0.8625], 
He US a plane stale Ue 
b = uicontrol('Parent',H JAN, 
"‘Units', 'normalized', 
"BackgroundColor, (0. 752941 0.752941 0.752941], 
'FontSize',16, 
'FontWeight', 'bold', 
'Position', [0.0604396 0.757143 0.434066 0.145714], 
'String','Welcome to JANRAD 98', 
'Style', 'text', 
“Tag', ' StaticTexti" ) ; 
Db SsuUrcenermor( Parent. ;H JAN, 
PUniTeSe normallzedigs =. . 
Petey erounacolc:: PaO. 752941 O. 752941) Ca e294 Li 
'BOntsi zean10, 
aromtwei one’) “pold’. 
"Pesitionwm|0.07585097"0.617284 0.40322680. 0925926], 
'String','Joint Army/Navy Rotor Analysis and Design', 
“Sevier, text', ; 
lage, “Staticlext?*) ; 
De] Ui conGrom Parent: ,H JAN, 
nies - WOrmallZzea' ;, .. 
Rae amenccolous lOc 752941 0247 529416.0).752941), 
'Position', [0.0766129 0.395062 0.403226 0.0617284], 
'String', 'Department of Aeronautics & Astronautics’ 
“POntolze 7,10, 
PEOnEWerant ,' bold”, 
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(Styvle' texte, . 
'Tag' *StaticText6' ae 
b = ulcontrol( Farents HlJAN, 
'Units', 'normalized', 
‘packgroumacoler Pdi iif Oe 752941 0.752941 0. 75294147 
'Position', [(0.0766129 0.333333 0.403226"0 061 /2can7 
'String', 'Naval Postgraduate School', 
"PONCE S1.Ze 2eu0- 
PPontwel ght asco lan, 
“Style! ~ textes : 
'Tag', 'StaticText5'); 
Db =suicontrol('Parent~,HAJAN, 
bUniCS )s nermalized . .. 
ppackgrounacenon: 7 LO. 752941 0.752941 0.752941], 
"Position', [0.0766129 0.271605 0.403226 0.0617284), 
'String', 'Monterey,CA 93940', 
yeontSize* a0; 
"FontWeight', 'bold', 
est vile’ *text", : 
'Tag', 'StaticText4'); 
b = uicontrol('Parent',H JAN, 
[Unaes! a nermallzea", 
‘BackGroeundColor*,102.7592941 0.7529410.752541), 
"Position ’, (0. 0766129 0.17284 0.403226 2040771605), 
EStLrings 7 March 2o,c 79s 
'Rontsize 710, 
"FontWeight', 'bold', 
SScyle* >  Cext 4.0... 
nlaqd:~@ otaticrexes”*); 
H EREF = OEE Oe Eee ee 
‘UnLtes » normalized : 
feaipsor Sivlameaeo Lemme ret jamie. 
TPOSPeE1On (Oro c7 792 00'.699172 0. 275362 0.103448], 
PSeeiniga,. bale, RUneExXisting File * 
'Style', 'radiobutton', 
"Tag’, Radiobuttonl’, 
pvarue: i); 
H CNF = uicontrol (‘'Parent’,H JAN, 
'Units', 'normalized', ; 
Peal lback”™ 7) janradgs cn chi”, . 
pPOsSLeton, s102.007992 UV.5l4d368 0. 275362 02106322) 
"String', 'Create New File' 
"Style’, ‘radiobutton’, 
‘Tag’, "Radiobutton2 '); 
be UlCOntroh( Parent , i JAN, 
‘Units', 'normalized', i 
"Callback*’, ‘janrad98 fcn stop', 
PEOMeoL Ze. 12, 
-EOneWelgne ; DOLG. 75.2: 
‘Position’ 700-567 992 07295977 0. 275362 0.112069), 
rStLring’, Ourt DJANRAD 9a"; 
"Tag', 'Pushbuttonl'); 
b = uicontrol('Parent’,H JAN, 
‘Unies 7 toermalized’, 
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"Caliback', 'janrad98 fcn cont’, 
PFontSize' ,12, 
PNoncewelght’, "bola’, .«.. 
boston .10- 590062 802149425 0.273292 0.112069), 
‘Strang’, * Continue >2a5 
bbag | PushbuttenZ 7 

b = uicontrol('Parent’",H JAN, 
‘Units’, normalized’,... 
"BackgroundColor, [0-. 752941 05752941 0.752941), 
Pyomtolze’, 12, 
'FontWeight','demi', 
"Post-taon’ > 10.591097 0.793103 0.269151 0.091954j- 
'String', 'Select Option and Continue', 
wpeyac , text’, 
‘Taq’, sStatielrextZ.); 

be Unecontrol( Frarent {7h oJAN, 
‘Units * Pe cmasced, 
Hae aecundcolior' Pere 752941 O77 52941 0. 752941], 
BeOSitLon , 70eoolL0/7 "0.445402 0.329193. 0.482759], 
"Style', 'frame', 
‘Tag', 'Framel'); 

b = uicontrol('Parent',H JAN, 
‘Units', 'normalized', 
"BackgroundColor', [0. 752941 0.752941 0.752941), 
PEOGTtLonewey0.0S842742 0.570988 0.479839 0.367284], 
"Style', 'frame', 
'Tag', 'Frame2'); 

be -sUteontrol, Parent =, JAN, 
PUT E Ss, eridlizean 
ebackqrounacolor ., (0. 752941 0. 752945680. 752941), 
PPositvonm , foe@ea 7596 0.149758 0.417349 0.345411), 
‘Style’, 'frame', 
'Tag', 'Frame3'); 


ascigmhint base | WebREF 7H EREFE) ; 
assignin('base', 'H CNF',H CNF); 
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APPENDIX E. JANRAD98_FCN.M 


Switchyard Callback function called by the janrad98.m GUI function. 


function janrad_fcn (Action) 


%$ Switchyard Callback function for janrad98.m 
% JANRAD 98 VERSION 4.0 


global H_EREF H CNF H_JAN NAME COUNT 


condl=get (H_ EREF, 'Value'); 
Sonac—det (HaECNE Value —)s, 


COUNT=0; 


we nargin, 
switch Action 
case ‘'cont', 
if condl== 
analysis 
close (H_JAN) 
else, 
Name it 
performance input 
close (H_JAN) 
end 
case ‘stop' 
efelay e fefbla 
case ‘eref' 
Set (heehure .alue 1) 
set (H_CNF, 'Value', 0) 
case ‘cnf' 
See (Heebner, Value"7.0} 
Set (HIeChib,, Value oi) 
end 
end 
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APPENDIX F. ANALYSIS.M 


This file creates the GUI to select a JANRAD 98 data file and select the analysis 


method. It is called in the janrad98_fcn.m Switchyard Callback function. 


function analysis () 


%$ GUI figure window to Select Type Analysis 
% JANRAD 98 VERSION 4.0 


%$ This is the machine-generated representation of a Handle Graphics 
object 

% and its children. Note that handle values may change when these 
objects 

%$ are re-created. This may cause problems with any callbacks written to 
%$ depend on the value of the handle at the time the object was saved. 

% 

%$ To reopen this object, just type the name of the M-file at the MATLAB 
%$ prompt. The M-file and its associated MAT-file must be on your path. 


load analysis 
global H_P H_SAC H_RD H ANAL H_ LB NAME 


H ANAL = figure('Units', 'normalized', 
PCOlor (0.8 0 670.0); 
'Colormap',mat0O, 

'Name', 'Analysis', 

‘Numperiutle”® ,*off*, 

"PointerShapeCData',matl, ... 
Prosition' , (-O2005125.0.0625 02954688 0.39625), 
hag, PiG2 ja; 

be uimenu ("Parent *,;H ANAL, 
"Label','JANRAD Options', 

'Tag', ‘uimenul'); 

c = uimenu('Parent',b, 
‘Callback','analysis fcn quit’, 
"Label', 'Quit JANRAD', 

'Tag', 'JANRAD OptionsSubuimenul') ; 

c = uimenu('Parent',b, 
(Callback jana yors Len return 7... 
'Label', ‘Return to Begining', 

'Tag', 'JANRAD OptionsSubuimenul') ; 

c = uimenu('Parent',b, 
‘Callback',*analysis {cn delta input",... 
‘Label', 'Change Input Parameters', 
‘Tag’, subuimenul }- 

c = uimenu('Parent',b, 
Caliback 7 sanalyors €¢n about’ ,... 
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Whabel., “Aboout Janradso0e...-', 
'Separator','on', 
'Tag', 'Subuimenul') ; 
b = uicontrol(‘*Parent’,H ANAL, 
'Units', 'normalized', 
eee roundeolo:' vai [.Oxc 752941 Oe 94 Le Ol OZOA 
"FontWeight', 'bold', 
'Position', [0.106033 0.896175 0.332724 0.0546448], 
"String','Type Path to Working Directory', 
"Style','text'’, : 
“TNag*; 'StaticText2"); 
H WORK = uicontrol('Parent',H ANAL, 
‘Units', ‘normalized’, 
Hee roroundcolor: nie ll, 
"Position', [0.104205 0.822404 0.330896 0.0546448], 
bereatercn ’ ;" a...» 
‘Stang * = 
"HorizontalAlignment', ‘left', 
Style) “edit, 


sGaliback |. ca (get (H WORK,” strings) ));7list=din( 9 ~.matee );str—{listene 
ie i aie et hays Vinee whe ia) 
blag’ , -caiclextl )-; 
De-SULecCnerol PabrenG.,H ANAL, 
‘Units', 'normalized', 
"BackgroundColor',[0. 752941 03752941..0.752941), 
PHOneES1ze*, 12, 
"FontWeight', 'bold', 
'Position', [0.102377 0.734973 0.340037 0.0628415], 
"String','Select Data File', 
wotyle’, *text', Bee 
‘ag’, staticlext3"); 
H_ LB = uicontrol('Parent',H_ANAL, 
'Units', 'normalized', 
ree ouneolor: Pa) (i ae il, 
TROS@ttion , 10-1017 cG 0.355742 0.342857 0.369746), 
SSL LNG oy, 


'Createren’ , list—adir("’*.mat'*);str={list.name} ;set (qcbo, , str Ast re 
cae: 
TistesrngSt, Mean 
"Callback', 'global 
NAME; Value=get (gcbo, ''Value'') ;String=get (gcbo, ''String'') ;NAME=String{ 
Value};' : 
Lad, ULStbOxl); 
"Value',1); 
b = uicontrol(’Parent',H ANAL, 
‘Unwes’> “normalized', 
\BackqroundColor’, (0. 752941 0.752941 0.752941], 
PRontsize’ p12, 
"FontWeight', 'bold', 
PEGsition, (0.098214 O-815I26 0.285714 0.0644258)], 
"String', 'Select Type Analysis' 
rStvyLe’, text’, : 
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'Tag', ‘StaticTexe®) - 
HP = uilcontrol( Parere im aiAL, 
'Units' ‘normalized 
BackgroundColor’ pO 752941 Oo Z9 4 Or o2941 iF 
‘Callback 7) analliyeis etcn ames, 
'FONESTZe fel Z eae 
"Position! (0. 598214" 02672269) 0.2657 140052437 ie 
'String', ‘Performance’, 
Style", radrobueton: , 
'Tag', 'Radiobuttonl', 
"Value',1); 
H SAC = uicontrol('Parent',H ANAL, 
‘Units’, ‘normalized. , 
"BackgroundColor POs 752941 O-7 52940. 752941 )\- 
bCaliback 7 sanalyois Tchemssae., 
PhOmts1.2Z6°%.L27 oe 
PPOSTELON +, |O. 590214 90. 5406016 0.285714 02092437 )]- 
PScraAng-, stabil vey Ana Control’, 
Lovie 7. Ladi! oouceon |, 
'Tag', 'Radiobuttonl1'); 
H RD = uieontrol (*Parent’,H ANAL, 
'Units', 'normalized', 
"BackgroundColor nLOs 752941 0.752941 0.752941], 
‘Callback", ‘analysis fen h_rd', 
“Hontsize', 12, 
\PoOsteron , (O.598214 0.408964 0.285714 0.0952381), 
VSELING! > ROLOE DY namics , 
vstyle!> tradiobuteton ,; 
nag ;) Radlobuctonl") > 
b = uicontrol('Parent',H ANAL, 
‘Unwes! se normalized’? : 
"Callback', ‘analysis fen back', 
"FOntsize , 12; 
phontWeirghnt. ,.bo1d", 
PeOsttvonw es 0 20S0zcl43s 0.176 0.210714 0.096), 
‘String 7-<< Back, 
‘lag 7 PUshbutEOon! ') > 
Pa Leone nol Eanrent | H ANAL, 
‘Units', 'normalized', 
peal pach analysis Chr enx 7)... 
PFOnuSiEZe. 7 UZ, 
rFOntwelgnt | “bold, 
"POSItion'; (0.391071 0.176 0.208929 0.096], 
foeling., Cancel”, : 
Plage, eusibuctonl }); 
b= uilcontrol (*Parent',H ANAL, 
‘Units', 'normalized', : 
'Callback', 'global COUNT;COUNT=0;analysis fen cont', 
‘FORESPZe., 2, 
'FOonkgWeignie 7 botd', 
"Position', [0.678571 0.178667 0.208929 0.096], 
 SerInge Continue >>", 
‘lag’; sPusnbuttonl") ; 
b= uicontrol ('Parent’,H ANAL, 
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UMhes, ano rma bi Zed! a... 

Peete hObounccolor ,(90.752941 0.752941 Oms2z941), 
Beosmedom- seem 143 .0.35294) 0.355357 0.557423), 
"Style’, 'frame', 

‘Rag’, *"Framel*); 


AS SEGtenn soase yon Po) HP); 
assignin('"base', 'H_ SAC',H SAC); 
assignin('base', 'H_RD',H_RD); 
assignin('base', 'H WORK',H_WORK) ; 
assign ase 7 no Lb | HO LB) ; 
assignin('caller', 'NAME', NAME) 
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APPENDIX G. ANALYSIS FCN.M 


Switchyard Callback function for the analysis.m GUI function. 


iunetion analysrs fon (Action) 


$ Switchyard Callback function for analysis.m 
% JANRAD 98 VERSION 4.0 


global H P H SAC H RD H ANAL 


eondl—get (HoUP,, Value}; 
cond2=get (H_SAC, 'Value'); 
eonas—get (HRD, Value}; 


if nargin, 
switch Action 
ease hap: 
Set (Hie. value 7) 
set (H_ SAC, 'Value',0) 
Set Herne value 70} 
Caco es hesdc” 
set(Her,. Value s.0} 
set (H_SAC, 'Value',1) 
set (HERD, Value), 0), 
Case h Ed’ 
set(iee, Value 70) 
set (H_SAC, 'Value',0) 
set (HeRD, value, 1} 
case ‘cont' 
if condl== 
Peefomance sin rur 
close (H_ANAL) 
elseif cond2== 
Sbabt Prey sdnccont ro. 
elseif cond3== 
rotor dynamics 
else, 
error('SomeThing is Wrong in Analysis Function') 
end 
case ‘cnx' 
analysis 
close (gcf) 
case 'back' 
janrad98 
close (H_ANAL) 
case ‘return' 
janrad98 
close all 
case 'quit' 
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qu ist quia: 

Sone ae hE ae put : 
performance_input 
close (H_ANAL) 

case ‘about’ 
about _janrad 

end 
end 
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APPENDIX H. PERFORMANCE INPUT.M 


This file creates the GUI that displays the 36 input parameters as loaded from a 
previously saved file or created new by user. It is called in the analysis fcn.m 


Switchyard Callback function. 


Lunction periEormance anput() 


% GUI window to display and/or edit input values. 
% JANRAD 98 VERSION 4.0 


%$ This is the machine-generated representation of a Handle Graphics 
object 


%$ and its children. Note that handle values may change when these 
objects 
%$ are re-created. This may cause problems with any callbacks written to 


depend on the value of the handle at the time the object was saved. 


To reopen this object, just type the name of the M-file at the MATLAB 
prompt. The M-file and its associated MAT-file must be on your path. 


we ® ® oP 


load performance input 
global COUNT NAME H_PERF_IN S_USER_INPUT S PERF_INPUT 


switch COUNT 
case QO 
if ~isempty (NAME) 
eval([{'load ',NAME] ) 
unstructurel 
else 
load create new 
structure 
end 
case l 
unstructurel 
end 


PeEERE IN = figure( Units", normalized", 
Beaton’, (OsG-0..6 Oss], 
"CreateFcn','', 
"Colormap',matdO, 
"Name', ‘Performance Input Parameters', 
‘NumberTitie), "orr 
"PointerShapeCData',matl, 
"Position -t-0.003125 0.0625 0-954688 0.86251, 
Nag" 77 ralg2  )-s 
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heeopes -muimenud Searent 7H PERE IN, 

"Label', 'JANRAD Options' 

'Tag', ‘uimenul'); 
Ce—euamenu( Parent 7h Opt, 
‘Callback, “performance input fen quit’, 
'Label', 'Quit JANRAD', ; 
'Tag', 'JANRAD OptionsSubuimenul1') ; 
uimenu(" Parent',h opt, 
"Callback", ‘performance input fen return’,... 
"Label', 'Return to Begining', ... 
'Tag', 'JANRAD OptionsSubuimenul'); 
uimenu('Parent',h opt, 
"Callback', 'performance input_fen delta_input', 
"Label', "Change Input Parameters' 
"Enable',‘'off', 
Tag td Souneimenuii iy: 
uimenu('Parent',h opt, 
'Callback', 'performance input _fcn about’, 
'Label', 'About Janrad 98 ...', 
"Separator','on', 
“lag ,usibuimenul") ; 


‘o) 
ll 


0 
It 


Q 
tI 


OF, 
tl 


UUConE rol sParent HereRe iN, 

‘Units', 'normalized', 

ec oeendcoler PLO 752941 029752941 0.792941). °-... 
Soom et olen 0702 12700 0.92029 0.175123 050493092), 
"String', 'Pressure Altitude (ft)', 

(Styvlet,' text’, : 

"Tag, ‘StaticText2*); 

uicontrol (*Parent’,H PERF IN, 

"Units', 'normalized', 

Mee charoundcolers lel ak 1], 

"Postteron 710.217670 0.922705 0.0981997 0.0483092) 7 
ieee eae 

PoOeLang ; PA 


2. 
lI 


‘Callback’, 'PA=get (gcbo, ' String’');S_USER_INPUT. PA=str2num(PA) ; ° 
'Tag', 'EditText1") ; 

dvarlneontroil(” Parent’ ,H PERF IN, 
"Units, ‘normalized; <<. 
"Packgrounacolor™, (0. 152941 O77 5294 On 1 52941), 
Weosit1on —7)0. 022760 0.059072 0.175123 0. 0433092), 
'String', 'Temperature (deg F)', 
“elev, “CExt,”, : 
(Lage, Stabiclrext2"); 

ae sutconit rol Parent’ 7H PERE IN, 
‘Units', 'normalized', 
Macho Lorna licx Peles VE ll, 
*Position', [0.217676 0.857488 0.0981997 0.0483092], 
PoE yle ww eart.”; 
CString +, emp, 25 
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*"Callback', 'temp=get(gcbo, ''String'');S_USER_INPUT.temp=str2num(temp) ;' 


t ® e ® 


‘Tag’, ’EditTextl ae, 


d = uicontrol(’Parent 7,HeEBRFSIN, 
‘Uniwes"s, ‘noma lee edae : 
"BackgroundGoteor |) 0: 752941 0.752941 0.752941], 
"Position’, (0.0212766 0.792271 0.1/512350.,0402092Z)— 
‘String’, ‘Airspeed (kts)'’, 
‘(Style text |, : 
'Tag', ‘Staticlext2 )-; 

d = uicontrol("Parent*,H PERF IN, 


"Units *, a 

'Backgroundcolor’, [1.1 ll, ee 

‘Position, (0.217676n00 792271 0.0981997 0.0483092)], 
rouyle', "ede; 

VOU LLG yy, Vint 7 


'Callback', *Vinf=get(qcbo, *’String'');S USER INPUT. Vinfi—=steZnum(Vint) ;* 


'Tag', ‘EditTextl'); 


ee ULcoOntrol( Parent 7oeEERE TIN, 
tUnwes*) normal rzeda 
"BackgroundColor, [0. 752941 0.752941) 0.752941), ae 
'Position’, [0.0212766 0.727053 0.175123 0.0488092]; 
VSELEING , Gross Wt... {ibs.) 7 
notyle', '"text’, 
‘lag, Staci clexeZ.), 

d = uicontrol ('Parent' ,H_PERF_IN, 


BUrL cS « eal cod" 

"BackgroundColor',[{1 1 1}, sae 

"Position', [0.217676 0.729469 0.0981997 0.0483092], 
fory le’ vedi tir 

‘Strang: , GW, 


‘Callback', 'GW=get (gcbo, ''String'');S USER INPUT.GW=str2num (GW) ;' 
"Tag', 'EditTextl'); 

gee CONCrOl(’ Parent. ,HePERr IN, 
Ponaigs ‘nome. 
"BackgroundColor',[0. 752941 0.752941 0.752941], ; 
*"Position’,(0.0212766 0.661836 0.175123 0.04838092)¢ 
'String', "Rotor Vel. (rad/sec)', 
‘Style’, 'text', , 
‘Tag’, ‘StaticText2'); 

geet concrol( Parent.) HOPERE IN, 
‘Units’, “norma Lien 
"BackgroundColor',[1 1 1], 
"Position', [0.217676 0.664251 0.0981997 0.0483092], 
"Style* "edit", 
‘String ,omega, . 


"Callback', 'omega=get (gcbo, ''String'');S USER_INPUT.omega=str2num(omega 


1 ps” aoe 
‘Tag’, ‘EditTextl1'); 
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om 
Il 


Uteontrol( Parent yi PERF IN, 
nits >. normalized’, 

(BackgroundColor 7 (02702941 0.752941 0.752941), 
"Position',{0.0212766 0.596618 0.173486 0.0483092], 
‘St&ing , No. Azmith Sectors", 

'Style','text', 3 

tag’, Stati clexuZ. 

UTeoOncCrol Parene. 7; HoOPERF IN, 

EOh owe eta). ee cea 

rae eanoundeoloxe Pal We ll, oye 

‘Position’, [0.2176/76 0.599034 0.0981997 0.0483092], 
YStyle'; edit’, 

 SEing 7, NaZ jr. >). 


QO. 
Il 


"Callback', 'naz=get (gcbo, ''String'');S USER_INPUT.naz=str2num(naz) ; ' 


"Vag > Eaverextl"); 

Uleonerol( Parent Ho RERF .IN, 

‘Units', 'normalized', 

Pe ecundcolor: lle 752941 Of 752941 0.752941); 
PRos@ereroen 5 (0.0212766 0.533816 0.175123 0.0483092)], 
TSerangu.«COollePatch @ .7>£/R', 

PStyle | text; ne 

‘Tag’ 75 Seatlerext2-)c, 

UEcCOnerEOl( Parent ~HoPERE IN; 

‘Units', 'normalized', 

mee oraundc olor: Pe i lel ll, ; 

‘Position', [0.217676 0.536232 0.0981997 0.0483092], 
'Style','edit', 

"Sering ;, thetao,.6. 


©. 
II 


Q. 
Il 


‘Callback', ‘thetao=get (gcbo, "*String’’);S USER INPUT.thetao=str2num(the 
tao);' ; 

wag , EditText is); 

Umeent Loli Parent 7H PERF IN; 

‘Units' ea Sa" | = 

meckoroundcolsn UR 752941 0.752941,.0. 7529417 
"Pos?tion’=10.0212766 0.468599 0.173486 0-0483092], 
rstring’, ‘Wing Area (ft*2)', 

‘style’, text; eh 

"Tad, StaticlextZ |, 

WreOneLOl Parene 7 HePERE oIN, 

"Units', 'normalized', 

Mpaenecoundcolor' paleo il, 

"Pos@trvon:, (0.217676 0.471014 0.0981997 6.04383092\, 
Tetyle'; edit'; 

VSotuind: , “Wing, 4). 


d 


O. 
Il 


'Callback', 'Swing=get (gcbo, ''String'');s USER_INPUT.Swing=str2num (Swing 
ee 
plagi) Edi erextl” }; 
AS sUteonerol( (Parent 7h PERE IN, 

"Units’, ‘normalized’, .. 

eoeteroundcol on: LO: 752941 0.752941 0.752941], 


56 


'Position', (0.02966 0.403382 O.173486 0. 0483092) 
‘String |. Windgwopan( £e)oo, 
votyle @ textes, : 
‘Tag 7.. > CablelexeZz aie, 
qd = uicontrol( Parent 7 Heri Re ein, 
‘Units' ome Le 
‘peckqrounde aioe: Pon ad Gat 1], 
"Position. (0.217676. 0.4057 9770,0968199770, 04830972], 
fst ylel, ede, 
Pot ring” ,;Owemegs... 


'Callback', 'bwing=get (gcbo, ''String'');S USER INPUT. bwing=str2num(bwing 
as A 
‘rag 7 BALETeKt!  ); 
> control Parene ier E he _IN, 
Urmaites eeeraiied 
Pana eoundcolont iO . 752941 0.752941 0.732941], 
TPOsteron (050212766 0.34058 0.173486 0. 04630S21~ 
'String', 'Expected Wing CL', 
'Style', 'text', , 
‘hag: , "StacicrextZ* )] 
ad = uicontrol (* Parent ';HeeERF IN, 
Units’, "normalized ', 
specrarouncolon er (dies 1], 
“Position*, [0.217676 0.342995 0.0981997 0.048309Z2/ie 
'Style','edit', 
sString , Ciwing,< 


'Callback', 'CLwing=get (gcbo,''String'');S_USER_INPUT.CLwing=str2num(CLw 
TAG) paces ese 
tag’, ‘EditTextl'); 
Bgoeul control ( Parent, H PERE IN, 
FUntes* Seo rmaie ed. 
seer roundeolonnn a: 752941 0279294170. /5294ale 
Position’, (0.0212766 0.275362 0.173486802049 3092)" 
PerIiga, Wind CLO, 
Sot vee. text", ‘ 
'Tag' | 'StaticText2' )3 
fea elecont Lolebarent' > H PERF IN, 
LUnits ‘ptepaien nee, 
ee ckGrontleonen 7 (ie 1). : 
Phosition ,ren,ziq¢go76 0.277778 0.0981997 G.04b3092]7 
Poeyvle. . edit’, 
Poering 4 CDowing,... 


‘Callback’, 'CDowing=get (gcbo, ''String'');S_USER_INPUT.CDowing=str2num(C 
Dowing);', : 

Tag’, BALeELexti '); 
d = uicontrol('Parent' ,H_PERF IN, 

Wiaieess: aie cd 

po inners: Pa Om 752941 Oe 7oZ94 10.7 5294S 

"Position', [0.0212766 0.210145 0.173486 0.0483092], 

‘String, Wing Eff. Factor =e! 

 oeyle: 7 ext, 
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JhAGg) yiotaticlext2 jy 
d = wuicontrol(*Parent',H PERF _IN, 
BUT ts *s, Pe eralicedt, 
"BackgroundColor', [1 1 1], ee 
EPostemon ,.60.21/7676 0.21256 0.0981997 0.0483092], 
tseyv ile", "edit*, 
"Seling , wing, .. 


'Callback', 'ewing=get (gcbo, ''String'');S_USER_INPUT.ewing=str2num(ewing 
dit pees 
‘Tag’, EditTexti"); 
Ul comtrolmy Parent *,H PERF IN, 
‘Units*;7 normalized", 
"DacKkgrounaGolor’, (0. 752941 OF 32941 Oe S294 sy 
"Position', (Ove5e429 0.922705 0.173436 °0.0483092) , 
'String', 'Blade Airfoil Type', 
Suv le* |; Cexte, : 
'Tag', 'StaticText2'); 
Pepe] Uteontro | Parent. 7H PERE =N, 
‘Units', 'normalized', 
PBackacounacolor ’ , (0. 752941 0. J/az2941 50. 752041), 
PLest elon 0.249018 0.922705 0.12556, 0483002I, 
PoE cIng: 70012 hH-O2 | VR-l2 | sClOo4r3)SCl09ars , 
‘Style’, 'popupmenu', 
'Tag', 'PopupMenul', 
'Value',afoil, 
(ealbsac , deer r—get(gcbo, "'Value*');S USER INPUT. afoil—-afoil;); 


d 


qd — sui cone rom( Parent |, PERE IN, 
SUnats |; ernaliseas, 
“BackgroundColor, [0- 152941 0.752941 0.752941), 
“Position » [0.358429 0.855072 0.173486 0.0483092), 
Pstring’, Blade@eaft Curve Slope’, 
PoeVie’ pttexe’s a... 
‘Tag';) Statiwerextz ' ) ; 

dy Uicontror@ Parente PERE IN, 
'Units', 'normalized', 
"Backgroundcolor', [1-1 1}, . 
‘Position’, [0.549918 0.857488 0.0981997 0.0483092], 
'Style', ‘edit’, 
Bie eG sy clpctes 
"Callback", 'a=get(gqcbo, '*String'');S USER INPUT. a=strZ2num(a) ;" 
‘Taq’, Ediatelext li; 

a>= uaconerol Parent ,H PERE IN, 
“Unees*. “normalized’, <«. 
BackgroundColor’, (0. 752941 Omg 52 94150275204 1). 
EPOSteton 7 10-550429 0.794686 0.175123 0.0458937)], 
"Soeriang , No. ~ Blades’, 
Seyle ,. texc", : 
lag, sctatlerexte ); 

Of— SU reont rol Parent tere eh, 


"Units *> eerie 

"BackgroundColor',[1l 1 1], 

"Position', [0.549918 0.792271 0.0981997 0.0483092]), 
SeyLe > Calt -, 
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‘Secing ,b, : 
‘Callback', 'b=get(gcbo, ''String'');S USER _INPUT.b=str2num(b) ;' 
"Tag", Barelextl).,; 
qd = uicontrol (2 barent. ,HUPERE SIN, 
"Units 3 eer cean 
"BackgroundColor', [0. 752941 0.752941 0.752941], 
'Position', [0.358429 0. 731884071 75123 "OM04SS0S2 7 
‘String , BladewRadius ire.) , 
"Style: 7 veext.. : 
‘Tag, Statrerexe2)).; 
@ = Ulcontrol(eearenmta, 1 PEREGIN, 
'Units', 'normalized', 
‘Rach yroundeemens | (ee 1}, 
'Position', [0.549918 0.729469 0.0981997 0.0483092], 
'Style','edit', 
(SEEING pj Ryacee 
'Callback', 'R=get (gcbo, ''String'');S_USER_INPUT.R=str2num(R);' 
'Tag', 'EditTextl1'); 
ee  OULtCenerom( Parent ,H PERF IN, 
MUnt ts > Seer ie 
pBbackgroundcolor _ (0. 752941 0.752941 0.752941), 
‘Position', [0.358429 0.664251 0.175123 0.0483092), 
Boering: , Hinge Ottset (te.)", 
pocvile*, *texe-| 
'Tag' Rcramiererts ie 
Gee Ul concerol ("Parent ,H PERP IN, 
'Units', 'normalized', ; 
“BackgroundGolor'AbL 1 1), 
"Position', [0.549918 0.664251 0.0981997 0.0483092], 
Pstyvles “edit, 
i eed EG) Gp oi 
Callback’, 'e=get (gcho, '"String’');S_ USER INPUT.e=str2numi(e) 7" 
tag, EBditText1'); 
Wee cOonterol(" Panent., 7 EERE IN, 
‘Units', 'normalized', 
aeacioroundenton: LO. 752941 0.752941 0.752941], 
"Position', (0.358429 0.601449 0.1751237070463092)- 
"String',Non-Aero Part (ft.)', 
rouyle", ‘text’; : 
plag-| staticlextZ ji; 
ge Uicontrol (~Parnent ,H PERE IN, 
‘Units', 'normalized', 
SBackgrouncdeolor ;, (1 1 1], 
"Position', [0.549918 0.599034 0.0981997 0.0483092], 
Poty let Sedi ty 
Sel Nng 7 OE 1. 


Callback", “grip—-gerigebo, String’ )7S USER INPUT Giep-StrZnum(grrp) 7 
‘Tag’, Eai tiext li) 

d = uicontrol('Parent' ,H_PERF_IN, 
“Unies” Wioeta lived. 
Baa ROEOUnEeeO r PEO. 752941 0.752941 0.752941], 
"Position', [0.358429 0.538647 0.175123 0.0483092], 


ao 


SSeringu. Blades Root Cha (ft.)', 
Poy leur tex tis : 
whag =, ~otatichext2 ‘); 
de -sUtecherol Parent", H PERE <IN, 
‘Units', 'normalized', : 
[Backgqroundcolor', (1 1 1], 
[pose nom yO, 54991800, 5536232 0.0981997 0.0483092], 
UStyile'|, edit’, 
Poerang ,rchord, 


‘Callback, renord=get (gebo,* "String'');S USERSINPUT: rchord=st@e2num(rch 
CEG); 


lage, EG1erext i”); 


d= waseontrol({ Parenti PERE IN, 
‘hits > notmalLezed- 
Tere roandealon 7 O:. 752941 0.752941 50. 752941), 
'Position', (0.358429 0.471014 0.173486 0.0483092], 
"String', 'Blade — Ratio’, 
YOeVile*;, “EeExt ; 
‘Tag' Beat eree De 

d= Wieconurol( Parent:, HeeeRe. IN, 


Unaes.. normalized: , 

‘BackgroundColor, ft) 1 1}, ais 

‘Position', (0.549918 0.471014 0.0981997 0.0483092], 
"Seyile*, edit’, 

wotrIngd, » ti, as 


pea Pipa acre (ch qassrr ing: 9,5 UsBRIINPUT tCr=stEeznum( tr)’ 
‘Lag, EdL eText laa 


ae—suteconerol( Parent. ,HoPERELIN, 
‘Unites. normalized‘. 
pacKGhoundeolon: ZALOe 752941 Oye 94170.752941), 
UeGstironeso.coo429 0.405797 0.175123 020483092), 
pString , taper Starts @ (r/R)", 
‘seve , text, : 
“Tag, Staticlext2Z*); 

Qa Uteontrom ne Parent ji TERE iN, 


‘Units', 'normalized', 

eer oundeolor’ malting: 1}, 2 

"Position', [0.549918 0.405797 0.0981997 0.0483092], 
PSotbvyilre*, “edi6 ye. 

StkEing ,trse, 


‘Callback’, 'trst=get(gcbo, ''String'');S_USER_INPUT.trst=str2num(trst) ;' 
Pe 


TLaguen caLcrextl"); 


do = Wreouerol(* Parent jH PERE IN, 
‘Units, , *mMemmalized’, .. 
ree er ounecolor: oO 752941 On oz 94t Oe Joos ll, 
bpPos@eaon . |0.35c429 0.342995 0.175123 0.6483092), 
‘String’, 'Blade Twist (deg)', 
“Style. GeExt’, : 
“lage, -seaciclext2Z” ) ; 

a= "Witeenctrolt Parent 7H PERE IN, 


60 


‘Units!’ “normalized 

'"BackgroundColor',{1 1 1), 

‘Position’, (0.549918 0.342995 020981997 G2 0483s09Zi— 
‘style’, edit 

‘SULing:, EwWise, 


‘Callback', 'twist=get (gcbo, ''String'');S _USER_INPUT.twist=str2num(twist 
| gers, 2s © 
‘Tag’, * Edit@exti"); 


d= Utcontro (gba remit, PERE IN, 
Use Se ee 
"BackgroundColor', [0. 752941 07752941-07752941); 
‘Position’, [O.s58429 02280193 0.175123. 0. 6493092), 
'String','Blade Wt-Aero (lbs.)', 
TSeyle .~ text. | : 
'Tag‘, 'StaticText2'); 

d = uicontrol('Parent',H PERF_IN, 


‘Units', 'normalized', 

‘BackgroundColor’, [la 1], 

"Position', [0.549918 0.277778 0.0981997 0.0483092), 
eile) ena 

*String',wblade,... 


'Callback', 'woblade=get (gcbo, ''String'');S USER INPUT. wblade=str2num(wbl 
ade);', 7 7 
'Tag', 'EditTextl1'); 

Uicontrol( Parent, leartkhra LN, 

"Units', 'normalized', 

meeckarcusdcelon POs 752941 0.752941 0.752941], 
bEOsicieon Fi. Soe429 0.21256 0..17512350.0463092]), 
'String','No. Blade Elements' 

Seville" 9 ' text, 

Sladen cvactlelext2) ; 

uicontrol ('Parent',H PERF_IN, 

pUniTes "aermali cea 

Pack aroundcoken Pe Wes aoa | 1}, 

"Position', [0.549918 0.21256 0.0981997 0.0483092], 
PSevle!, “edit ‘, 

 SEEIMMGweIbe,. .«. 


d 


Q. 
II 


‘Callback', 'nbe=get (gcbo, ''String'');S USER INPUT.nbe=str2num(nbe);', 


bhag 7, EQLeETextl”)> 

uicontrol('Parent',H_ PERF _IN, 

Units “y Seema lie sae 

*BackgroundColor', [0O. 752941 0.752941 0.752941 ]2 
“Posieion ,70.711948 0.916264 0.175123 0.0483092], 
rotLing’,; Auxillary Thrust (lbs)', 

Style", ‘text! ye... 

‘Tag. p.- oCalucrextz2:): 

Niconerol( barent) HePERE IN, 

'Units', 'normalized', , 

“Ba€KGLOunacolor. gamle le i] > hoa 

"Position', [0.90671 0.917874 0.0981997 0.0483092], 


Q. 
H 


oF 
H 
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Pobyle*, edit, 
Potrdine, LAUK).<.. 


"Caliback'", “Taux=get (gcbo, “*String'');S5_ USER_INPUT.Taux=str2num(Taux) ; ' 
Pee: 


‘tag’ , EaQrtTexti"); 


Q = uscont rowm~erarent',H PERF IN, 
WnLtes’ ,-normalized’, 
Wee ckqroundecwier’ PLO. 752941 0.752941 0.752941], 
"FontSize',6, . 
PeOsition- 710. 711948 O. 851852 02173486 0.0483092), 
"String', ‘Flat Platte Area (ft*2)’, 
~SeVie us Gext ", F 
hag; Stakicrext2 ) ; 
Q7=suacencerel( Parent 7 HoPERE IN, 


"Units' Soares 

Wescreroendcolor: a oe 1}, Bouts 

"Position', [0.90671 0.852657 0.0981997 0.0483092], 
"Style’,#edit"', 

Potrina: ,ALh? 


"Callback', 'Afh=get (gcbo, ''String'');S_USER_INPUT.Afh=str2num(Afh);' 


Tag’, Ed2ztTexti'"); 


Ae aeitheoneroin Parent ,H PERE IN; 
'Units', 'normalized', 
eee eroundGolion | Oe 752941 0.75294] 0.752941), 
"Position', [0.711948 0.78744 0.173486 0.0483092], 
“string”, 'Vert. Proj Area (ft*2)*, 
oevVlew text-, : 
‘Tag’, 'StaticText2'); 

av sutconerou, Parent ,HOePERP TN, 


‘Units', 'normalized', 

ieee candeolon: va) (i Coma 1}, : 

‘Position', [0.90671 0.789855 0.0981997 0.0483092], 
'Style','edit', 

“Strang !,ALV, 2 


"Caltimeack |) tAtv-get (gcbo, °'String’’);S USER INPUT Afv—=str2num{Afv); ¢ 


"Taig', ‘EdwtTextl'); 


Q = urcont rem Parent’,H PERF IN, 
‘Units', 'normalized', 
eeitcroumieslor' At ON: 752941 O..752941 027524)... 
"Rostevlon 7 (0.711948 02723027 0.174577 0. 0483092], 
‘String’, 'Vert. Tail Area (ft%*2)', 
“Seyle'; *text*; , 
"Rada, Staticrext2’ ); 

Ge swlecontrolkerarent ",H PERE IN, 


‘Units', 'normalized', 

ieee eroundcolor' plod il, erste 

"Position', [0.90671 0.724638 0.0981997 0.0483092], 
‘style’, 'edit', 

“Se ringd.,oVeLt , ss. 
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'Callback', 'Svert=get(gcbo, "String ~);S USER INPUT -Svert—-strZaumisvert 
ear 2 © 
‘Tag’. Edittexcl ~)- 
d= uicontrol() Panent steps ner 
‘Units "> normalized, 
‘naetenoundecion Ole 752941 O- 7529416077 5254 1 \- 
"Position', {0.717948 0.658615 0-1 7457750. 048309712 
‘String’; Verto. fails spans te) 
“Stvile*,. text’, 
‘Tag’, 'Statlcrlext2”); 
ae= uicontrol(’ Rarene 7 uerERn SIN? 
‘Units }; normalized, 
Eaearoundeolor: corey. a 1}, bang 
‘Position', [0.90671 0.65942 0.0981997 0.0483092]j, 
potty le 2, edie’, 
Lob ring |, DVeELe, 


'Callback', 'bvert=get(gcbo,''String'');S USER INPUT.bvert=str2num(bvert 
Were oy vs ote 
‘Taq, 7 ECQierext ll,’ )> 
de sureoneLoh( Patetie: PERE. TN, 
"UnLts-, normalized’, 
pac haroundcolen: FeO 752941 Oo 75294 1 0.752941) 7a 
PPosirtron.  (O-. 71949 10. 59420370 .174577_0.0483092)- 
(oering ., Vere. “Tall cL’; 
‘Style | Gext '; 
Tag otacticrext2?). 
ge ULControl@rarcne, H PERE IN; 
PUMtese se NOrmal: Zed, 
ipAehgrounacelons Pm (lea I 1], 
ePOSiItE1LOn 7) (0.9060 IO. 596618 (0,0981997 0.0483092)- 
Potyvle!, “edit 
'Sering , CLVert, 


‘Callback', 'CLvert=get (gcbo, ''String'');S USER_INPUT.CLvert=str2num(CLv 
SUSIE) A 
'Tag', ‘EditTextl1'); 
d = ULGenELOI() Parent, 7 Web ERr EN, 
‘Units', 'normalized', : 
"“BackgroundColor', [0.752941 0.752941 0.752941], 
‘Position’, [O.711948 07529791 0.174577 0.0483092], 
Poierang 4. Vere. ball (GDo', 
PoeyLe > text’, 
Lag’ StatLecrTextZ > 
Cea sUrcontrol( Parente mM eERE IN, 
Omitts*,*normalized*, 
Teg ounaeo lor: Py camel 1], 
"Position", [0.90671 0.531401 0.0981997 0.0483092)], 
Style’; ear. , 
String ;CVOVerc, 


‘Callback’, 'CDovert=get(gcbo, ''String'');S USER_INPUT.CDovert=str2num(C 
DOVertu) 2. 


63 


Placa BALGText1 *).; 


dj-sumcontrol(,Parent*,H PERF IN, 
"Units', 'normalized', 
eee Per oundeoler: A eele 752941 On 329415 On 7529417 
Positron jt0.71194870.465378 0.1745a7 0.0483092], 
PoOELIng , HOrlz.. Tail Area (ft%2)*, 
‘Styles, text; ; 
'Tag', 'StaticText2'); 

dJ= Uteomterol( Parent’, ha PSRE_ IN, 


‘Units 7 nobmalazed’, 

Ie Larcundcslor’ po ke all 1], 

"Position', [0.90671 0.466184 0.0981997 0.0483092], 
‘Styae', edit"; 

PSEcCInGg 250051 2,6 


'GCal¥back |, Snorrz—get (gcbo, "String'*),;S_USER INPUT .ShoriZ~strZ2num(Sho 
biz) eee 
‘bag: -EGLETeXt1™) ; 
uicontrol('Parent' H_PERE_ IN, 
Wmits*, normalized', 
Snecraroundcolon LOR 752941 0.752941 0.752941], 
*POsSition’, [0.711948 0.400966 0.174577 9070435092], 
ZoeLIng . Horiz. Tail Span «(£t)i"; 
vot yle., text’, 2 
hag, Static lextZ'); 
Cia eterecome ho era Genin ter rnr LN, 
MWhtes so normalized’, 
"BackgroundColor", |] 2 1), 
"Posieron:, (0.90671 074033682 0.0981997 0.0483092), 
‘Style’, 'edit', 
Rrs)c ca liate (ergo) pYoy aalyrs seneren 


d 


ae neem Jee oy ering’) 7S _USER_INPUP.enoriz-strenBm Phe 
riz) 
‘Tag! 'EditText1' 1 


Ga sUtCOnEEGm@erarent ,H PERF IN, 
Unites’ normalized’, 
TemenaroundGelor PLO. 752941 O27 52944 0s 7 o2o4isl, 
"Position. , (0.711948 0.336554 0.174577 0.0483092], 
TSeLIng.- HOGiz. tarl Cl 4 
“OeyLe, text, ; 
‘Tag’, 'StaticText2'); 

e- USconeroljeeranemt 7h PF ERE IN, 


"Unit s ¢ Mee iiced’. 

BREE EOUHACOLOr: ie 1), ie 

"Position', [0.90671 0.338164 0.0981997 0.0483092], 
Style", “edit, 

‘SeErng 7CUVIOrlZ, i... 


"Galiliback';, Clhoriz—getigebo, 'String'*);S USER INPUT .CLhoriz=str2num (Cc 
LhOrasz op; : 

‘Tage, EGLtrextl" ); 
Qieeteontrol(, Parent ,;H PERF IN, 

"Unites -, normalized’, 
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'BackgroundCollommai0i 529414037 5204 204 
'Position',; (0.711948 "0.272141 0.174577 0.046305 
'"String'; Horuz a lasechan, 
'Style','text', 
'Tag', 'StaticText2'); 
d = uicontrol( "Parent 7HoPERr SiN, 
'Units', 'normalized', 
Racker ounaeanas: julelepel: 1}, ‘ 
*"Position’, [0.90671 0.272947 0.0981997 0.0483092), 
iotyle*., ‘edits, 
Siang, ,AcDOMmoOrdez a. 


peal lback*;, *CDehoriz=get (gcbo, “string ), ssUSERELNEU TT Clone -1Z-s5er2num 
MEDONORTZ) fe 4. 3. 
hag 9. BQgthe xed, i); 
Hepick = uicontrol( Parent [he ERP min, 
'Units', 'normalized', 
Pe engreundcolor® Pip a 752941 0.752941 0.752941], 
"POSItion 770% 71194820. 207729 02297872 0.0463092)', 
‘String’, 'Horiz. all Wander Main Rotcr Disk"; 
"Style*, "checkbox", 
‘Tag’, Checkboxl”, 
‘'Value',0O, : 
'"CreateFcn', 'taildisk=2;' 
'Callback',[. 
eii gee (Hepisk;, 7 value) —--17 = 
"taildisk=1;' 
‘else,'". 
"taildisk=2; ' 
"end, ' 
'S_USER_INPUT.taildisk=taildisk;'}); 


em ULCoOntrol( @rarent | WH PERE GIN, 
‘Units’, 'normalrzea’ 
"Callback", 'performance input_fen back', 
PE ONCS1LZe yelc, : 
"FontWeight', 'bold', 
~POSlCiIon FU. 293, 00555906 0.163666 0.0805153)], 
PStLing’ | ——<eepack , 
“kag’', Pushbutconl”); 
ee sULeCONt roll(s barent , HOPERP SEN, 
PCalllback | periotmance Gmput fen print’, 
'Units', 'normalized', 
r£OneSL Ze vel, 
"FontWeight', ‘bold’, 
“Position’7 (0.51860370.0550556 0. 163666 0,.0805153}, 
PSELING’. Pelmex screen , 
"‘Tag*’, Fushbuteonz'); 
d= UI Control Farent jie ERe arly, 
Callback Fa pe rlormance Inome ten enx 7... 
'Units', 'normalized', 
*Fontsize*s iz. 
"Fontwelgnt.,) bold’) ... 
*POostelom 7 (0.525914 0.0555556 0.163666 0. 0805153), 
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Est ting), Cancel)", : 
'Tag', ‘Pushbutton3'); 
de eurconcror ( barenm 7H PERE IN, 
lUntcs = normalryzed' , 2. 
"Callback', 'global REGIME PICK 
S_PERF_INPUT, PICK=0; REGIME=0;S PERF _INPUT=S_USER_INPUT;performance inpu 
tl ienecont , : 
"FontSize’,12, ‘ 
'FontWeight', 'bold', 
SEOs UETOM 7mlO 7 s5225  0.0559900 0.163666 0.0805153], 
foerindg 7 Ccontinue.>>'7 
'Tag', 'Pushbutton4'); 


dost gmimnwease 7 HeDLOk ;HeDion), 


assignim( base ; Ss USER INPUT’ ,S USER INPUT); 
assignin('base','S PERF INPUT',S PERF INPUT) ; 
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APPENDIX I. PERFORMANCE _ INPUT_FCN.M 
Switchyard Callback function for the performance_input.m GUI function. 


Puncelon per tOmianecen! Nou eEen (1c etc 


$ Switchyard Callback function for performance input.m 
% JANRAD 98 VERSION 4.0 


global H_PERF_IN H IT METH S PERF INPUT 


if nargin, 
switch Action 
case ‘'cont' 

fieeisempey (get field (SsPeRr INEUlL | eA) |... 
Geetieldits se ERF | INPUT, 'temp')|. 
getfield(S PERF INPUT, 'Vin£') |... 
FE ee ura ot atiae 
getticldis PERE INPUL, omega.) i-.. 
jerrrelate PERF EN POT trace a) [oes 
jetevelatoy PERF INPUT, “thetrao”’) |... 
Gettielce. PERF INPUR owing: jell. 
getfield(S PERF INPUT, Vbwieng.) lees 
detii celal PERE Set Ciawang:  ) slime 
geetreld(sS PERF LENPUT,  CBowing -) lee. 
getfield(S PERF INPUT, 'ewing')|. 
Geericla(seeeRe iNeUl, atoid.)a|i- 
dec Ere lich( Sm PERF SEN PU c 1). mee 
gett eld soPERr INPUT, ae... 
getfield(S PERF INPUT,'R')]|... 
getfield(S_PERF_INPUT, 'e') |... 
Gebrretd(s PERE SINEUL, Grip) |.o. 
Jercere lation PERF INPULS renora )) |p. 
Gectrelaicn PERF_ ET oe re alles 2 
Getivelata PERF PVeU TT “trst* ) lin. 
getfield(S PERF INPUT, 'twist')|... 
feeirela(se PERF INPUT, 'wblade')|... 
dJectield(s PERF _ PNEOR, noe) |. 
getfield(S PERF INPUT, aod Gf 0 bic i eee 
getfield(S PERF INPUT, 'Afh')|... 
Getivellays PERE INPUT, *Afv*) |... 
getfield(S PERF INPUT, 'Svert')]|... 
Jetrrveold(s PERE INEUL, “bvert ') |). 
getmreld(S PERE INPUT? "Civert") | >.- 
Gettield(SePERr INPUT, “Ebovert ) |... 
Gerrveld (SePERRRINEUT | Shorirz:) |... 
Getricla(s PERNEINEUD, bhoriz') |... 
Geer relic PERE SINPUL | Cineriz ))| 22. 
Ggetivelld(s PERF INEUE CDchoriz’)); 

empty boxes 
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else 
itebat lon method 
elose mw Henke EN) 

end 

case ‘cnx’ 
performance input 
close (gcf) 

case 'back' 
analysis 
eTose (Hh _PEREN) 

case, oprante 
set (gcf, 'PaperOrientation', 'landscape') 
Setuiger, PaperPosition’ , (oe. 5 10e7.5)) 
pring —awine 

case ‘return' 
janradg8 
close all 

case ‘quit' 
Gud tegud 

case ‘about' 
about _janrad 

end 

end 
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APPENDIX J. PERFORMANCE OUTPUT.M 


This file creates the GUI to display the calculated results from a previously saved 


input file or newly created user input. It is call in Perf.m. 


function per formance soucpur,) 


%$ GUI window to Display Janrad perfromance output. 
% JANRAD 98 VERSION 4.0 


%$ This is the machine-generated representation of a Handle Graphics 
object 

$ and its children. Note that handle values may change when these 
objects 

are re-created. This may cause problems with any callbacks written to 
depend on the value of the handle at the time the object was saved. 


% 
% 
% 
%$ To reopen this object, just type the name of the M-file at the MATLAB 
$ prompt. The M-file and its associated MAT-file must be on your path. 


load performance output 


global COUNT H_PERF_OUT S PERF OUTPUT S_USER_INPUT 
H datain H dataout H vecdata H check] H_ check2 Hecheer > 


COUNT=1; 


HO PERE OUT = figure *Units', 'normallized', 

'Color' PO eS OreeeOr sore 

"Colormap',mat0O, 

"Name', ‘Performance Output', 

"'NumberTitle', 'off', 

'PointerShapeCData',matl, 
PpOstetone 1 -On008125 0705625 0.954688 70ncocr 5, 
lnel Ghages) a ay lee, 

HEMehU( Parente 7h PEREsoUT, 

‘Label', 'JANRAD Options', 

'Tag', 'uimenul'); 

uimenu('Parent',b, 

(WallbDack') pextormance output fen quit’, 
"Label', 'Quit JANRAD', 

'Tag', 'JANRAD Optionssemaimentm ee 

role Parent',b, 

'Callback', i peRcOnnames WOULpMC Len Tetum... 
‘habel", ‘Return to Begining", . 

'Tag', 'JANRAD OptionsSubuimenul ' ye 

une Parent',b, 

‘callback’ "per fommancemeurpar ton deliCawinpur 7 
"Label','Change Input Parameters', 


oO 
tl 


Q 
II 


Q 
lI 


Q 
I 
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It 


tl 


'Tag', ‘Subuimenul'); 

uilmenu('Parent',b, 

"Callback', ‘performance output fen about', 
'Label', 'About Janrad 98 ...', 
'Separator','on', 

'Tag', ‘Subuimenul1'); 


Hikconitaoh (Parent ',H PERE OUT, 

‘Unies. Pena li ceo. oes 

"BackgroundColor', {0. 752941 O.7 5294170. 752941), 
“POntolze =i, .«<% ; 

“PoOsttelen. , | 0205310966 00-925659 102327 332.0.0383693], 
'String', 'Fuselage Drag (lbs.)', 

fSeyle"} text, ; 

'Tag', seaticlextl"); 

ulecentrel( Parent ,;H PERF OUT, 

“Unises 9 Mommalicea”, 

"BackgroundColor', {0. 752941 0.752941 0.752941], 
"Position', [0.376432 0.925659 0.0981997 0.0383693], 
Long ~smE ERE OUTPUT Druse, 

J stewie? peas o 

BLagdi, Sctactterext ll’) > 

Preentret barencte | 1  PERMBOUd, 

 Ubolses)y eee ee 

meoreroundcolor Pe ie Ole 752941 OV 752941 0.752941 }- 
“FoneSize", 10; 

POST eon. alOn0s LOVGGrOec7 290270. 327 532 OF 0363693], 
‘string ; ROLores Dragmelbs.) *; 

nocyle > text: ; : 

lag, otaticrextl” ) > 

nicontrol(” Parent”, H PERE TOUT, 

Unites‘; ‘orem e ‘ 

"BackgroundColor', [0. 752941 OTfo29al 0.732941), 
HPOosition., (0.570432 0.8/2902 0.0981997 0.0383693)], 
(oeEing =, 5 PERF (OUTPUT. Hrotor, 

'"Style' meena : 

tTag j. SsStaticlextl "); 

Wicome rel Parente aero (Oui 

‘Units *, normalized < so... 

"BackgroundColor', [0. 752941 0.792941 GesjoZ94biy 
'FontSize',10, : 

‘eesitron .|0. 0310966 0261/7746 0. 327332 0. 0scs673]7 
POEcang. . Wang bite. (lbs...) ", 

Tpey le. text", 

‘Teg’, StaticTextl'); 

He Ourno le Parent, Hv EBRE OU, 

‘Unites’, ernaiaeed!. 

"BackgroundColor', [{0. 752941 02752941 O275Z294|- 
‘Position’, [0.376432 0.817746 0. 0901997) 0.035s4693)], 
"String',S PERF OUTPUT.Lwing, 

‘Style’, ‘text, : 

"Tag, Statieclrextl’ ); 
HEcCentErol(] Parent jHeeekr OUT, 
‘Unies > mormalized', 
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'BackgroundCol or aie 529410 2941 Ce a2 94 
"Fontoize ' 7 107 nae 

‘Position’, [02.0310966..0.762590.327332 0. 03383652), 
‘String’; Wingeprag ibs.) s 

Motyle! tex =, : 

‘Tag’, Stataclextl ); 

UlLCOonEEOU( Parent 7 Hwee ReLOul, 

‘Units ', norma lizea = 

‘ packeround@olen: Puke 752941 0.752941 0.752941), 
‘Position, [0.376432 0.76259 0.098997 0.0sea509osie 
'String',S PERF OUTPUT. Dwing, 

“otyle 7 text. 

‘Nag’, ‘StaticTextl’); 

MecOnCEOM Farent 7H PER OU, 

‘Units’, normalized’ 

TeReLerounde sien: , GOR 752941 0.752941 0.752941], 
"FontSize',10, 

'Position', [0.0310966 02709832 0.327332, 0.,0383S6e93), 
Moe Ling we Hhorezontal Tall litt (1bs.)"; 

PoOtVle, text, 

‘lad jy Stacie lrextl'); 

UTEOnL EOL, Parent ,He@EBRn OUL, 

(Unkes  , normaliazea 4 

Tene a coUndcoNor: p10: 752941 0.75294 IWA0e 75294); 
"Position', [0.376432 0.707434 0.0981997 0.0383693], 
PoGilng oo EERE OUDPUL. Lhoraz, 

petvile’™, "text', : 

‘jag’, 'StaticText1"); 

UacOnLEol( Patent ,H PERE QUT, 

Wnts; normalized 

Beerercundcanen| Pa (ole 752941 0e@752941 0.752941], 
sHonts1ze'* , 10, ; 

'Position', [0.0310966 0.654676 0.327332 0.0383693], 
oering’ , 'Horizontal TaaieeDrag (lbs .)', 
psuvile'’,'text'; ; 

lag: ;,; Gtaclelextr )- 

MEcontrol( Parenu: pH sesh r OUly 

"Units', 'normalized', 

een roundesilor? LO. 752941 0.752941 0.752941], 
‘Position', [0.376432 0.654676 0.0981997 0.0383693], 
'Style', 'text', 

Pieri ngs Een Oud PUT phoriz', 

"Tags. staelcrexctl' ); 

ULcControl ("Parent’,H PERF OUT, 

"Units eld eam ws 

 BacKgroundcoiees plc0:. 752941 02752941 0.752941), 
PHONES 1 Ze 20; 

*Position’, [0.0310966 0.59952 0.327332 0.0383693], 
'Steing. = Vertical fall thart (ibs. ),", 

‘stvle", texe., : 

"Tag,  otatrerext) ); 

ulcontrol( Parent’, HUEERE OUT, 

‘Units’; ‘nermalized’, 

‘Bacharouneenlor’ AKO. 752941 Og o2941°0. 752941),5 
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beocmeron .» (0-376432 0.59952 0.0981997 0.0383693], 
'String',S PERF OUTPUT. Lvert, 

poevlew, text; F 

lade, statlerextl ); 

Mmcomeurol Parent: | H PERF OUT, 

hUmnites', normalized", .. 

meee grcundeslor® , (0. 752941 Oa o294950.752941)- 
“FontSi ze’, ie, 

'Position', [0.0310966 0.546763 0.327332 0.0383693)], 
String’, ’ Vertical@Tailebprag (iiss.)"; 

‘Style | text’, 2 

‘Tag, Staticrextl*),; 

NEeconerol(sraréenk., Hh PERPSOUr, 

‘Units', 'normalized', 

Dea chorounde alos" pile 752941 0.752941 0.752941), 
'Position', [0.376432 0.546763 0.0981997 0.0383693], 
Potiing 7 oo bRoseuULPUL, DVert, 

‘oe VLe*, text", : 

‘Tag 7, sStatLierextel'); 

uicontrol ("Parent*,H PERF OUT, 

‘Units', 'normalized', 

Menergrcunacelar: mILOs 752941 OR 1529418077 529415 
Ponto ze’ ; 10; 

‘Position, , [0;0310966 0.491607 0.327532 0.0383693), 
"String', ‘Tip Path Angle Sees 

VouvyLe , text', 

‘lag |, sovaticrTextl'); 

Mucont rol Parente ,H PERF OUT, 

Whew ei Sed! . : 

NBacharounc Colon LO 752941 0.752941 0.752941], 
'Position', [0.376432 0.491607 0.0981997 0.0383693], 
'String',S PERF _OUTPUT.alphaT, 

"SEVLe’ 7 text. , : 

‘Tag’, Staticlexetl ); 

uLeconErol( Parent ,H PERF OUT, 

sUnLts 3° normalized’, 

pacha roundcewer’ ples 752941 0.752941" 0.732941), 
'FontSize',10, 

"Position ’ "0. 03510966 0.436451 0.327332 0.0383693], 
'String', 'Rotor Coning Angle (deg)', 

“Sie ylen, text 7 aise. 

lage ocaticlextl") 7 

HEGCOnE LOI Parent 7H PERE OUT, 

‘Whits', mormalized', 

PeeelenoundcolGr: LO: 752941 0.752941 0.752941], 
'Position', [0.376432 0.436451 0.0981997 0.0383693], 
‘String',S PERF OUTPUT.betao, 

‘Styvle%,"text', : 

‘Tag’; staticrexel”) ; 

uicontrol ('Parent’,H PERF OUT, 

"Units', ‘normalized’, 

Geter Gr oundeciior' Po 752941 0.752941 0.732941 \, 
"FONE SI Ze, 10, 

'Position', [0.0310966 WesesoIse0.s2 1852°0-,0305095), 


a2 


"String, 'Vocatiionmot Mainulhn ius te 2 meee 

Style" rae eee 

'Tag' Diceaevere sea 

uicontrol ('Parent’ ,H PERE OUT, 

‘Units', 'normalized', 

‘'BackgroundColor', (0.752941 0.752941 0.752941], 
‘Position’, (0.376432 0.383693 0.0981997 02038369siy 
"String | sere RE OUsPUT ro, 

"Style, texte, : 

(Eags) Stake rexed ); 

Udconerol( | Parent ji PERE Our, 

‘Units' noire cau 

"packarounaeolen ah). 752941 0.752941 0.752941], 
THOME S1 ze 70 =. 

peOs i Gi1on (O07 0s5U09G6 0.328557 0.327532 0.03336o31, 
Pouring | lst Gat. .cyellile Term = Al”, 

UStylier text). ten 

flag ,; otabiclextl.); 

uicontrol('Parent' ,H_PERF OUT, 

“Units Weta cee 

Bee ar ouncenion: Or. 752941 0.752941 0.752941], ‘ 
'Position ,10.376432 0.328537 0.096199) 0.0scseos)] 
boeing jo PERRSOULEUL.thetalc, 

mSeyle', text’, ‘ 

hag, Statler ext); 

WicoOnt rout barene  7HerPERE IOUL, 

“Units ' aac 

Demo igroundeollens 7 bOe 752941 @.752941 0.752941), 
BECOMES Ze!. 10; 

epOstelon —swOe0st0966 O327338) 0.327332 0.03es59s), 
Boer InNG' | 1stehong. Cyelie Term = Bilt, 

Style! “text', , 

Wiad OcatLerext1l' jys 

Uteonr rol ( Parent 7H EERE SOUL, 

“nats, ee : 

mene Ground@alor: a (Os 752941 Oe, HES) 0) 152,94 le 
BPOSTENON 71053764352 0.273381 0.0981997 0-0333693)- 
/St¥ing’,S PERF OUTPUT.thetals, 

’Style' icerem 

wuad-—, otatcrcrextl!) + 

HvcOuELoi Parent ,H PERE OUT, 

Lontes! pre ease 

meee crounacoon Aen 752941 0.752941 0.752941], 
'FontSize'’,10, 

epost tiem, (0.02275 0.92326! 0.327332 0.0383693)], 
'String','Collective Pitch @ .7 r/R (deg)', 
eoeyvle’, text’; : 

Taga StaticTexte' ); 

UTeconeroR@Paremic. | Heb mnr OUT? 

‘Units’, ‘normalized’, 

Tach arounaeelon Oke 752941 0.752941 0.752941], 
"Position’, {0.859247 0.925659 0.0981997 0.0383693)" 
(oErIng, 7ouPeRr OUTPUT. Cthetao, 

YStyle’, "text"; 
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It 


phag > StaticTextl*),; 

Peont rol Parent: 7H PERF OUT, 

Mnivts*, ‘normalized’, 

Mpackercunacolor' POs 752941 0.2.75294140-75294 1), 
MBbonts1 ze", 10, »... 

Position’, (0. 512275 Omer 0904 0.327332 (0. 0383693), 
‘string’, solidityeisigma) +, 

vecvyle", (texG, 2 

Tag", otaticTexel)y 

icone rol Parent’ ,H PERE OUT, 

‘Units', 'normalized', 

Pome rorounceotion« a Or 752941 O. 752941 0.752941) 
*Position', [0.859247 0.872902 0.0981997 0.0383693], 
Pot ng Ss obeRe sOUIPUT (solidity, 

rstyle eee ca 

'Tag' i ieeationeer) ; 

WIGORCLON( Farent 7H PERE OUT, 

Units: er malized". 

eae rGEoundcolen! Or 752941 0.752941 0.752941], 
'FontSize',10, 

ppesition'. —0. 512275 902815348 0.327332 0.0383693)- 
‘SGering , isk Loading (lbs. /£E72) "| 

"Style'; 'text', : 

"Tag', 'staticTextl'); 

UPeont rol Parente ji EBKE OUl, 

‘Units', 'normalized', 

Ppackgrounacolor', (0. 75294150. 752941 0.752941], 
"Position', [0.859247 0.817746 0.0981997 0.0383693], 
Poe rind eo oe ke oUCPUT DL, 

potyle text | : 

lag, “staticlextl" ); 

MMeoneroh( Parent 7H PERE OUT, 

"Units', 'normalized', 

eer Oemaeoliot | 10. 752941 0.752941 0.752941), 
*FontSi ze", 107 2: 

"Posation’, (0. 512275 O2 760192 02327332 0.038369]; 
'Sring', ‘Figure of Merit', 

"Sty#e', 'text', : 

whacG = Stacicrextl”); 

Dreoue Lome rarene ,Horshh sour, 

‘Units', 'normalized', 

ere rcundcoLer Pe (8 752941 0.752941 @7752941)5 
‘Position', [0.859247 0.76259 0.0981997 0.0383693], 
suming’ 75 PERE OUTPUT .EM, 

Style, text- : 

-fagey Staticlexel"); 

UleGenerol(@ Parent’ ,H PERE OUT, 

Waves, normalized’, 

memener cundcoron 7 lO. 752941 0.752941 0.752941], 
tHOMNtol ze ~lOe. . sk. 

“POsleron 40.9227 500.70/434 0.327332 0.0383693)], 
"String’, ’ Cl/Ssiagmae 

"Style", ‘text, A 

‘Tag 7 StaulcTexGis; 
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ulcontrol (Parenter ERE .OUT, 

‘Units’, ‘normalizea ; 

‘Backarounaeolons Oe 752941 0475294 12.027 52944) ) eee 
‘Position’, [0.859247 0.709832 0.0981997 CO 03s3693|— 
'String',S PERF OUTPUT.CT sig, 

‘Style "| Eexe 7, : 

'Tag', 'StaticTextl1'); 

UlcCont-rOIl( Parent | HalhRE seu, 

'Units', 'normalized', 

Bacis couneenite ey Pe Os 752941 0.752941 0.752941], 
'FontSize 710, 

"Position? (0-512275 0.652276 0.327332 0-0383693)- 
‘String: , €O/ Siama - 

"Style: , text; 3 

‘Tag, Staticlextl ©); 

icone rol ( Parent Har sie sOUdr 

"Units', 'normalized', 

pacha rouna@oller: eiiO 752941 O2752941s037529 4 ia\7 
"Position ,»p(02859247° 0.654676 60981997 0.0383693)], 
UString), SaPoREwOuUreurl.€O sig, 

Sotyvle" > "Cexe., ; 

"Tag: etaticrexe lt”); 

UTCOnoEC Werarent nme ok JeUmy, 

‘Unres. Noticed 

mpaergroundealor: PO 752941 0.752941 0.752941], 
FOnMtS2 Ze" 7.0, 

PrOSITtLON 70. 51227570-.597122m0Rs7 7332 0. 0383e°s |, 
Focrang.,. Ch, orgma 7 

PoLVle: 74, text, : 

'Tag', 'StaticTextl'); 
MuGonerol Parent “7HuEPERE IOUT, 

"Units', 'normalized', 

‘BackgroundColor = 0. 752941, 0. 732941 02752941), 
“Position , [0.859247 6.59952 0,0981997 0-0383693)] , 
‘'String’,S PERF OUTPUT.CH sig, 

"Style' ,'text', ; 

tlag =, otatichext ll"); 

Leone ro) | Parent” ,HeeERE .OUT, 

fUniese ore lizedis 

ipachynoundcolons LOR 752941 0.752941 0.752941], 
LFontoize’ ,l0, 

BeOste Lom Oso 2 275 02544365. 02327332. 0..0383693), 
Paering 7 LiseMacheNc. of Advancing Blade, 
"Style', 'text', 

“lag; statied extl* ); 
BLeonerol( Parent 7H ePERE OUn, 

"Units', 'normalized', : 

“packgroundcolor* [0.752941 0.752941 075294117... 
“POSLELON 7 0.859247, 0.546763 020981997 020383693]; 
String | eben eeoUrPUL.Machtip, 

‘style; texe = ; 

“Mad, octatdelextl .)7 
Utcontrol( Parent 7H EERE GOUT, 

“UNLESS”, NOrmalrzed’, 
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pate nmorTouncdcolon, (0.792941 0.752941 0.752941], 
'FontSize',10, 

Ubposa tion), (Oe0l2275°0.4289209 0.327332 0.0383693], 
‘String’, ‘Advance Ratio', 

Toeyvie 7 text’, : 

Tago, StaticTextl’); 

Hicontromt Parent ,H PERF OUT, 

Luntes*, normalized’, .. 

mee arouideolar! mitOe 752941 O..752941.-02..752941)|> 
'Position', [0.859247 0.491607 0.0981997 0.0383693], 
'String',S PERF _OUTPUT.mu, 

Poe vle. 7. text 7.0... 

'Tag', 'StaticTextl1'); 

UteGont rol Parent 7H PERF OUT, 

‘Unites 5 Pema seat, 

"BackgroundColor', (0. 752941 0.752941 0.752941), 
'Fontsize':,.10, .. 

"Posh tion’, (0.5122 7/5 07434053 0.327332 0.0383693),; 
‘String’, "Rotor Thrust Required STPP (lbs.)', 
"Style', 'text', P 

‘Tag’; Staticlextl"); 

Heconceol(. Parent 7H PERF OUT, 

Minacs > ’nmormalizeda'; .. 

Acker oumicollon: 20). 752941 0.752941 0.752941], ‘ 
‘Position’, [0.859247 0.436451 0.0981997 0.0383693], 
Petra), Se bene eCUTPUL a, 

ncoyle', “text 7, 

‘Tag’, ‘StaticTexel”) ; 

MECONEEOMN\: Parent. 7H PERE OUT, 

ntes Weciaincea”. 

"BackgroundColor', (0. 752941 0.752941 0.752941], 
"FontSize',10, 

ME@sitelOn 0.02275 O50 1295) 0.327332 0.0383693], 
'String', ‘Rotor Power Required (hp)', 

'Style', 'text', ; 

‘lag’, StaticTexti*); 

HMiGentEoin barent 7H PERE OUT, 

‘Units', 'normalized', . 

eee roroumicalor®, (0: 752941 6.752941 0-75294 17, 
Position’. 102059247 sO4seso9s 020961997 0.03383693)], 
Strang ~oebeReecUrPUT.Protor, 

'Style','text', 

'Tag', 'StaticTextl'); 

UlecheEol\ banent «, 1 PERE OUl, 

‘Units', 'normalized', . 

‘patter cundcolor Plo 752941 0.752941 0.752941), 
‘Fonte ze, 10, 

VWPOswe hon lO co b227580.520139 4327332 020333695), 
Voereing=," Rotor Lorque (ft.—lbsa)i, 

Sty ke 7Cext. , : 

'Tag', 'StaticTextl'); 

NPeConerol( PanrcnG ,H PERF OUT, 

Witts-.  nNorknallZea*, .. 

"BackgroundColor', (0. 752941 06752941 0.752941), 
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"Position", (0. 859a4mm@0.. 328557 20. 0981997 507 0scs569sl7 
String 7s PEREPOUTEUT Orotor, 

‘Style 2" texts, . 

‘Tag. staeieclhexting.,; 


b = uicontrol( Parent ~H PERE OCU 
'Units', 'normalized', 
 BaekGmounieeien: Pao 752941 Oma Z941° 0. 7a794 tie 
'FontSize', 10) 2. 
Position |) (Ovo 2 1 S027 O9ten 0575727332 ..0, OScools, 
‘String Auxiliary fnruse (tps), 
"Style’, *texe”, 
'TaGi a ctattoerext ll”) ; 
Da WUlcemerol( harent |, HOPE RE eOUL, 


‘Units’ ep Nala e a ee 
pacraronneeie LOY 752941 02752941 0. 752941 )]- 
‘Position -,l0.859247 0.273581, 0,0961997 Oe0seso9s], 
'String',S USER_INPUT.Taux, 
woevie text’, : 
‘Tage 7. Statlelexti); 
neeneckie = UicontEol(” Parent ,H oPERE OUT, 
Unde Se a NOrmales 260 see. 
Pe aeaeoundcolort oe) Oe 752941 OF 70294 lO. 79294 \- 
“Positi1on', |[0.0310966 0°1774535 0.266776 0-04/961G), 
Peering’; “Save Input Data vas J. .:", : 
"Style', 'checkbox', 
'Tag', 'Checkbox1'); 
H datain = EE CINE EO ao tag EE OU 
Unies |. normalazea’— 
Tee a moUndcolen: eivie 1}, 
*"Position', [0.302782 0.177458 0.0981997 0.04796 lg), 
phOntsl ze’, 12, 
Potyle', "edit", 


Potrrng. ; see 

Hcaliback’ | [a.. 
PSCtilgcero, “String “,get(gqebo, "string = 7. : 
"set (H_dataocut, '"String''’,get(H datain, ''String” 9) 7)" 
rset (H veedata, '*sString’',get(H datarn,’ "Sterngue je 


eSee(Wecheecr Value! !) lite 
PSet(hecheck2, “Value 7): ," 
“see (H checks, “ Value’) | lee, cer 
“HoOrizontalAligqnment > right’, 
Pag 4 bOLtTextil™): 
Pee Ulcontrol({’Parent*,;H PERF OUT, 
“Unies , normalized", 
mpeercreundcolok Oe 752941 0.752941 0.752941), 
PEONCSITZe , 2224. 
EPOSTETON aOR 405892 0.179856 0.0981997 0.0479616], 
"Style', 'text', 
TSOeLInG. ~ o.mat. 
‘'HorizontalAlignment','left', 
"TaG’, stabaiclext. ); 
H check2 = UEecOnErCimiweatent pH PERE SOUT, 
‘Units "| ‘normalized’; 
‘Backgroundeslor: Ag Loe 752941 0.752941 0:°75294 ley 


(a 


bPosit 1 Ono ler oloowO.117506°0.266776 0.0479616], 
Potting, cave Output Data as ....', , 
'Style', 'checkbox', 
'Tag', 'Checkboxl1'); 
H dataout = uicontrol('Parent',H PERF OUT, 
‘Units’; normalized’, 
Dee ercmdcoior: Pe Oe 152941 Om7o2 941.0. 752941)), 
'Position', [0.302782 0.119904 0.0981997 0.0479616], 
'FontSize’,12, 
F Shicienlyals jars 
PStyle , text, 
'HorizontalAlignment','right', 
'Tag', 'StaticTextl'); 
b = uleontrol( Parent {HO PERE OUT, 
‘Units' ence 
ppeay ar ounce vo cu anon 752941 02752942 09752940; 
"PONG S1 Ze. > too 
'Position', [0.405892 0.122302 0.0981997 0.0479616]), 
PStrEing | .prti', : 
'HorizontalAlignment','left', 
"Style', 'text’, 
rag: Pesca crete ig 
Hechecks — wacont rol@eravent  /HerPERE OUR, 
'Units', 'normalized', 
‘BackgroundColor’, [0. 752941 0. 75294) 0.752941), 
"Position’, [0.0320966 0.059952 0.266776 0.0479616], 
"String’, ’Sawe Matrix & Vector Data aS ....', ... 
"Style', ‘checkbox", 
'Tag', 'Checkboxl1') ; 
H vecdata = uicontrol (*Parent',H_PERF OUT, 
'Units', 'normalized', 
Beeorounacolor’ Be Ok 752941 Oe goZz94) 0.75294) )\, 
'Position', [0.302782 0.0623501 0.0981997 0.0479616]), 
LBOne Size yilZ, 
‘Style’, text’, : 
"HorizontalAlignment','right', 
"Tag', "StaticTexite '); 
b = uicontrol('Parent',H PERF OUT, 
Unies. mOrhalicza’, 
Peo meroundeaitar: Pe Cele 152941 O.75294 1305752941), 
CEONCSI Ze 1212, 
'Position', [0.405892 0.0647482 0.0981997 0.0479616], 
SEnIng.. op .mac., : 
"HorizontalAlignment','left', 
woeyle',; "text ; 
‘Tag, , StaticlextZ’ ); 
be = Uecont rom Parent, PERE TOUT, 
'Units', 'normalized', 
“HOMES LZe Valz es 
"FontWeight','bold', 
"Fostti10on ylO-5 7263) 0.146283 0.140753 0.0815348)), 
YStrang’,'<< Back’, 
'Tag', 'Pushbuttonl', 
'Callback', ‘performance output fen back'); 
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b = uicontrol( Parente PERF OUT 
‘Units', 'normageazeda”, 
"Fontsize* 72, de 
'FontWeight','bold', 
‘Position’, [0.749591 0.146283 0.740753 30. @0853451) 
‘String’ 7 Optionsa— | 
‘Tag’, Pushbuteonss 
‘Cailback’;, pertormance Oucput fen ope), 
By="uLeontrol (Parent. | Here Rheoul, 
‘Callback’, ‘performance output fen print §;-.. 
‘Unvpes, normaliezea, 
“FOntS1 ze 71 2- 
PEOntWelgnt. | DOL. js... 
"Position';, (0.574468 0. 059952 02314239 -0. 0527576 )4 
'Stking', Print sereen 
‘Lag, Pushbutecen!™); 


assignin('base','H datain’,H datain); 
Assigmium@ base ,'H dataout',H datacut) ; 
assignin('base'’,'H vecdata',H vecdata); 
assignin('base', 'H check1',H checkl); 
gosignim( base’), Heycheck2 ,H check2); 
assignin(*base', 'H check3',H_ check3) ; 
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APPENDIX K. PERFORMANCE _OUTPUT_FCN.M 


Switchyard Callback function for the performacne_output.m GUI function. 


function) per foemance oulpucmicn (Aecion) 


6 Switchyard Callback for performance output.m 
% JANRAD 98 VERSION 4.0 


global 5 ae PERF _OUT S _USER INPUT Se PERE > INPUT S = PERE _OUTPUT Ss eMATR _ VEC. 
ee datain He: ~ dataout H_ vecdata H_ checkl H_ check2 H_ check3 
H outputfile H vecfile H inputfile 


if nargin 
switch Action 
case ‘'back' 
elose (Her eRE VOU) 
S PERE = INPUT=S _USER_INPUT; 
iteration method 
case ‘opt' 
at get (H eneck), Value” )—— 
S_USER_INPUT=S_PERF_INPUT; 
S_USER_INPUT.Vinf=S_USER_INPUT.Vinf/1.68894444; 
Ss _USER_ INPUT. thetao=s _USER_ INPUT. Gnhetac so -3; 
S _USER_ INPUT. twist=-S _USER_ TNEUT. Ewist*5 7-3; 
filenamel= =get (H datain, 'String' 3 
eval({'save ',filenamel,' S_USER_INPUT']) 


end 
Hf get(H check2, Value") ==1 
Pimlendmel Gee (Hdacadin, string )7 
eval(["!copy print templ ', filenamel,'.prf']) 
end 
v£ get (H check3,; 'Value')== 
unstructure3 
Pac mame ss vee enamer : ah a 
eval(['save ',filename2,' r psi vi theta betat alpha Tpsi Mpsi 
DMpsi dT dM dD cblade CL SDeie 
end 
options 


set(H inputfile, 'String', [filenamel,‘'.mat')) 
sect (aCe pUEL ile string. 1dee(H datacut, String’), *.prt*)) 
Ssen(Hevectile, Vstring), (get (Hh veedata, "String’),' p.-mat")} 
close VEE RE FOUT) 
case ‘print' 
set(gcf, 'PaperOrientation', 'landscape') 
SCL (OQCER, Paperrosition’,{.o «5 10.7.5]) 
Prink -dwine 
case ‘return' 
close all 
janrad98 
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case ‘delta anpute 
closewn PERRZOUT) 
performance input 

case ‘'quit' 
gulbe gui 

case ‘about’ 
about janrad 

end 

end 
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APPENDIX L. ITERATION METHOD.M 


This file creates GUI to select iteration method and display the status of JANRAD 
98 computations. Status comments are set in Trim.m and Perf.m. When computations 


are complete, this window 1s closed in Perf.m. 


function iteration method () 


% GUI window to select iteration method, start computational 
routines, 

%$ and display clock and performance method status. 

%$ JANRAD 98 VERSION 4.0 


% This is the machine-generated representation of a Handle Graphics 
object 

% and its children. Note that handle values may change when these 
objects 

% are re-created. This may cause problems with any callbacks written to 
depend on the value of the handle at the time the object was saved. 


% 
3 
% To reopen this object, just type the name of the M-file at the MATLAB 
% prompt. The M-file and its associated MAT-file must be on your path. 


load ieeraeGlon method 


global H_IT_METH H_NI H_AS H_AL H_GW H_BT H_BTR H_SOT H WSA ... 
oe STATUS H _STATUS1 H _STATUS2 
H_ GO H_ RUPT Baek ne RES H _MEN 
COUNT S _USER _ INPUT S_ PERF_ INPUT REGIME PICK 


COUNT=1; 
S_USER_INPUT=S_PERF_INPUT; 


H IT METH = figure(’Units’, 'normalized', 
'Color' FiO. 6 0.0) , 
enicrmap | mace, 
‘Name', ‘Iteration Method', 
“VNUMNDerTiLelre’, "off", 
"PointerShapeCData',matl, ... 
"rosurron , |—-0.003125 0.0625 0.954600 0.8625), 
ehag-,. Fig. jG 

HeMEND >] ulmenu( Parent") hn TT METH, 
"Label', 'JANRAD Options', 
‘tag’ , Wamenul'); 

c= Rate oe ee , H_MEN, 
'Callback', Teerarien method Eenmqueu., 
‘Label', 'Quit JANRAD', 
'Tag', 'JANRAD OptionsSubuimenul'); 
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Q 
II 


uimenu('Parent',H MEN, 
[ca liback 7) iteratien method fon return', .. 
"Label', ‘Return to Begining', ; 
'Tag', 'JANRAD OptionsSubuimenul'); 
Ceeeumenuy Parent 7 MEN, 
"Callback"; "iteration method fcn delta input’, 
"Label', 'Change Input Parameters', 
‘Tag’, 'Subuimenul') ; 
uimenu('Parent',H MEN, 
"Callback', Hiner SEion emethod fen about’, 
‘'Label', 'About Janrad 98 oa 
"Separator','on', 
'Tag', 'Subuimenul') ; 
ULCOnEEoON, Parent | Hl METH, 
'Units', 'normalized', 
Mpackqroundeoler Pa Os 752941 0.752941 0.752941), 
PEOnt Ss) Zeural2 - 
"FontWeight', 'bold', 
"ROSadtion:, (O20 do014 0.664 0.333333 020583333), 
"String', 'Choose Iteration Method', 
"Style', ‘text’, : 
“lag, Seabiclexcl* ); 
HONS =surcontrol Parent’ 2no lt METH, 
LCalblpack., steration method fengh ni';... 
'Value',1l, 
hUnats ‘normalized’, 
oe eer oundcoler Or 752941 0. 752941 O.7aZ2941)), 
'Fontsize',;12; 
eos melon 710. Ovoo7 | 4 0.760 -0,552145 -0..0533333), 
rotrang', Nowiteration’, 
roeyle’, ‘radiopuctioen * ; 
'Tag', 'Radiobuttonl'); 
H AS = pepe ponent. He METH, 
‘Callback’, ‘iteration method fcn h as’ 
A059) ote ese areeetr 
"BackgroundColor', [0. 752941 O.752941. OF 752941) 
PROnNtSI Ze’. 12, 
weGsitlon -0.0765714 0,69533370.332143 0.0553333), 
"String’, 'Airspeed', 
‘ocvle >, LadioOnbucton , 
‘lag’, Radiobutton2 ') ; 
H AL = Be Cuenca: «HEE METH 
ECatlbachk ia accrativon Mechod fen naad. | 
"Units', 'normalized', 
BackgroundColor’, (0. 752941 0.752941 0.752941], 
‘Wontsize’,12, 
PPOsit tron, (0.0785714 0.605333 0.332143. 020533333), 
Poerang ; Altitude’, 
'Style', 'radiobutton', 
Plage 7) Ragiobuttons |) ; 
Ho GWo= Uicontroel( Parent’,H tT METH, 
"Callback', ‘iteration method fcn h gw',... 
"Units', 'normalized', 
ac Groundeolor: moi Ore 752941 O275294) 70. 752941), 


Q 
I 


Fh 
II 
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 PODEST Ze ez 
‘Position’ 40-7 0065/1450. 52266760. 5324s e0n Uo sss soe 
'String', 'Gross Weight', 
'Style','radiobutton', 
'Tag', 'Radiobutton4'); 

H BT = pipe oo = Ue eremngtpg img oe 
'Gallback 7 iteration method@renen bt | 
"Units',; normalized, 
Background@oileee alee: 752941 O23 752941 02752941), 
'Fonteize jaz, 
"POsSition: -0.073857147 02.44 0.332143. 0.0533333)]- 
USELING.,- blade Twist: 
"St ville") raditonpucton. , 
'Tag’,'Radiebutvcon> ).; 

H BIR = CIRO a eRM eee mn 
Heal back jaa iteraderon Meehioa shehehwoti., .. 
Units" ,“normalized', 
Tecgecmenies: 7 LO 752941 0.752941 0.752941), 
SPOntsSi ze" 12: 
'Position', [0.0785714 O2s6-0- 332143. 0. USSss333)- 
'String', 'Blade Taper Ratio', 
POeCVLe. LradlLOouston, ; 
'Tag', 'Radiobutton6'") ; 

H SOT = Bene Ot ecee eet? METH, 
Callback’, iteration method fen h sot", 
PUnLES = normalized; 
pbackgroundColor’ 710. 752941 0.75294 m0r7 52941 |e 
PEONtS1 Ze" — 127 aa 
'Position', [0.0785714 0.277333 0.332143 0m@53ee33), 
mirting yo sotaLrenoL Llaper, 
"Style', 'radiobutton', 
'Tag', 'Radiobutton7') ; 

fevon — UConn rol Parent ,H TT METH, 
Callback”. ‘iteration method fen beweay 
“Units’, ‘normalized’, 
"BackgroundColor', [0. 752941 0.752941 0.752941], 
uBentSize*, 12 7mee 
SPOS Hone s| Oe OT vor la G.205355 0.552 L4300,05s5530), 
'String', 'Wing Span Area', 
Potvile* > radlobutton: 4 
'Tag', 'Radiobutton8') ; 

fet CONer ol ( Parent 7 Hlil MELH, 
Units 2 normalized’, . 
ecocounnenlec! 7 LOS 752941 0.752941 0.75294 ie 
SHOMeoDZe 7 i2, 
TRanitieiigine | Mere, 
EPOSLUlomMe 10. 402143 0.866667 0. 439206 .0 0533234) 
roaring’, Analysis Status Box"; 
Poevyle*,. text", : 
‘Tag, Statbiclext.. ); 

HeSTALUS = wicontrol( Parenc 7H Il sMETH, 
‘UNnLEs "normalized", 
‘Bachguemceeloe: Or, 752941 0.752941 0.752941), 
"Position',[(0.482143 0.186667 Om439286 0.64), 
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H STATUSI SUreomerOlmnaLeme i Ll METH, 


Hest 2 aircon rol Parent.) iH lh eMELH, 


Boieyte ., 'CExt 
Bneritsize!, 12, 
Lronewelgnt 4. bold’, 


eee oument, center’, -- 


“Selene fuer ; 
"Tag", 'StaticText2' )3 


‘Untes ) "normals zed', 


eee reeoundcolon: ZO? 752941 0. 75294170. 752941), 
“‘Pesation , (0.466163 0.411764.0.436029 0.217195], 


'Styiies, "text", 
“penesize',12, 
'FontWeight','bold', 


"HorizontalAlignment', 'center',... 


Roe EEN: 7° * A 
‘Taq ; "ctabiecrexts” ); 


"Unies, normalized’, 


spachenoundcolen mus 752941 0.752941 0.752941], 
"eesttien , (0.486szZ29 O- 191554 0,435398l 02205128), 


cstyle” , cext', 
NEOnMEOlZe' > 2, 
"FontWeight', *bold', 


"HorizontalAlignment', 'center',... 


PSEGIng 7° : 
'Tag', 'StaticText4'); 


HUBKe=sUlecontrol (" ParenG ser MEG 


‘Un@es *, “normalized, 


‘Callback’, ‘iteration method fcn back’, 


‘Honest ze 2, 
MHOneEWelguic , .DOld* ; 


PeOoltlenwmslee07 G7 057 O07 064)70.1785 7/1 0.072], 


string 7° << Back’, 
Wage eushbutceen!'); 


H GO = We Glee l mea renc yk IT METH, 


HORU Eee enero Parent: ) iT MEP, 
(Gabliback 7, 1teratien method) fen sinterrupe’, 


hUnaits  —' Mormalizeda’, 


"Callback', ‘global PERF _ OUTPUT; REGIME=0;iteration method fcn 


UBOntLS1 ze 412, 


"FontWeight', 'bold', 


BeC=te1 cna 0.501706 0.0613333 0.178571.0.072), 


‘String’ ;, Analyze‘, 
'Tag', 'Pushbutton2*); 


‘Unies , normalized', 
rFOHES1zZe4, 12, 
WHonewer ght", *bold', 


 feomiedone, (0.526071 0.0613333 0.178571 0.072), 


TSELINng  -INnteLeupte', 
"Enable! Soft", 
‘Tag. EPushbutten3s”* ) ; 


Hehe oteeetconterol( Parent 7 i IT METH, 


pCallback |, iteratien method fen resume’ ,... 
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anal = 


Lastest, 
'FontSize 


*FontWeight’; 


'Position' 7 [OF Jo535/_0) 0ols333 
, 'Resume' 
pon a) Galen 

', 'Pushbutton4'); 
Hold METH, 


‘Strings 
‘Enable’, 
'Tag 


‘Wrnats * 


WICONEEOL@ Parecne.. 
"normalized', 


"normalized', 
pore 


(boLd!, 


O24 7 ear tO. O2)\y, 


eee clon nlLOs0678571- 0. 176 OS asso 7 (Ole Ue, 


YSEYLe* 
‘Tag’ 


‘Units 


,'frame' 

, 'Framel') ; 
UleEONCEOUm@eatent ~HelleMETH, 
oR ried 


eee ronnie oilom: 70. 752941 O77 9294 VE05.752941), 
“Position';, (0.476706 0.17670.45 705762]; 
style’? £rame * 
'Tag', 'Frame2"); 


assignin (‘base 


assignin('base' 
assignin('base' 
assignin('base' 
assignin('base' 
assignin('base' 
assignin('base' 
assignin('base' 
assignin('base' 


assignin('base 


assignin('base 
assignin('base 


Uy Henn. 
ease 
pe iy 


f 
f 
f 
f 
f 


f 


, 'H_GO' 
pores sh od) 8 bay 
assignin('base' 


,H_NI); 
t H_ AS) P 
,H_AL); 
'H_ GW',H GW); 
eet ae BT) ; 
Poe Reel R):, 
SOT 7 SOT); 
'H WSA',H_WSA) ; 
,H_GO); 


'H_BK',H_ BK); 


', 'H_RES',H RES); 
',"H MEN',H MEN); 


,H_RUPT) ; 
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APPENDIX M. ITERATION METHOD _FCN.M 


Switchyard Callback function for the iteration_method.m GUI function. 


fuUncCTCILON Tterpacrtonsmethod fen (Action) 


$ Switchyard Callback for iteration method.m 
%$ JANRAD 98 VERSION 4.0 


global H_IT METH H_NI H_AS H AL H GW H_BT H BTR H_SOT H_WSA . 
H STATUS H_STATUS1 H_STATUS2 
H GO H_BK H_RES H RUPT H MEN 
H HIGE H_IT BOX H_ASPECT H_ASPECT EDIT 
S PERF INPUT S USER INPUT S PERF OUTPUT PICK REGIME 


at Nnargin, 

switch Action 

Case ena es 
set (Heh ie Values} 
Serhan varue | 0) 
Set(ivAl Value” 70) 
set (H_GW, 'Value',0) 
set (Heb, values.) 
Sel (iebihy value 710) 
Set(heoOl 7. value.) 
set (H WSA, 'Value',0Q) 
PICK=0; 

Case “Nias: 
Ser CHa Val Wen} 
SEGGQHeAs, value" mil) 
Sel (HEAL Value 70) 
set (H_GW, 'Value',0) 
set (H_BT, 'Value‘, 0) 
SGelteerR,. value , 0} 
set (H_SOT, 'Value', 0) 
set (H_WSA, 'Value',0) 
PICK=1; 

CaSeue hc 
Set (a Nive /altue, 3} 
Seu(Hen., Value’, 0} 
set (ieel, Value", 1) 
set (H_ GW, 'Value',0) 
Set(Hebe, Valuet, 0) 
Seu HeeiR, Value’, 0} 
SCe(deosol, Value’ , 0} 
Sel (HeWoe, value’, 0) 
PICK=2; 

case th gw' 
set (i Ni,’ Value’, 0) 
SEC(Heie. Value’; 0) 
Scllngee, value 70) 


co 
ND 


seu (i GWyayvaluec iy) 
Serle: ,. Valves, 0} 
See(HiebiR, Value 770) 
Sette sO aavalue |, 0) 
set (H_WSA, 'Value', 0) 
PICK=3; 

ease ~h Dt 
Seu (Hei. Value. 7.0) 
ser i(HeAs, value’ ,0} 
Sete yd Lue 770) 
set (H GW, 'Value', 0) 
Sees iaevalue | 1) 
See (i eSih value, 0) 
SCUln sOl, »Value ; 0) 
set (H_WSA, ‘Value’, 0) 
PICK=4; 

eacen abit 
set(H Ni, value 70) 
Seu(Heao, wa lue , 0) 
Seu (heal, value”) 0) 
set (H_GW, 'Value', 0) 
Sem sl ay alue: 710) 
Set(H BTR, Value 71) 
set(H SOT, *Value'’,0)} 
set (H_WSA, ‘Value’',0) 
PICK=5; 

Case = esoe 
set(H NI, Value ,v) 
sev(H AS, Value; 0} 
See(H AL, Value 70} 
Seuviacw, valuc’ , 0) 
Sermxn Bi, value , 0} 
Seu(HeBriR, Value , 9) 
Sec(HeoOn value’) 1) 
Sec(HeWwon, Value", 0) 
PICK=6; 

Case fewsa 
Sette Nil,. value, 0} 
Ssec(H AS, ‘Value', 0) 
SCG (hens value 70) 
set(H GW, 'Value’, 0) 
Semi 8). va Lue, 70) 
Sem(UHebik a value 70) 
Seuuiaoon, value. 0} 
Sec(ti WoA, Value’, 1} 
PICK=7; 

case 'back' 
close (H IT METH) 
performance input 

case ‘anal' 
set (H_GO, 'Enable', 'off'); 
Set Hervlh, Enable’, “on );; 
Setar, onamle 7 Off jy, 
SGC (HeRES, “Emable*;, orf); 


set(H MEN, Bnaptes, Off.) 
if get (haa alwe |) —— 
Perf 
elseif Get iias . Value.) ——1 
iteration parameters 
Set (HPL EebOx ote nga, Seek oeEED 
ellsSereedet Give, avalue.)—— 
PEC EAtTOn, PaLanebers 
set (H HIGE, ‘Enable’, ‘on") 
Set (Heh) (BOs as sth lig. (eA Liviu Dis) 
elisei tage rwieol,, | vde wi 
iteration parameters 
set (Het) BOX, String’, “GROSo WEIGHT ') 
elseif get (HB BI, “Value’)—-1 
iteration parameters 
set (Hell BOX, string. BEADE iwisi) 
elseif get(H BIR, “Value® )=-= 
PECL at SoMmparalieters 
set(ieis BOx String: , “BLADE TAPER  RATEO™} 
eliseir gecuiHeael, Valtey) —— 
Fteration parameters 
Set (HaldbOx SiO) 4 sod eOr sly CEs) 
elseif get(H WSA, ‘Value*’) ==] 
iteration parameters 
set (H IT BOx, String’, "WING SPAN AREA”) 
Sethe ASLEC?, Eaaple’ >) on: } 
Se€t(H ASPECT EDIT, "Enable", ‘on, ) 
end 
case ‘interrupt' 
sel (iGO Enavile | OLf. ); 
sel (HeRUE Ma Ehavle 7 OLE’), 
Set (HeBr, Enable’, off); 
Set (HeRES, Enable on); 
Set (HEMEN boa le’ > on}; 
uiwait; 
case 'resume' 
SEG HeCO,aunalle 7) /OLr.) ; 
sec (H ORUER Enable )on”); 
Seti Hub enaowe | Orr")! ; 
SEE(Henio, Cnacle , OLE”); 
Seti MEN enable,  Oft.), 
ulresume; 
case ‘guit' 
Gu ea gu 
case ‘return' 
close (H_ IT METH) 
janrad98 
Cases Geltaginpuc. 
close (Hal) METH) 
performance input 
case ‘about' 
about _janrad 
end 
end 
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APPENDIX N. ITERATION _PARAMETERS.M 


This file creates GUI to enter iteration parameters. It 1s called by the Switchyard 


Callback function iteration_method_fcn.m. 


function iteration parameters (} 


% GUI window to enter iterative steps. 
% JANRAD 98 VERSION 4.0 


% This is the machine-generated representation of a Handle Graphics 

object 

% and its children. Note that handle values may change when these 

objects 

%$ are re-created. This may cause problems with any callbacks written to 
depend on the value of the handle at the time the object was saved. 


% 

3 

%$ To reopen this object, just type the name of the M-file at the MATLAB 
% prompt. The M-file and its associated MAT-file must be on your path. 


ROoadeiterdt bone baLdmeters 
global H_IP H_HIGE H_IT BOX H_ ASPECT H ASPECT EDIT H MEN 


Hele figure ( Units, normalized”, 
Meolor., |Oes 035.0. 38)., 
*Colormap',matO, 
‘Name’, ‘Iteration Parameters', 
BNUMpDer rat Le" Off. 
"PoINTERShapeCData',matl, ‘ 
PPesttion , (0.04375, 0-089 5633 0.G7o- 0.78125), 
UU oa) alia pa 

pee urmenu( Pareniu, HEP, 
"Label', 'JANRAD Options', 
‘Tag’, uimenud”); 

c = uimenu('Parent',b, 
(Callback Viteracion parameters fen Guat 7 
"Label','Quit JANRAD', 
'Tag', 'JANRAD OptionsSubuimenul') ; 

c = uimenu('Parent',b, 
Peaiback weeterdtion palameters fen retwin , 
"Label', 'Return to Begining’, : 
'Tag', 'JANRAD OptionsSubuimenul'); 


c = uimenu('Parent',b, 
‘Callback iteration parameters fen delta iiput "7 
‘Label', 'Change Input Parameters', 

‘lad. SubuLMenuwL | > 

c = uimenu('Parent',b, 


"Callback', 'about janrad', 
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PLhabel’, "“Abovteoanmrads98.....', 
"Separator','on', 
'Tag', 'Subuimenul'); 
Pe sUdcOnrE boll sbanent jo LP, 
‘'Units', 'normalized', ‘ 
Mee ectndcolor: ehO.. 752941 O27 o2 940752941), 
HROMEGLZC 712, 
‘'FontWeight', 'bold', 
PeoStewvon ./0.26007 14 OF 6860, 476786 0.0533333)], 
"String', 'Performance Analysis' 
PStyle 7 bext", : 
[bade Statirelrext Ll”) +} 
Hence we su teOntron( Parent 7H EE, 
‘Units’, 'normalized', 
"FOmMeo) Ze, Lz, 
"Fontwelgnit “mold; 7... 
"Poesat LOnw.10.2600/14 OnecZovOn4avb7e80 0.0533335), 
POoCVYLG ay, EeKE ‘ 
Lagden, ctatLierext)”’ ) ; 
De— UbeGMErol| Parent ,Helk, 


'Units', Pioemaic cd’. 
bPOSitiLon., (0.2625 0. 705778 O30 1 14320 ,/05553533 15 
ocr net chart Prerationvat.:., 


wSEVLe: 7: Lexe., : 
"'Tag', 'StaticTextl') ; 
Po -siecon trol earent "7, Hole, 
'Units', 'normalized', 
Peackoroundeae es = tL ale al 1], 
Ub@sigeron - (Oe625.0.704 0,107143 0.0533333), 
"Style','edit', 


"Callback','global MINUM;MINUM=str2num(get(gcbo,''String''));' 


lag, EQLEText ll); 
Pe eee oOnur on (Parent, alr, 
PUES: > normalized’, ««. 
Brest tan: PLO. 26250): 634667 0.357123 020532233) 
ESELInNG | End iteration at <", 
Foyle 4 text; : 
‘hag, staticlextl*) ; 
Pe eUrcone rola Parent. ,; Hale, 
"Units', 'normalized', 
rene heroundcoler Pe (el ah ll, : 
“MPOsaeron 7 (O.02 5 (0.634667 0. 107 1437 030535333] 
‘SGVLe", “ede, 


"Callback', 'global MAXUM;MAXUM=str2num(get(gcbo, ''String''));' 


Lade ys, Bae lrextl *); 
b —sutecene rol ( Parent Hl iP; 
'Units', 'normalized', 
"Pos tevon  - (0.26025 (0.563556, 0.357143. 0 0n33333), 
Sering-, iteration, Interval, 3", 
“style”, ’ text", ‘ 
‘lag: 7. Stabierext |’ > 
Pa -ewheonl rou Parent, i  lE, 
"UNEES*, "nOormalized*, 
mae Srcundcolor Pid (ie ae 1], 
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'Position', h0e625. 0.562667 -Oa00 143 (ORs oeeecc 
‘Style edit 


"Callback', 'global INTER; INTER-=str2num(get(gqebo, “sering . 07, = 


'Tag';, Eailerextin)c; 
H ASPECT = uicontroli@rarent. , HIP, 
"Units', 'normalized', 
Ppositiont miLOs2o2s 0. 492444 02355357. 02 05555331, 
String’; (Aspect. Ration: | 
Style’ text | 
‘Enable’ >= ort; 
‘Tag, Stavrlerexti )- 
eS eect  EDit = suiCcont role Pagent 7H. iP, 
'Units' ‘normalized', 
eeraroud eno rag leh 1], 
EPOSI ETON iL 0e.o25) 0. 49066) 70. O04 3. 0 205s3555 15 
rstyle',*edice*, 
“Enable, 'oLrk*, 
"Callback', 'global AR;AR=str2num(get(gcbo,''String''));' 
PTagde bale hextl *):; 
BeotGhe—suicontrol ('Parent 7 HOIeP, 
‘Units’, *normalized', 
"BackgroundColor',[0. 752941 0.7 5294126. 52941), 
"Enable. OLt “+ 
EPOS Teton a) Or 2c Ocoe -4205 5G 0.550929 6.053533 3 ]/; 
Pocking ) IncludesHIiGE Calculations ?*, 
roeyle™, checkbox’, 
pValue * 0-9. 
PCallia Cla. lf 
get (gcbo, ''Value'!')==1, REGIME=1; else, REGIME=0;,end', 
'Tag', 'Checkbox1'"); 
WECOMELOIn we bAarene., ti TP; 
HUnles , nNornateced , 
Bed) pac tee Lacon Pardtecers | fem dich. 
PEOMES1LZe 7 9h2,  ... 
FOoneweiqnt , bold’, 
yeoscitcion 7 10.260714 0.245333 0.196429 0:.0986067]), 
'String', '<< BACK’, 
had, Fushouttonl') ; 
PeaeuLcOnerol(=Parent*,H IP, 
MUnbes. normalized’, 
pea PEDACk teeta ton Parameters fen anal", 
Peontoize zs 
HOMEWCIGNE GF. DOld', 
“MPOSition 7 10.54265/ 0.245333 0.196429 0:0986667], 
'String', 'Analyze >>', 
Lag.) cushbutton!') > 
Ce srcont cols Parenc ,H IP, 
"Units', 'normalized', 
Mpackoroundeclon LO: 752941 0.752941 0.752941], 
"FOHneotlzey, L2- 
Tae) Deo : 
"Position', [0.180357 0.088 0.646429 0.106667], 
'String', 'Warning - Exessive Iteration Limits May Increase 
Processing Times!', 


b 
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It 


rotydle , text’, , 

‘tog, ‘StaticTextl'); 

Meconcrol( Parent 7H IP, 

Unmet S eee 

ap eeaeiont LOn ty 67 386.0. 0746667 ON Sys a pSyoyy ion Cre ies) SPC y Cy 
'Style', 'frame' 

seagate amel'); 

Uaeomenei| Parent ‘7H ir, 

‘Unhres. eae reat, 

poention: pO weooo) 1 Ur, 810667 0.498214 0.146667], 
'Style', 'frame' 

‘Tag’, 'Frame2'); 
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APPENDIX O. ITERATION_PARAMETERS_ FCN.M 
Switchyard Callback function for iteration_parameters.m GUI function. 


Lunction 1teratton paramevers fon eeeten, 


% Switchyard Callback for iteration parameters.m 
% JANRAD 98 VERSION 4.0 


global H_IT_ METH H_IP H NI H AS H AL H GW H BT H BTR H SOT H WSA H HIGE 


H_GO H_BK H_RES H_RUPT H MEN H STATUS H STATUS1 H_STATUS2... 
S USER INPUT PICK 
MINUM MAXUM INTER REGIME 


if nargin, 

switch Action 

case ‘back' 
SG (Heh, Enable’) on j7 
SceieGe -senable 7 one), 
Scllehvil, sEndble | Ort )7, 
Set(ierss, Enable”; oft); 
set (REMEN, Enabte 7 on 
close (HEIE) 

case ‘anal' 
Se Cee Enable oti. 
SseliieGO, Enable’) "otk: }; 
Sem ioRUPT, Enable { en’); 
setiH RES, Enable @7 Oft )y 
SeeqH MEN, Enable’) off"); 
eloseqHelr)) 
Perf 

Case ‘quit’ 
quae qui 

case ‘return' 
janradg98 
elose Hite) 
close (H_IT_ METH) 

edse del taeinput 
performance input 
eroce (Her) 
close Ope sLab isis) 

case ‘about' 
about janrad 

end 

end 
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APPENDIX P. OPTIONS.M 


This file creates the GUI to select additional analysis methods and print input and 


output files saved from the performance output window. 


function options () 


% GUI window to Select user options at end of performance routine. 
% JANRAD 98 VERSION 4.0 


% This is the machine-generated representation of a Handle Graphics 
object 

%$ and its children. Note that handle values may change when these 
objects 

are re-created. This may cause problems with any callbacks written to 
depend on the value of the handle at the time the object was saved. 


To reopen this object, just type the name of the M-file at the MATLAB 
prompt. The M-file and its associated MAT-file must be on your path. 


oP © o& &® oe 


load options 


global H_OPTIONS H_PSCA H_ PRDA H CIM H_CID H RTB H_EJANRAD NAME 
H datain H_dataout H_ vecdata 
Bepwicin He printout Heprintvec 
Bee npuet ile HW Sutpult rules vyecrr lc 
Hecmeckl HT ieheckZ checks 


H OPTIONS = figure('Units', 'normalized', 
VEOlLOG: ,( 0.6 0. cord |, 
‘Colormap',mat0O, 

"Name', 'Options', 

"NumberTitle','off', 

"PointerShapeCData',matl, ... 

"Position"; (-O7003T2Z5e0" 0625 0.954685" 0.8625], 
lah, PLO lees 

Peeve nue Parente, HOT LONG, 

"Label', 'JANRAD Options', 

Lag, ULmenul” ); 

uimenu('Parent',b, 

peal loaew, . OPelons LeCnequit , 

‘"Label', 'Quit JANRAD', 3 

'Tag', 'JANRAD OptionsSubuimenul1') ; 

uimenu('Parent',b, 

(Cal lpacke ye Opelons: ECh ect urn jc. 

"Label", "Return to Begining’, ... 

"Tag', 'JANRAD OptionsSubuimenul') ; 

uimenu('Parent',b, 

-Callpach |) ©petons Ech del taginpue jc0. 


e) 
H 


Q 
Il 


Q 
II 


2 


'Label', 'Change Input Parameters', 

'Tag', 'Subuimenul'); 

uimenu('Parent',b, 

'Callback', ‘options fcn about', 

‘'Label', 'About Janrad 98 ...', 
"Separator','on', 

'Tag', 'Subuimenul'); 

ULeCnEEOly Parent 7HrOPLIONS,.../- 

‘Units', 'normalized', 

npaeRoLOundcolor miro 752941 Oe S294eto. 752941), 
[POmeaEzZe: 71.6), 

"FontWeight', 'bold', 

"POSlelon , (0.07 1565. 0.006 0.376531 0.109333], 
‘Strang: , oclect Option’, 

"Style", ' texte, ‘ 

'Tag', 'StaticTextl'); 

H PSCA = uicontrol('Parent',H OPTIONS, 


Q 
II 


o> 
HI 


'Value',1, : 
feat loaciie LODE Ons EChah psca' 7... 
ntti s Mecndliced’. 


epacrarcundcolor' Pon) ( Oh 752941 O72 79294170. 752941), 
Eeocit2om’ , |O.0715631 (0.927 0437893 1.0. 0533333], 
'String', 'Perform Stabilty & Control Analysis', 
tseyvile 7. radi obutton , 
'Tag', 'Radiobutton3', 
"Value',1); 

H PRDA = ec aay CR B_OPTIONS, 
[Carirbach  Cpciens Tenenaprda’ ; 
tUnits', ‘neormalizea’, . 
Wpoeneroundcol or aOR 752941 O52 Ode On oc ai. 
BPOSttlon 7 p00 7 eoOoleO. O52 "0-37 003) 0,0533553)]- 
'String', 'Perform Rotor Dynamics Analysis', 
Potvlew, Ladiopuccon ; 
'Tag', 'Radiobutton4'); 

H CIM = ae eas PIONS? o:2' 
pea llback, ,Opelons fen h Pe TG sate, 
‘Units' 'normalized', 
eee Geen LO 752941 0.752941 09752944) ae 
Position’ 7 (05071563 1 e@m541323m0me7 S531 0.0533333],; 
'String', 'Change Iteration Method', 
“Style; radlobutton", 
'Tag', 'Radiobuttonl1'); 

st (CM) Pease e eee TONS 
iGalilback | options fen hicid’, 
LUNES. 'normalized', av 
Spee to roandco lc EP LOs 752941 0. 75294100. 752941, 
Positron 10.0 715631 0.453333 0.376531. 10505305551, 
'String', 'Change Input Data', 
‘seyle", radlobucton , 
'Tag', RadiobuttonZ'); 

Heeb = ULcontrol( Parent 7H OPTIONS, 
(callback 7) options fen h Ttb:, 
TUnLES |oonGtinaliazed', .. 
TpaGhgpouneeoile c yaleO:. 152941 027752941 0.752941], 
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‘Position ",. |0\ O7isseeie 0. 365333507 378531 20, 0533334)5 
‘String,’ Returnseo Begining 
"Style. “Gadivoburton., 
'Tag’, ‘Radzobuttona )); 
H EJANRAD = ppeninic ia Magenta oo Cam 
'Gallback: (). OputoOnomecnahwe s cimisacha, 
‘Units; normalazea’, 
CBaclerommeeior: a Oe 752941 0.792941 0.752941), 
‘Position :, [0207 15631)0-2/4667 70. 376531) 0-050 422). 
‘String’, ‘Exit JANRAD', 
Losey le.,  GacdloObUecwzon , 
‘lag., RadLobucvons. )- 
De= sua CONtEO!( Parene 1 OPErONs, 
‘Units', 'normalized', 
eae GrCundccion: LOR 752941 Oe792941 0.752941), 
'FontSize',16, 
VHontWelgnt-; bold", 
‘Position’, (0.546139 0.805333 0.376643 U-.106667)- 
Poeblng., Print selection , 
bocyle-," EExE" | so... 
plag 7 StaklereaxEz ); 
Bee eiacan = Ulcontrol( Parents, Hf OPTIONS, 
‘Units, "normalized", 
ie eto ee F (ial 752941 0. 7229410, 752941, 
*“Position’, [0.545548 0.71644 0.169906 0.0527903), 
tlrng 7g Prince input file ye, : 
‘Style', 'checkbox', 
plag>s. Gheckboxi, 
‘Value', 0); 
Peupucil le — Uiconcroh|) Parente ,H OPTIONS, 
‘Units', 'normalized', 
EPOntol ze, 12, 
Baa belie od Lak a RM a oi 0.752941 0.752941), 
mote eG sc, : 
peer [0. 748209 O27 1644 071698908 0.0527 7503)- 
bouvyle text’, : 
PHOrizontalAlignment’,'left", 
Phage ye obaclrelexto.) ; 
BEpLINtOUt = Uiconerol( Parent 7H OPTION, 
‘Units', 'normalized', 
Wack GESuNdColLer: Pe ig OMA 752941 07752941 0, 75254lie 
Position 7 (0.545046 O.626959 "07 l69905 00527 20c1- 
Fostring 7 Print Cutpue Pile 2", 
'Style', 'checkbox', 
Hlag’, Checkbox’, 
"Value',0); 
pececeet ttle = uicontrol('Parent',H OPTIONS, 
“UnLES* ; “normalized” , 
Reh aLouniColor: 5 [ele 752941 O,752941 0.792941), 
'"FontSize' wee 
Y Sie aa ovess - 
postion [0. 748209 OTe28 959307 16990s 070527 7dan, 
“SCvyle a text 
"HorizontalAlignment','left', 


101 


bbage, ccaticTexts we 
Rep Eeievece= Ulconcerol({ Parent’,H OPTIONS, 
‘Untes', ‘normalized’, 
ieee eoundccllor' peu. 752941 Ore7 5294120. 2941 ee 
PPosttion 710.94 05495 0.55060953m0m1 69908 0.0527903], 
PSebiag 7 Print Matrix & Vector File :*', 
'Style', 'checkbox', 
'Tag', 'Checkboxl', 


'Value',0); 
H vecfile = uicontrol(‘Parent‘,H OPTIONS, 
‘Units’, 'normalized', 


Ee eroundcolon: mis 752941 0.752941 0.752941], 
WHOneoIZe.” 2, 2. 

‘Position’, [0.748209 besso9%53 0.169905 0.052/790s1" 
PotLanGg 4g yes 

'Style', 'text', 

"HorizontalAlignment', 'left', 

‘Taig’, 'StaticText9*}); 

Due onc Lol ararent 7H OPTIONS, 
peailiback |) Opttons FCn print «, 
‘Unate 7 normalized", 

“ROneoi ze’, 12, 
"Fontweight’, “bola*, .#: 
"Position', [0.595104 0.402667 0.288136 0.072], 
‘Stbingd: , senda to Printer’, 
‘Lage, Eushbuttonz~ ) ; 
bys Ureontror( Parent, i OPTIONS, 
"Unwes | normalized", 
een (0.599774. 0. 362667 02435028 02570667 |, 
'Style', 'frame' 
'Tag', 'Framel'); 
be] uteontrol ( Parent’ 7H OPTIONS, 
‘Units', 'normalized', 
Teed ction! rlO-O4eatS 77 0. 0826667 0544957650 co0GGviL, 
"Style', 'frame' 
‘Tag', 'Frame2'); 
b = uicontrol(*Parent*,H OPTIONS, 
"Callback', ‘options fen back', 
"“Units;)’, ‘normalized', 
“HontSize'", 12, 
'FontWeight', 'bold', 
“Position ,|0.0809793 0.112 0.178908 0.088], 
BString., << Back, 
‘Tag, ‘Pushbutton! *) ; 

b = uicontrol(*Parent*,H OPTIONS, 

‘Gaiback, ; Voperons £Cn Cont , 

‘Units’, 'normalized', 

“Ponto. ze", 12, 

‘FontWeight', 'bold', 

PPosttaon 7 loezcocoe O.112 0.177024 0.088], 
Poe ring), COntinue >>", 

"Tag’, Pushbutton! *) ; 


MEgqese i cheer, Value*)— 
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set (H printin, Envopilew oft") 
end 
HE get (Hh cneck2 ss Vaiie = 

set (H printout,’ Enable “7 ort t) 
end 
cf Get (heehee s Vatieu = 

Set (Ho print vec, Ewaolen ole) 
end 


assignin( ‘base’, 'H PSCA'",H PSCA); 

aSSiGnin ease oeenmonO a ee RDA i, 
assignin('base','H CIM',H CIM); 
assignin('base','H CID',H CID); 
assignin('base','H RTB',H RTB); 
assignin('base', 'H_EJANRAD',H EJANRAD) ; 
asSagnint base |, Me plintim , fH primein ); 
assignin('base', ‘H printout',H printout); 
asain ene base |, H printvee’,H printvec); 
aeeagnin( base =, "H inputftile’,H imputfile) ; 
aacignin( base” , "H outputfile*,H outputfile); 
as Siqmint base se ohevect le’, H veectite); 
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APPENDIX Q. OPTIONS FCN.M 


Switchyard Callback function for options.m GUI function. 


LTunetion Options fcn (Action) 


% Switchyard Callback function for options.m 
% JANRAD 98 VERSION 4.0 


global H_OPTIONS H_PSCA H_PRDA H_CIM H_CID H_RTB H_EJANRAD 
H PEIMC Ip Hapraneeue Hoprintveec. - 
JEP Re eet NAME oo MAtTReaveC Hivecti le prince cemnps 


eendl—get (WH PoCA, Value’); 
Gond2—get (H PRDA,' Value’); 
Gends—get (H CIM, Value"); 
cond4=get (H CID, 'Value'); 
Cendo-get {(H RTB, ‘Value*); 
cond6=get (H_EJANRAD, 'Value'); 


if nargin 
switch Action 
Easeu up Sca | 
set (HoPoCA,) Value’; i) 
set (H_PRDA, 'Value', 0) 
set (H CIM, 'Value', 0) 
Sel eel, value ).0) 
Sec (HeRIB | Value 7.0) 
set (H_EJANRAD, "Value', 0) 
Gase “hoprda’ 
set (H_ PSCA, ‘Value’, 0} 
Seu (Hee RDA value...) 
set (H CIM, Value ,0)} 
see (HeClp, | Value e)) 
sevC(HeRIB, Value 7.0) 
set (H EJANRAD, 'Value',0) 
ease eelin: 
set (H_PSCA, 'Value', 0) 
set (H PRDA, ’Value', 0) 
See Haein Value>, 1} 
Set (HiCrp, |\Value*, 0} 
Sec (HORTB, value 7.0) 
set (H_ EJANRAD, 'Value', 0) 
case herd: 
set (H_PSCA, 'Value',0) 
set (H PRDA, 'Value’, 0) 
Seu (i Cin value 77-0) 
Ssecu(HecEp value’, 1) 
Seu(HoRIB, value; 0) 
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set (H_EJANRAD, 'Value’, 0) 
exsist) i)ai Jere}: 
set (H_PSCA, 'Value’, 0) 
set ({H PRDA, ‘Value’, 9) 
set (H CIM, 'Value’,0)} 
set (H CID, “Value”, 0) 
Seq (Henne value ; 1) 
set (H_EJANRAD, 'Value',0) 
Edoee Mme dialitdc | 
set (H_PSCA, 'Value', 0) 
set (H_ PRDA, 'Value', 0) 
set (H_CIM, 'Value', 0) 
See (Heel Dry Valve 77.0} 
set (H RTB, 'Value',0) 
set (H_EJANRAD, 'Value', 1) 
case 'back' 
close (H_ OPTIONS) 
performance output 
case ‘print' 
tedet (He printin, Value’ }——1, 
eval GCOpyab re Bnew ctemD lgot das) 
delete prine temp 
end 
Hi get (ieprrnecout, Value )==!, 
ehieidl | UY deen? “jeves joie eeniiewh = akon th 3) 
delet eupri nes empl 
end 
Bide rie ee iibnewecy oVcehle 9) —— 
evclitieekedau, decim Vecti le, String’) )}, 
diciay spilt seempZ 
diary off 
delete print temp2 
Gracy prinien cempZ 
Pepsi, Vi, theta, betat, alpha, Tpsi;9Mpsi, DMpsi, a1, dM, 
aD, cblade, CL, CD, 
diary off 
evo copy, Do, printstempe, Iptl*}) 
delece print. tenp2 
end 
case ‘'cont' 
if condl== 
Scabi lily vana control 
elseif cond2== 
rotor dynamics 
elseif cond3== 
eloses (HeOLTIONS) 
S PERF _INPUT.Vinf=S_PERF INPUT. Vinf/1.68894444; 
Sobber INPUT. CWISt=-SUEERP INPUT. twist ov. 3; 
SeeSRE NEO T. thetao—Serhar INPUT. thetac’ 57.37 
iteration_method 
elseif cond4== 
close (H_ OPTIONS) 
pentormance Input 
elseif cond5== 
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close (H OPTIONS) 
janrad98 
elseif cond6== 
eWelies (etal 
else, 
error('Something is wrong in Options Switchyard Callback 
Puneet von) 
end 
case FECUrn 
close all 
janradg8s 
CaASemeeetta Inpucy 
elose (HeOETICNS} 
performance input 
case ‘quit' 
quit gua 
case ‘about' 
about janrad 
end 
end 
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APPENDIX R. STABILITY_AND_CONTROL.M 


This file creates figure window indicating the stability and control functions have 


not been incorporated. 


function stability and control () 


$ GUI Window to notify user that this module not installed. 
% JANRAD 98 VERSION 4.0 


fonts 15 the machine-generatea mepresencation Of a Handle Grapes 
object 

$ and its children. Note that handle values may change when these 
objects 

$ are re-created. This may cause problems with any callbacks written to 
depend on the value of the handle at the time the object was saved. 


To reopen this object, Just type the name of the M-file at the MATLAB 
prompt. The M-file and its associated MAT-file must be on your path. 


oP 0 of of 


ligagd stability and control 


a = figure('Units','normalized', 
VCOlor (0.30. 0.8 0.3), 
"Colormap',mat0O, 
"MenuBar', 'none' 
'Name' , "Stability ar Control Not Installed', 
TNumberTitle*,"oLk*, 
mee Carshapeepacat mati: ‘ 
Pee@sttr1on';,[0. £90625 0.383333 0.446875 0.34375) - 
Plage, bagi); 
b = uicontrol('Parent',a, 
'Units', 'normalized', 
'Galliback., close(ger)*, 
PHonesize., 14, 
PphOontwerghte |, "bold',. -..« 
EPoOstETOn On s0g112.0¢ 109091 O22Z02797 0. 18ree |), 
PSering' | OK, 
whage, Fushbucton|*) ; 
We Urecontrol (*Parent",a, 
‘Unmts*, "normalized', .. 
iP chovounoconer: BuO: 752941 0.752941 0.752941), 
VE OMG ST Ze 0l2- 
“POHEWeCrGhe ~oDOkd', ..:. 
"POostetom, (0. 06s9161l O.515152 0.911169 07345455], 
"String', 'The Stability and Control Function is not yet Avalilable 
in JANRADS8'!, 
TSeVle , text’, : 
"lag: | ovaererext |’) > 
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pa=eutcont rol ('Pardntwijna, 
'Units', 'normalized', 
PEaengroundeotor’ 7710-75294 0.752941 0.752941), 
LROn EST Ze" ..12, 
"FontwWweight 'jaybold’, 
PPOs heron Oss 2 Moo 0.359576 0.318182 0.127273], 
Potuingd, fu sORKY.  , 
'Style', 'text', : 
'Tag', 'StaticText2'); 

b = uicontrol('Parent',a, 
'Units', 'normalized', 
"BackgroundColor', [0. 752941 On o294)1 0.752941), 
'Peosttron , (0.015 900, 0,0363630 0.972028 02933333), 
"Style', 'frame', 
'Tag', 'Framel'); 
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APPENDIX S. ROTOR_DYNAMICS.M 


This file creates figure window indicating the rotor dynamics functions have not 


been incorporated. 


function rotor dynamics () 


%$ GUI Window to notify user that this module not installed. 
% JANRAD 98 VERSION 4.0 


%$ This is the machine-generated representation of a Handle Graphics 
object : 

%$ and its children. Note that handle values may change when these 
objects | 

are re-created. This may cause problems with any callbacks written to 
depend on the value of the handle at the time the object was saved. 


To reopen this object, just type the name of the M-file at the MATLAB 
prompt. The M-file and its associated MAT-file must be on your path. 


oP &® oC CH 


load rotor dynamics 


a = figure('Units', 'normalized', 
ncolLor 710.5 0.650.931], 
'Colormap',mat0O, 
"MenuBar', 'none', 
"Name', 'Rotor Dynamics Not Installed', 
‘NumberTitle't, "oft, 
"PointerShapeCData',matl, ... 
-PoOcuviton . (0.190625 0.3875 0.45625 0.339585], 
tage, Figi:), : 
be — Uuicontnrol (‘'Parent', a, 
‘Units’, ‘normalized', 
"Callback', ‘close(gcf)', 
'FontSize',14, 
'FontWeight', 'bold', 
“POSition!710.389078 0.153374 0.204778 0.184049], 
wisn eal) ole ere 8 cue 
tlag: , Pushbutton! ) ; 
b = uicontrol('Parent',a, 
“Untts. ) normalized’, ~.. 
‘BackanoundColor', (0.752941 0.752941 0.752941), 
MEONG SIZE! 4127 
'FontWeight', 'bold', 
'POSition! woe OsaZI9I6 0.595092 0.849315 0.319018), 
'String', 'The Rotor Dynamics Function is not yet Avalilable in 
JANRADS8', 
‘Sotyle",; vex, 
“tagG’, Statlelextl’ ) 


11] 


be= ua Control (‘Parenme a, 
"UMmacsS 2’ mormalized’, : 
ipdemgmounmacolor', (0./52941 0.752941 0.752941], 
hrontesrze', 12, 
'Feontwevegnt', ‘bold’, Eas 
"Position', [0.334471 0.429448 0.317406 0.122699}, 
tString = ,;,.oORRY!*, 
Toevle. , CexXt ', , 
"Tag’, ‘StaticText2"); 

b = uicontrol('Parent',a, 
"UnweES jo nerMmali 2760", <a. 
"“BackGrounadGolor*, |0./52941,.0./752941 0./5294,L1, «wu... 
"Position -(0.0205479 0,03806748 0.962329 0793865), 
'Style','frame', 
"Tag’', 'Framel'); 
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APPENDIX T. QUIT_QUILM 
This file creates GUI to verify the users intention to quit JANRAD 98. 


PuncE LOM Cu isG eG uime) 


% GUI window to verify the users intention to quit Janrad98. 
%$ JANRAD 98 VERSION 4.0 


$ This is the machine-generated representation of a Handle Graphics 

object 

% and its children. Note that handle values may change when these 

objects 

% are re-created. This may cause problems with any callbacks written to 
depend on the value of the handle at the time the object was saved. 


% 

% 

% To reopen this object, just type the name of the M-file at the MATLAB 
$ prompt. The M-file and its associated MAT-file must be on your path. 


moad quit qui 


a = figure('Units', 'normalized', 
BcCoLor > (0.6.05 650-0), 
"Colormap' Do heae 
"MenuBar','none', 
"Name', '''Quit SANRAD peo am 
‘NumberTitle' moO r » 
"PointerShapeCData',matl, 
'Pestteron. ,)Oeczso9s0 0.433333 0.389063 0.266667], 
plag. rag); 
b = uicontrol('Parent',a, 
‘Units', 'normalized', : 
"callback", “close (get) ', 
“POM od Ze es oe. tare 
PPeGnEWeclGnht oe bold', “2. 
“Posteren ;, (0.188755 0.273438 0.240964 0.15625), 
UStiunge ov NO’ sack 
lage, = Ushpaceon |” )i; 
b = uicontrol('Parent',a, 
Whits’, normalized", 
"Callback', 'close all,clear', 
WROntws1ze" 7 127 
"FontwWeight*, 'beld", 
*"Position'’, [0.566265 0.273438 0.240964 0.15625], 
“SEpUnG us YEO; 
“laGege Eulsheuttonz  ) > 
b = uicontrol('Parent',a, 
"Units', 'normalized', 
‘pact rounacoiee: , 108 752941 OO. 75929410. 7529417) 
"POnEoLZe! Giz, 


Inks 


PHonkeweignte', 'bold', 

Secoitieronn| O.19239 4) O.59375 0.618474 0.289062], 
'String','Do You Really Want to Quit JANRAD 98?', 
'Sstyle*, "text', ‘ 

‘Tag’ @ StatLeText ier, 

WCcontrolit Parent’ ,a, 

“Umies normalized’, ....% 

Backo Loundcolor', (0.752941 0.752941 0.752941], 
UPostelon 7il0.0401606 0.0625 0.907631 0.898438], 
‘Style’, 'frame'’, 

'Tag', 'Framel'); 
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APPENDIX U. TRIM WARNING.M 


This file creates GUI to inform user that the performance routine did not calculate 


a valid solution. 


function trim _warning() 


% GUI window to notify user that conditions will not trim. 
% JANRAD 98 VERSION 4.0 


$ This is the machine-generated representation of a Handle Graphics 
object 

%$ and its children. Note that handle values may change when these 
objects 

$ are re-created. This may cause problems with any callbacks written to 
depend on the value of the handle at the time the object was saved. 


3 
3 
%$ To reopen this object, just type the name of the M-file at the MATLAB 
% prompt. The M-file and its associated MAT-file must be on your path. 


load trim warning 


a = figure('Units', 'normalized', 
Teolorn’ ; (056 0-58) O61, 
"Colormap',mat0O, 
"MenuBar', 'none', 
"Name', 'WARNING!', 
"NumberTitle', 'off', 
*PointerShapeCData',matl, : 
PeOsition® -1 0.16437 540.352083 075375. 02404167 ]— 
lads. Baga.) ; 
uicontrol('Parent',a, 
"Units', 'normalized', 
ppaeroroundcoltent Lor 752941 0579294170. 752941 |, 
BONES Ze” , b47 en... 
MPOsStETCh: 0.07 59014 ©. 551546 0 5e255c1) 0. se Ck 
roeiaing' this Conltiguration will not trim '", 
EStyle*;, "text, 
PlaG  » StaticTextl"); 
UPeonc role Parent’ ja, 
‘Units *, *normalized' ; 
pcaliback 7 trim warning fen ; 
"FontSize',14, 
VEOMEWerGht ’; sao ld 7-2. 
"POsiTt10On + \L0.s573947 0-092 78035m0.22965!] 0.16494 3)~ 
Touring , Che. 
‘Tag’, Pushbutton! ):; 
b = uicontrol('Parent',a, 
'Units', 'normalized', 


oO 
tt 


log 
H 
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SSehgrounacoloms| 0.992941 0.752941 0.752941), 
wEQemelon 710,07 o90514 0.840206 0.82558] 0.170103), 
seiang ,;mMatZ, 

noovile , text, : 

iTtag ;, "seatrerext ly"; 

uicontrol('Parent',a, 

‘Units', 'normalized', 

peeetgeoundeolor: eal. 752941 Of 752941 0..79294 1), 
“Ponta aze', 2, 

"FontWeight', 'bold', 

"rOsteron slo. 0 7550614 05695876 0.625581 0.104948), 
"String', 'Performance Analysis Routine Terminated!', 
*Style', 'text', : 

‘lad js scatLLelext 1’); 

uLcontrol('Parent",a, 

‘Units. normalized”, : 

'BaekgroundColorm | 0. 752941 G2 7529 40Oe7 529411; 
*Position’, ([0.0348837 0.0463918 0.924419 0.886598], 
Esty le' "frames; 

'Tag', 'Framel'); 
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APPENDIX V. TRIM WARNING FCN.M 


Switchyard Callback for trim_warning.m GUI function. 


function trim warning fcn() 


% Switchyard Callback function for trim _warning.m 
% JANRAD 98 VERSION 4.0 


global H_GO H_RUPT H_ BK H_RES H MEN 


Se ma heeGO pn oi aoe Or vale, 
Be GienUbT  / enable lotta) 
set (H BK, 'Enable','on'); 
fea(dekeos, Enaole', "ori'); 
Set (H MEN, Enable’, ‘on'); 
close(gcf) 
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APPENDIX W. EMPTY_BOXES.M 


This file creates GUI to inform user that all input edit boxes must contain a entry 


to properly execute performance evaluation. 


EUNCELen enptyeboxes () 


% 


GUI window called if Empty edit boxes exist when continue button is 


pressed. 


% 


% 


JANRAD 98 VERSION 4.0 


This is the machine-generated representation of a Handle Graphics 


object 


% 


and its children. Note that handle values may change when these 


objects 


& ® © ® O° 


are re-created. This may cause problems with any callbacks written to 
depend on the value of the handle at the time the object was saved. 


lomreopen this ObjJeCcE, gust type the name on the M-file abeene Mares 
prompt. The M-file and its ass@ciated MAT-file must be on yveur, path. 


load empty boxes 


a 


= figure('Units', 'normalized', 
Peolor';, 10.8) 0.8 O48)7 
*Colormap'’,matO, 
"MenuBar', 'none', 
'Name', 'ERROR', 
eNumberritile*, “off, 
"PointerShapeCData',matl, 
PPGsieLon., (0.204686 0235970451565 0.52875)., 
ea Cieye Fl Galen jes 
= Wwcontrol ("Parent ',a, 
WUmies”, normalized", ... 
Pcabiback 7. closem (ger) ", 
'FontSize',14, 
MEOnUWCLOhNG.. Dold” 7 Us -.. 
"Position', [0.401384 0.0942029 0.207612 0.202899], 
POtringe. OK, 
*Tag’, 'Pushbuttonl1"); 
= Uicontrol('Parent',a, 
Units’; normalized", .«. 
"BackgroundColor', (0.752941 0.752941 0.752941], 
"FonGsize 16, 
‘FontwWeigqhe "7 bold", 
'"POSITCION #102107 260 0.65217470.795ie4s 06217391), 
'String', 'ERROR !', 
‘Style’, 'text', ‘ 
‘Tag’, Staticresxtl"); 
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be=surcontrol ("Parent 7a 
'Units', 'normalized', 
Dee aEecundcolor' Als 752941 Oe 152941 02/5294: |, 
PROntSL Ze >» 1Z, 
'Position', [0.107266 Om333353 0../959048 "0. 371594), 
'String', 'All edit boxes must contain a value to create a new 
Files; : 
PStvyile’ . text’, : 
'Tag', 'StaticTextl1'); 
b = uicontrol('Parent',a, 
‘Units', 'normalized', .. 
Pach aroenacelion (0: 752941 02752947 0.752941], 
"Posttiron ,([U0.03806Z23 0,.050/246 0.930796 0.905797), 
"Style', 'frame', 
‘lad je ELamet "jy 
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APPENDIX X. JANRAD98_PLOTS.M 


This script M-file creates plots for performance airspeed iteration evaluation. It is 


executed in Perf.m. 


$ Script file called from Perf.m to plot airspeed iteration 
$ performance results. 
% JANRAD 98 VERSION 4.0 


EESS OUTPUT - change figure nbrs as needed ESSS 
SESS plotted only if airspeed iteration done %%%% 
SESS add/comment out figures as needed SESS 
if PICK==1] 


$%%% Main Rotor plots - Speed vs T/RHP/TPPangle/Liftpercent - 4 on 1 
page %%%% 

figure (11) 

bitle(* Lift to Drag Ratio for Entire Alreratty 

subplot 221;plot(speed,LoverD); ylabel('W/De'); xlabel('Airspeed 
(kts) ') 

subplot 221; plot(speed,thrust); ylabel('Thrust') 

subplot 222; plot(speed,RHP); ylabel('Rotor Horsepower') 

subplot 223; plot(speed,angle); ylabel('Tip Path Plane Angle (deg)') 
subplot 224; plot(speed,Lperc); ylabel('Wing Lift Percentage') 
ssupplot 223; plot (mu, ctonsig) 


$%$%% Exact numbers for above graphs 
disp(' Speed W/De RHP TPP Angle Lift Percent") 
final=[speed LoverD RHP angle Lperc]; 


$$%% Total Pwr reqd/avail plot %%%3% 

$% first compute pwr avail - 0.91 accts for 9 percent losses fm eng pwr 
output %% 

$% THIS CONTAINS ENGINE SPECIFIC PWR AVAIL NUMBERS %% 


SHPO = 0,.91*{6252; 5436; 4300); sols \ 

SHEO = 0.91*(4037; 3511; 2777); % SL Trop(i08)\ "> use appropriactc ogee 
SHPO = (4226; 3674: 2907]; % 8000 ft ISA / depending on 
conditions 

gamma = 1.4; 

M = speed.* (6080.2/3600) ./sqrt (gamma*1714* (temp+459.6)); % Mach nbr 
delta = (1+(gamma-1)/2*M.%2).%*( (gamma-1)/gamma) ; 

SHPav = SHPO*delta; % Single Eng pwr avail 

aSHPav = 2*SHPav; %$ Dual Eng pwr avail 

Sxmsn = 5040*ones(1,m); %$ Single Eng Xmsn limit 

Dxmsn = 7000*ones(1,m); % Dual Eng Xmsn limit 


Eagure (12) 
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plot(speed, Totpwr);%hold on; 

Ploe(specd, SHPavi(l,:),'--') 

DlOe(specad, csHlavi2,:),speed, dSHPav(3,:),'-.°*) 
plot (speed, Dxmsn) 

title('Total Power Required/Available') 
xlabel('Airspeed (kts)') 

ylabel('Shaft Horsepower' ) 

axis([0 200 0 10000} ) 

axis (‘square’) 


$%$3% Fan-in-Tail Plots - Thrust vs Speed & Pwr vs Speed - on 1 page 
$SSB 

figure (13) 

subplot 211; plot(speed, Tfan,speed,Tfin, '--*) 

title('Anti-Torque Thrust Required and Vertical Fin Thrust Provided') 
yYlabel('Thrust (lbs) ') 

subplot 212; plot(speed,pwrfani, '--',speed,pwrfanp, '-.',speed,pwrfant) 
title('Fan Power Required' ) 

xlabel('Airspeed (kts)') 

ylabel('Power (hp)') 


$%$%% Fan-in-Tail Plot - Speed vs Tot Fan Pwr %%%% 

figure (14) 

plot (speed, pwrfant./RHP*100) 

title('Tail Rotor Power Required.as Percent of Rotor Power') 
xlabel ('Airspeed') 

ylabel('Percent Rotor Power') 


$$%% Plot of Coeff of Pwr vs Coeff of Thrust at Max Airspeed %3%% 
$figure (15) 

$plot (Tcoeff, Pcoeff) 

Stitle('CP versus CT for Maximum Airspeed’) 

$xlabel ('Thrust Coefficient (CT)') 

$ylabel('Engine Power Coefficient (CP)') 

$gtext('SLS Conditions') 


$333 Plot of Coefficient of Thrust vs Figure of Merit %%3%% 
$figure (16) 

S$plot (Tcoeff, figmrt) 

$title('HOGE FM versus CT for SLS') 

$xlabel ('Thrust Coefficient (CT) ') 

$ylabel('Figure of Merit (FM)') 


end 
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% 
% 


APPENDIX Y. STRUCTURE.M 


This script M-file creates the input structure S USER INPUT. 


SEEUCEUre Cons ruction 
JANRAD 98 VERSION 4.0 


SUSE Re UN PUTS SiGe |i... 


‘oye 


UPR Ais ene 
‘temp: , temp, ..... 
‘NATE Sy Vite, 
'GW',GW,... 
"omega',omega,... 
waz. , Naz, 
'thetao',thetao,... 
SOWING o, OW LNG pues; 
EWAN GC, v7 Wllll Gig -sne 
“CLwing ,CLWwing,s.. 
"CDowing 7CDOWING,... 
‘ewing',ewing,... 
ator’, afOilm a... 
Wee ap ono. s 

ey © ae 

72a no eee are 

er pees 6s 
(6 Sal on Mio A aly oa yarar 
Blade 
MECNOLd .t Chord... i 
pte oc by ew 

Btrsie. tests. on 


(r/R) 


BEWISt “twist ,c2%.. 
"woblade',wblade,... 
rnbpe* ,nbe, >. . 
Lai 2 haw 5 cars 
PAL Eli: 4 Give 


Plate Drag Area 


"Afv',Afv,.. 


under disk) 


"Svert',Svert,... 
POVeErt SDVeErt,... 
MeLVere ,CLVEEe os 
"CUOVELrt.  CDOVEEFE, .... 
BONO EEZ Since) arse 
'DHOEEZ “<b horize ss. 
NChLNG Piz * | Ciltotei 2... 6. 


“CDONGHiZ: . ChDOnOnLIZ, .-.. 


Tail 


rEarldisk* ta ticdisk) > 


for JANRADS8 Performance input.m 


oP o Oo oP oP OP oP OP OP cP OH OP OP OP HH CH OP 


oP oO oP 


oP om ® ® CH HH oe oP © © OP OP 


oP 


Pressure Altitude 

Air Temperature 
Airspeed 

Gross Weight 

Rotational Velocity 
Number of Azimuths 
Collective Setting (deg) 
Wing Area 

Wing Span 

Wing CL 

Wing CDo 

Wing Efficiency Factor 
Airfoil Type 

Blade Lift Curve Slope 
Number of Rotor Blades 
Rotor Radius 

Effective Hinge Offset 
Length of Inner Non-Aerodynamic Portion 


Root Chord 
Rotor Blade Taper Ratio 
Rotor Blade Taper Ratio Starting Point 


Blade Twist 

Blade Weight 

Number of Blade Elements 

Aux Thrust 

Fuselage Equivalent Horizontal Flat 
Vertical Projected Area (fuselage area 
Area, Vertical Tail 

Span, Vertical Tail 

Lift Coefficient, Vertical Tail 

Profile Drag Coefficient, Vertical Tail 
Area, Horizontal Tail 

Span, Horizontal Tail 

Lift Coefficient, Horizontal Tail 
Profile Drag Coefficient, Horizontal 


Tail Under Main Rotor (1-Yes, 2-No) 


bea: 








APPENDIX Z. STRUCTUREI.M 


This script M-file creates the output structure S PERF OUTPUT. 


$ Strueture Consruction for JANRADSS Performance output sm 
%$ JANRAD 98 VERSION 4.0 


SEPP RE OULU LSS crc c (a). 


“DLuse 7 DLUSe , 2% 
DH EOtO rier OtOr. <5. 
Viwind Sf uwWingd, .. . 
"Dwing',Dwing,... 
yA) oleh arly ae Fete) qn kA 
iD oven au kya) os) qa ly. 
“IVE saver 2 2. 
 DVErE Dvert..... 


ralphal *,alpohal* 57.3.3... 
HbeGac ,DeCaA0* a7. 3, ..« 


Miah Zar 2 exe: 


MEnetao- .chetao* 5/23, ... 
btehetale thetalc* 57.37. - 
‘thetals’”thetals*57.3,.4.. 
isOliaiey | sOlLidity,;.a.. 


Pe Dlbe sli... 
J7 GS 3 eae 
Bel Sig ', Clasctd, =. « 
‘CQO sig',CQ sig,... 
peHesig¢ 7 eiHics iG nn. 


(Machteip* ;Machtaip, ....< 


ST TN er 

3) Raed tere eae 

Le EGQEOY{.P GOtO E,... s-. 
SORGCOE ,OrOCOL) >? 


$Fuselage drag 

$Rotor drag 

$Wing lift 

$Wing drag 

$Horizontal tail 11ifT 
$Horizontal tail drag 
$Vertical tail side force 
$Vertical tail drag 

$Tip path angle 

Rotor coning angle 

$Location of mean thrust 
S$Rotor Collective pitch at .7 r/R 
$list lat cyclic term-Al (deg) 
$ist long cyclic term-Bl (deg) 
$solidity (sigma) 

$Disk loading 

$Figure of Merit 

$CT/sigma 

$CQ/sigma 

$CH/sigma 

$Tip mach of the adv. blade 
tAdvance ratio 

SRotor thrust required (TPP) 
$Rotor power required 

Rotor torque 
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APPENDIX AA. STRUCTURE2.M 


This script M-file creates the vector structure S MATR_VEC. 


%$ Structure Consruction for JANRAD98 Perf.m 
% JANRAD 98 VERSION 4.0 


S_MATR VEC=s CEUGCE (2. 


a ap a 
(DSay 7 PsilGu... 
Nalin NJ eee 


Steheta atnetag. 
"beta betat,.... 
"alpha',alpha,... 


“TDSi pi psi7 ss 
'Mosi.' 7 Mpisdiy sss 


/DMpsi 7 DMpsi1, . .- 


viel AloUlAia Ss 
elu! (chy ae 
“1s 1D)" (ol BF 


'cblade',cblade,... 


(Gityy AAs 
"CD. CD); 


$radial distance 

$azimuth 

S$induced velocity 

$collective pitch 

$coning angle 

Sangle of attack 

$Thrust at azimuth position 
$Thrust Moment at azimuth position 
Drag Moment at azimuth position 
S$incremental Thrust 

S$incremental Thrust moment 
$incremental drag 

tblade chord 

$incremental CL values 
S$incremental CD values 


WAL 








APPENDIX AB. UNSTRUCTURE.M 


This script M-file decomposes the S PERF INPUT structure into 36 separate 
input parameter variables. These variables are corrected to proper units used in Perf.m 


computations. 


$ Structure De-Consruction for JANRAD98 Perf.m 
% JANRAD 98 VERSION 4.0 


PA] SPR the dre, 
temp=S_ PERF INPUT.temp; 

Vin -Seeh RE veer. Vink: 

GW=S_ PERF INPUT.GW; 

omega=S _PERF_ INPUT. omega; 
naz=S PERF_ INPUT. naz; 
thetao=S _PERF_ INPUT. thetao; 
Swing=S _ PERF_ INPUT. Swing; 
bwing= S PERF INPUT. bwing; 
CLwing=S_PERF INPUT.CLwing; 
CDowing=S_PERF_ INPUT.CDowing; 
evinga osc hR be INPUT. ewing; 
afoil=S _PERF_ INPUT. afoil; 

a= S_PERF_INPUT. a; 

b=S PERE SINEUT.b; 

R=S _PERF_ INPUT .R; 

e=S _PERF_ INPUT .e; 

grip=s _PERF_ PN PU iG ip 
rchord=s _PERF_ ENEUT. rchord; 
tr=S _PERF_ INPUT. er, 

trst=S _PERF_ TNEUTErSt; 
twist= ae PERF_ INPUT. twist; 
wblade=S _PERF_ INPUT. wblade; 
nbe=S _PERF_ INPUT. nbe; 

Taux=S _PERF_ INPUT. Taux; 

Afh=S PERF_ INPUT. Afh; 

Afv=S_ PERF INPUT. Afv; 

SVveTE=S 7PERECENEUT overt; 
bvert—S PERRCINEUL  bVert, 
Civert=S_PERF INPUT.CLvert; 
CDovert=S_PERF_INPUT.CDovert; 
Shoriz=S_PERF_INPUT.Shoriz; 
bhoriz=S PERF _INPUT.bhoriz; 
CLhoriz=S_ PERF_INPUT.CLhoriz; 
CDohoriz=S PERF INPUT.CDohoriz; 
taildisk=S PERE INPUT. calldrsk; 
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APPENDIX AC. UNSTRUCTUREI.M 


This script M-file decomposes the S USER_INPUT structure into 36 separate 
input parameter variables. These variables are the actual values the user types or loads 


from a previously saved file 


%$ Structure De-Consruction for JANRAD98 Perf.m 
% JANRAD 98 VERSION 4.0 


PA=S_USER_INPUT. PA; 
temp=S_ USER _INPUT.temp; 

Vint =se USER INEUT Vint; 

GW=S_ USER INPUT.GW; 
omega=S_USER_ INPUT.omega; 
naZ=S_ USER _ INPUT. naZ; 
thetao=Ss _USER_ INPUT. thetao; 
Swing= S_USER_ INPUT. Swing; 
bwing=S_USER_INPUT.bwing; 
CLwing=S_USER_INPUT.CLwing; 
CDowing=S_USER_ INPUT.CDowing; 
ewing=S_USER_ INPUT.ewing; 
atorl—S USERD INPUT. atoit; 
a=S_USER INPUT.a; 

b=S_ USER INPUT.b; 

R=S_USER INPUT.R; 

o— > Ustke@ENPUT =e, 

grip=S_ USER INPUT.grip; 
rchord=S _USER_ TNEUL] tocnord- 
tr=s _USER_ INPUT. Ce 

trst= S _USER_ INPUT. Crse 
twist=S _USER_ INPUT. twist; 
wblade=s _USER_ INPUT.wblade; 
nbe=S _USER_ INPUT. nbe; 

Taux=S _USER_ INPUT .Taux; 

Afh=S _USER_ INPUT. Afh; 

Afv=S USER INPUT. Afv; 
Svert=S USER _INPUT.Svert; 
bvert=S USER INPUT.bvert; 
CLvert=S USER INPUT.CLvert; 
CDovert=S_ USER INPUT.CDovert; 
shoriz=S USER INPUT.Shoriz; 
bhoriz=S_USER INPUT.bhoriz; 
CLhoriz=S_USER_INPUT.CLhoriz; 
CDohoriz=S _USER_ ENPUT. ChohnoLiz.- 
taildisk=S _USER_ INPUT. taildisk; 
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APPENDIX AD. UNSTRUCTURE2.M 


This script M-file decomposes the S PERF OUTPUT structure into 25 separate 
Output parameters. These variables are displayed in the performance output figure 


window. 


* structure De-Consruction for JANRADOS performance output.m 
% JANRAD 98 VERSION 4.0 


Dilsc=oagLenre OUTPUT. DEUse; 
rotors EERE IOULEPUL: HrOtor, 
Lwing=S_PERF_ OUTPUT. Lwing; 
Dwing=S_PERF OUTPUT. Dwing; 
Lhoriz=S _PERF_ OUTPUT. EnOri Z >; 
Dhoriz=S _PERF_ CULTEUT - LHOEYZ; 
Lvert= S PERF OUTPUT. LYert:; 
DVvere—o PERE LOUTPUT- Dvert; 
alphaT=S_ PERF _ OUTPUT.alphaT; 
betao=S_ PERF OUTPUT. betao; 
rT2=S _PERF_ OUTPUT. rT2; 
thetao=s _PERF SOUTPUT. tnetao; 
thetalc=S _PERF_ OUTPUT. thetalc; 
thetals= =S PERF OUTPUT. thetals; 
solidity—-s PERF OUTPUT. solidity; 
DL=S PERF OUTPUT. DL; 

FM=S_ PERF OUTPUT. F&M; 

Siect Oo se Ee REOUlPUT Clisig; 
CO sig=s_PERF OUTPUT.CO sig; 
CH sig=S_PERF OUTPUT.CH_ sig; 
Machtip=S _PERF | OUTPUT.Machtip; 
mu=sS _PERF OUTPUT. mu ; 

T= =S_ PERF OUTPUT. de 
ErOcCot—> PERE ZOUTPUT . PEGrOE; 
Qrotor=S PERF OUTPUT.Qrotor; 
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APPENDIX AE. UNSTRUCTURE3.M 


This script M-file decomposes the S MATR_VEC structure for use in the print 


and save commands. 


$ Structure De-Consruction for JANRAD98 Perf.m 
&  JANRAD 98 VERSION 4.0 


ES > MATES. Gy 

psi= 3S MATR VEC.psi; 

= > MATR EC. vir; 
Eeta= 5 MAER VEC. theta, 
betat= S MATR VEC.betat; 
alpha= S MATR VEC.alpha; 
hos Fork evel psi; 
Mpsi= S_ MATR VEC.Mpsi; 
DMpsi= S_MATR VEC.DMpsi; 
dT= S MATR VEC.dT; 

dM= S MATR VEC.dM; 

aD= S MATR VEC.dD; 
eblade= 5 MATR VEC. cblade; 
CL= S MATR VEC.CL; 

CD= S MATR VEC.CD; 
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APPENDIX AF. PERF.M 


This script M-file is the primary computational routine for JANRAD 98. It is 
launched from either the iteration method fcn.m or iteration_parameter_fcn.m 


Switchyard Callback function. 


% Perf.m 
$ Main Performance computaion routine. 
% JANRAD 98 VERSION 4.0 


global S PERF INPUT S MATR_VEC 


feo eek ENEU Legmip. — le= 10) 
SePEeRE@eNPUT getp—le-l0; 
end 


oe > PERE @ONPUT- Swing <92e=167 
S_PERF_INPUT.Swing=le-10; 
end 


Peo PORE INR oowingea. Le—1 0, 
S_ PERF _INEUT.bwing=le-10; 
end 


Peso ERE INPUT sewing = fea lG), 
Se EERE eINEUT-ewing=le-0; 
end 


if S_ PERF_INPUT.Shoriz < le-10, 
o PERF INPUT.horiz=le-10; 
end 


Peo ERE LNPUL bihoriz, < le-i0, 
© PERF INPUT. bhoriz=le-10; 
end 


sie Sc eRe INRUG overt were 10, 
3 PERF INPUT.Svert=le-10; 
end 


dbf S PERE INPUT. bvert <sle=10;, 
S_PERF_INPUT.bvert=le-10; 
end 


S_PERF_ INPUT.Vinf=S PERF INPUT. Vinf*1.68894444; 


St PERF _ INPUT. Cw¥et= abs (S_ PERE INEUMaewist)/o7 <3; 
5 _PERF_ INPUT. thetao=S _PERF_ INPUT. thetao/57.3; 


ey 


unstructure 


Pid -LOpen\ print temp , wt); 


fprintf(fid, '\t 
ferimer (ir cyeey tc 


kts\n',Vinf/1. 69); 


foraner(fid, Vt 
ENE Gem) > 
forint. (bid, vt 
oncaghere (e.s)) Ye 
Pprinet (tid, \t 
foeimet (fia, Ve 
fprrner (tid, *\t 
fEyM 7 benord):-> 
former (era, © Gt 
PeEWise™ 57 o5)> 
Porminee ( fads! Vt 
Eorinin fla,  \t 
lbs\n',wblade) ; 
le ahenlio eye (atilole se 


rads/sec\n',omega) ; 


forint t (fied Yt 
fpeineh (cid Vt 
forint’ (Lid, Vt 
Geac Nn. Ath) 

forint e (ila; \t 
Pee ne ALy); 

Perimer (fra, Vt 
hEoZ Nn , OWING) | 
PpLriner (fia, \ Cc 
ft\n',bwing) 
foriner (fra, Ve 
Poa Miia £2 cya yt 
\n',CDowing) ; 

EOrInt £( tid,” yt 
POEinte (fid,: Vt 


se 


fee2 in, shoriz) 


FOEMMeEr (ria, \t 
lee \ ale fle) aleviosbya)| 
forint e (fid,: \c 
Ma ebhor 12) 
EpPrEantr (fia. Ve 
Mie DOnOr1Z )- 
Eoriner(fid, *\t 
EE 2 Nie overt): 
POE Ine rrr, \e 
ft\n',bvert); 
forIner (rid, * \t 
ED rintttcLra, \t 
\n',CDovert); 
EpGante( ria,” Vt 
ies \n', Taux) ; 
fporaner (rid, “\t 


MOD Leiner ne ENEUT 


Forward velocity 


Semperacure 


Pressure altitude = 


Gross weight 
Number of blades 
Rotor radius 


Blade root chord = 


Blade twist 


Blade lift curve slope 
Blade weight 


Rotational velocity 


Blade grip length = 


Hinge offset 


Equivalent flat plate area = 


Vertical projected area 
Wing area 


Wing span 


Wing CL = 
Wing CDo = 


Wingsetttcleneyerde tors 
Horizontal tail area = 


Horizontal tail span 
Hora zontal taal ch 
Horizontal tarmac De 
Vertical tail area 


Vertical tail span 


Vertical tail CL = 
Vertical tail €po = 


Auxiliary thrust 


Number of Azimuths 
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$6. 
$6. 
$6. 
$6. 
$6. 
SO. 
$6. 
iA S 


$6. 


$6 


$6. 
$6. 
$6. 
ao. 
36. 
$6. 
36. 


$6. 
$6. 


$6. 
8) 


$6. 
$6. 
$6. 
$6. 
$6. 


$6. 
oo. 


$6. 


$6. 


panes Dl. Tinka) y 


Of 
Or 
Of 
Of 
Of 
Ze 
Ze 
Z£ 


Ze 


2f 


fae 
Zee 
Ze 
/a3€ 
2G 
Ze 
Zo 


Ze 
4f 


Pani 
ane 


Ze 


Pane 


4f 


oe 


Zt 


Zi 
4f 


Of 


Of 


degs 


£E\N gee: 
lbs\n',GW); 
\n',b); 
ft\n',R); 


degs\n',- 


\n*,a); 


EU\n Ue gri ps 


ft\n',e); 


\n*; CLwing): 


\n',ewing) ; 


\n', CLVert); 


Xi haz) ; 


Pp Pine fb Collective Pitch =s%6.25 
degs\n’*,thetac*5/™3s 77; 


prints (hada Airfoil @ype = %6.Ofeime, afcil) > 
Bovenhehaee (resisl, © ste Taper Ratio = %6.2f \n',tr); 
isopaabayans ((ag2hole Nie Taper Ratio Starts At = $6.2fieft\n',trst); 
MG ewerhonare (eirel, Mc Number of Blade Elements = %$6.0f \n',nbe); 
fprint£(fid, *"\t Tail Under Main Rotor (l-yes 2-no) = %6.08 


Nee tad dad sk): 
fel ose (fid)-> 


switch PICK 


case l 
if isempty (MINUM) 
MINUM=0; $ sets min airspeed default to 0 kts 
end 
if isempty (MAXUM) 
MAXUM=160; %$ sets max airspeed default to 160 kts 
end 
if isempty (INTER) 
INTER=20; % sets INTERval default to 20 kts 
end 


if MINUM > MAXUM % ensures INTERval is the correct sign 
INTER=-abs (INTER) ; 

else 
INTER=abs (INTER) ; 

end 

doitl='airspd=itervar;'; 

doit2='Vinf=airspd*6080.2/3600;'; % converts kts to ft/s 

m= (MAXUM-MINUM) /INTER+1; % computes reqd nbr of elements for 

VeECEOLrs 


speed=zeros (1,m); % vector for each airspeed values 
mu=zeros (l,m); %$ vector for adv ratio values 
Lperc=zeros (1,m) ; %$ vector for wing lift percent values 
ctonsig=zeros(1l,m); % vector for CT/sigma values 
LoverD=zeros(1,m); % vector for tot acft W/De 


theone='airspd'; 


case 2 
if isempty (MINUM) 
MINUM=0; % sets min altitude default to O ft 
end 


if isempty (MAXUM) 
MAXUM=15000; % sets max altitude default to 15000 ft 
end 


if isempty(INTER), 


INTER=1000; %$ sets INTERval default to 1000 ft 

end 

if MINUM > MAXUM % ensures INTERval is the correct sign 
INTER=-abs (INTER) ; 

else 
INTER=abs (INTER) ; 

end 
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Gdolel=— PA=i1tervar;'; 

%$ sets correct ISA temp 

% for given alt (deg F) 

%$ note: S9=SLS, 103=trop 

doit2='temp = 59-1.981lle-3/.5555*PA;'; 

m= (MAXUM-MINUM)/INTERt+1; % computes regqd nbr of elements for 
VEGEOLrs 

if REGIME== 

Ptige=zeros (l,m); 

end 

altitude=zeros (l,m); 

theone='PA'; 


case 3 
if isempty (MINUM) 
MINUM=10000; %$ sets min GW default to 10000 lbs 
end 
if isempty (MAXUM) 
MAXUM=20000; %$ sets max airspeed default to 20000 lbs 
end 


we isempty (INTER), 
INTER=2000; % sets INTERval default to 2000 lbs 
end 
if MINUM > MAXUM % ensures INTERval is the correct sign 
INTER=~-abs (INTER) ; 
else 
INTER=abs (INTER) ; 
end 
doitl='GW=itervar;'; 
doitZ=—' i; 
m= (MAXUM-MINUM) /INTER+1; % computes reqd nbr of elements for 
VECCors 


wt=zeros (1,m); %$ vector for GW values 

Tcoeff=zeros (1,m); %$ vector for Coeff of Thrust values 
Pcoeff=zeros (l,m); %$ vector for Coeff of Power values 

figmrt=zeros(1,m); %$ vector for Figure of Merit values 


theone='GW'; 


case 4 
if isempty (MINUM) 
MINUM=0; $ sets min blade twist default to 0 deg 
end 
if isempty (MAXUM) 
MAXUM=-12; %$ sets max blade twist default to -12 deg 
end 
if isempty(INTER), 
INTER=-2; % sets INTERval default to -2 deg 
end 


if MINUM > MAXUM % ensures INTERval is the correct sign 
INTER=-abs (INTER) ; 

else 
INTER=abs (INTER) ; 

end 
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doitl='TWIST=itervar; ; 

doit2='twist=TWIST/57.3;'; % converts degrees to radians 

m=abs (MAXUM~MINUM) /INTER+1; % computes reqd nbr of elements for 
vectors 

thetat=zeros(1,m) ; %$ vector for blade twist values 

theone='TWIST'; 


case § 
if isempty (MINUM) 
MINUM=1; %$ sets min blade taper ratio default to l 
end 
if isempty (MAXUM) 
MAXUM=0.5; $ sets max blade taper ratio default to 0.5 
end 
if isempty (INTER), 
INTER=-0.1; % sets INTERval default to -0.1 
end 
if MINUM > MAXUM %$ ensures INTERval is the correct sign 
INTER=-abs (INTER) ; 
else 
INTER=abs (INTER) ; 
end 
doitl='tr=itervar;'; 
doit2=' '; 


m=abs (MAXUM-MINUM) /INTER+1; % computes reqd nbr of elements for 
vectors 

taper=zeros(1l,m); % vector for blade taper ratio values 

theone='tr'; 


case 6 
if isempty (MINUM) 
MINUM=0.1; %$ sets min blade twist default to 0.1 r/R 
end 
1f isempty (MAXUM) 
MAXUM=0.9; $ sets max blade twist default to 0.9 r/R 
end 
if isempty(INTER), 
INTER=0.1; %$ sets INTERval default to 0.1 r/R 
end 


if MINUM > MAXUM % ensures INTERval is the correct sign 
INTER=-abs (INTER) ; 


else 

INTER=abs (INTER) ; 
end 
doitl='trst=itervar;'; 
aoit2=" *; 


m=abs (MAXUM-MINUM) /INTER+1; % computes reqd nbr of elements for 
VecCcoOrs 


start=zeros(1,m); % vector for taper ratio start position 
values 
theone='trst'; 


case 7 
if isempty (MINUM) 
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MINUM=50; $ sets min wing area default to 50 sq ft 


end 
if isempty (MAXUM) 
MAXUM=8 0; %$ sets max wing area default to 80 sq ft 
end 
if isempty(INTER), 
INTER=5; % sets INTERval default to 5 sq ft 
end 
if MINUM > MAXUM % ensures INTERval is the correct sign 
INTER=-abs (INTER) ; 
else 
INTER=abs (INTER) ; 
end 
if isempty (AR) 
AR=6; % sets aspect ratio default to 6 
end 


doitl='Swing=itervar;'; 
doit2='bwing=sqrt (AR*Swing);'; % computes wing span given AR and 
Wing area 
=abs (MAXUM-MINUM) /INTER+1; % computes reqd nbr of elements for 
vectors 


area=zeros(1,m); % vector for wing area values 
theone='Swing'; 


case 0 
MINUM=1; %\ 
MAXUM=1 ; (a UsaCetnese Values Ehe for loop is cffectively 


INTER=1; %/ non-existent - will go through only once 
doitl=' '; 
doiltZ—") "> 


end 


$%% common vectors for compilation of iterative data %%% 


if PICK~=0 
thrust=zeros(1,m); 
RHP=zeros(1,m); 
Totpwr=zeros (1,m); 
angle=zeros(1l,m); 
Drag=zeros(1,m); 
end 


$%% account for vertical drag on wing and horizontal tail %%% 
Afv1l=Afv; $ puts the vert profile area into temp variable Afvl 
if taildisk== 


Afv=Afvl1t+3* (SwingtShoriz) ; % this assumes a vertical Cd of 1.2 
for the 


elseif taildisk== $ wing and horiz tail and a vertical Cd of 
Afv=Afvl1l+3*Swing; % 0.4 for the fuselage (i.e. 1.2/.4 = 3) 
end $ thus making wing/tail effectively larger 


% when hvr thrust calc using Cd=0.4 
$%% NOTE: All total power calculations are based on a Fan-in-Tail %%% 
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SSS anti-torque device. These equations must be modified %%% 
SSS if using a conventional tail rotor or NOTAR system 33% 


$%% Fan-in-Tail Parameters %%% 


Afan’="pi* 3.522; %$ Area of fan; radius=3.5 ft 
lfan = 41; % Distance fm M/R hub to fan hub 
sigmafan = 0.143; % Solidity of fan 
cdfanm—m0 7015; % Cdo of fan blade 


$%% vectors for Fan-in-Tail calculations %%% 


if PICK~=0 
Tfin=zeros(l,m); 
Tfan=zeros (l,m); 
pwrfani=zeros(1,m); 
pwrfanp=zeros(1,m); 
pwrfant=zeros(1,m); 
end 


$%% Beginning of Iteration Loop %%% 
p=0; 
for itervar = MINUM: INTER: MAXUM 
eval (doitl1) % reads and evaluates the string 'doitl' assigned 
above 
% assigns itervar to the reqd var needed to complete iteration 
eval (doit2) % reads and evaluates the string 'doit2' assigned 
above 
% used if input variable needs some type of manipulation 
p=pt+1; % used to move to next vector element 
set(H STATUS, 'String', 'ROTOR PERFORMANCE ROUTINE’ ) 
cle 
Set(HvolAlUsSl, -seuing.,. ctaresh lapsed aame.} 
pause (3) 
trim 


SSS *** Calculation of output parameters *** SSS 


load perftemp % Eccles addition - program was not recognizing 
% mchord and DMpsi. 

Protor=mean (DMpsi) *b*omega/550; % Rotor pwr reqd, NOT total pwr 
Qrotor=mean (DMpsi) *b; % Rotor torgue reqd 
solidity=b*mchord/ (pi*R) ; 
CQ=Qrotor/ (Adisk*rho*Vtip*2*R) ; 
CH=Hrotor/ (Adisk*rho*Vtip%2) ; 
CT sig=CT/solidaty; 
CQ sig=CQO/solidity; 
CH sig=CH/solidity; 
Machtip=(Vtip*cos (alphaT)+Vinf) /(49.05*sqrt (temp+460) ); 
Li Vint 16..9- 

DL=T/ (pi* RZ); 

FM=(T*sqrt (DL/ (2* rho) ))/(S50*Protor) ; 
else 

DL=0; 
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EFM=0 ; 
end 


$%%% Compute Fan and Access/Xmsn Power Reqd in Fwd Flt %%%% 

vifan = abs(sqrt((Qrotor/lfan-Lvert) /(rho*Afan) )); %$ use abs in 
case complex 

(Qrotor/lfan-Lvert<0) ; 


Thrustf = Qrotor/lfan; $ Thrust reqd for anti-torque 

Pfani = 0.5*rho*Afan*vifan%3/550; $ Fan induced pwr 

Pfanp = rho*Afan*Vtip*3*sigmafan*cdf£an/4400; =» Pan pororele pwr 
Pace — 106. 25+0.012/5*Frotor; % Pwr for accessories 

Ptot = Protor+PfanitPfanptPacc; %$ Total pwr reqd 

CP=Ptot*550/ (Adisk* rho*Vtip%’3) ; %$ Coefficient of Pwr based on 


PEOE 


$%$%% HIGE total pwr %%%% 


if REGIME== 
viave=sqrt (DL/ (2*rho) ); % Induced vel - hoge 
delpwr=T* viave*0.23/550; $ Decr in pwr due to grd effect 
Prtrige=Protor-—delpwr; S$ Rotor pwr hige 
Tige=Prtrige*550/ (omega*lfan) ; % Rotor thrust - hige 
vifige=sqrt (Tige/(rho*Afan) ); % Induce vel - hige 
Pfanige=0.5*rho*Afan*vifige*3/550; $ Fan pwr ~- hige 
Paccige=106.25+0.01275*Prtrige; % Access pwr - hige 
Ptotige=Prtrige+PfanigetPfanp+Paccige; % Tot pwr - hige 

end 

WonDe=GW/ (550*Ptot/Vinf) ; $ Total acft W/De (lift/drag) 


$%% Collecting and vectoring all the calculated data %%% 
%$%% vectors specific to desired iteration 


switch PICK 

case l 
speed(p)=airspd; % vector of airspeeds 
mu(p)=Vinf/Vtip; % vector of advance ratios 
Lperc(p)=Lwing/GW*100;% vector wing lift percentage 
Etons1qige)—-Cl sig; % vector of CT on sigma 
LoverD(p)=WonDe; % vector of W/De 


Case =z 
altitude(p)=PA; % vector of altitudes 


Casene 
wt (p) =GW; % vector of GW 
Teoeff(p)=CT; PEVCGCCOL "OF COCTrLGloncror slice 
Peoeff (p)=CP; % vector of Coefficient of Power 
figmrt (p)=FM; % vector of Figure of Merit 


case 4 
thetat (p)=TWIST; % vector of blade twist 


aise 2 
taper (p)=tr; $ vector of taper ratio 
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case 6 


start (p)=trst; %$ vector of starting position of blade taper 
case 7 
area (p)=Swing; % vector of wing wetted area 
end 
if REGIME== 
Ptige(p)=Ptotige; % vector of total hp in hige 
end 


$%% vectors common to all iterations 


if PICK~=0 
thrust (p)=T; % VECEOE OL Main erotor Ehrust 
RHP (p)=Protor; %$ vector of rotor hp reqd 
Totpwr(p)=Ptot; % vector of total hp reqd 
angle (p)=alphaT*57.3; % vector of tip path plane angle 
Drag (p)=Drotor+Dfuse+Dvert+DhoriztDwing; % vector of tot acft drag 
Tfin(p)=Lvert; eS veceor Of vertical fin Litt 
Tfan(p)=Qrotor/lfan; $ vector of fan thrust 
pwrfani(p)=Pfani; % vector of fan induced pwr 
pwrfanp(p)=Pfanp; % vector of fan profile pwr 
pwrfant (p)=Pfanit+Pfanp; %$ vector of fan total pwr 
eval('theone'); % displays current value of iterative variable 
end %$ ‘'end' needed to complete the 'for' loop 


end $ this is the ‘end’ needed to complete the 'for' loop 


$%$$%$ Save iteration data for future processing %%%% 

$%$% format: save <filename> varl var2 var3 %%% 

$%% (mote: all variables must be valid or will get error) %%% 
$%$% Works well to create short m-file to graph this calculated data %%% 
SSB just use ‘load <filename>' at the beginning of the file to %%% 
SSS read all the vectors which are stored in the .mat file %%% 


if PICK~=0 

Save output Totpwr angle RHP thrust Drag 

save tailop Tfin Tfan pwrfani pwrfanp pwrfant 
end 


switch PICK 
case l 


Save extral speed mu ctonsig LoverD Lperc RHP 


case 2 
Save extra2 altitude 


case 3 
save extra3 wt Tcoeff Pcoeff figmrt 


case 4 
Save extra4 thetat 
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Gase 5 
Save extra5 


case 6 
save extra6é 


case 7 
save extra7 
end 


if REGIME== 
Save extra9 


end 


if PICK== 


taper 


start 


area 


Ptige 


janrad98 plots 


end 


Set(H STATUS, (String 77,COMPLETING CALCULATIONS 
Sset(H slAnus), String jae hun slapsed Time as © 


seconds']) 
pause (3) 


% *** output to diskee(tex tefl le) 


kkx* 


fid=fopen(*print templ", 'wt"); 


POEInee fad, 
Por net (fac, 
fprintt(fid, 
EDEIMeL (ela, 
forint. | fac, 
fprint& (iid, 
foeiner (iid, 
EDEING HE Id, 
fpEIneiHi1a, Vt 
TF CWS ous); 

fprint Méfid,"\t 
moves polieoe (mex tely ! ic 
TDEiner (Lid, \c 
rads/sec\n', 
Poriner (fia, 
fPomints ( fic, 
fporiner (id, 
Loewe Lid, 
EPrIner (hia, 
EPEINeL (fia, 
Eprinte(tid, 
fFoOLrIneE ELA, 
Fournier (fid, 
Eprineiltid, 
fpHLiner (rid, 
fPprintete (fad, 
Eprinwe(fid, 


"\t 
Vt 
Ne 
Se 
EN te 
"\t 
"\t 
ove 
ete 
"\t 
NE 
"\t 
A 


omega) ; 


Equ 


Bee) clo) Oy Br laces 


Forward velocity 
Temperature 
Pressure altitude 
Gross weight 
Number of blades 
Rotor radius 
Blade mean chord 
Blade twist 


Blade lift curve slope = 


Blade weight 
Rotational velocity 


Blade grip length = 


Hinge offset 
ivalent flat plate area 
Vertical projected area 

Wing area 

Wing span 

Wing CL 

Wing CDo 

Wing efficiency factor 
Horizontal tail area 
Horizontal tail span 
Horizontal tail CL 
HOrizontal™ tail CDo 
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") 


Num2eser (fix (toc). 


ae Nigel 


$6. 
$6. 
oO. 
$6. 
$6. 
36. 
evr 
SO. 


36. 
SO. 
$6. 


$6. 
56s 
oC: 
60. 
(8) 
$6. 
$6. 
oC. 
hoe 
Soe 
50. 
$6. 
oC: 


Of 


kts\n'; Vint/1. 6097" 


degs F\n',temp); 
ft\n', PA); 
lbs\n', GW); 
\n',b); 

fev eee 
feE\n ' ~mchord); 
degs\n',- 


\n',a); 
lbs\n',wblade); 


£E\N 7 OElD) 
BEV pele 

Bt 2X Ne, Poe 
f£t*2n ' PAtV) > 
£E°2Z2N\n Swing je 
fE\n = bwing) : 
\n , Ciwing)e 
Vn, CLDOwing):, 
\n',ewing) ; 
feo2 ne Sno rss2)) 7, 
fE\M | DnOrI Zz), 
Xn 2ennoriz)s, 
Nae! | CDeme@riz) ; 


EprimtG (hie \t Vertical tail area = %6.2f ft°2\n 7 Svert); 


Porn: tira NS Vertical tail span = %6.2f ft\n',bvert); 
FOErI Ne tric Ne Vertical tail (Cls= 46228 Vn civerte 
EDrint feta ee Vertical tail CDo = %6.4f \n',CDovert); 
EPrints (eid ave Fuselage drag = %6.0f lbs\n',Dfuse); 
Eprintth (sada Rotor drag = %6.0f lbs\n',Hrotor); 
Ep rint th Lice Wing lift = $6.0£ 1Jbs\n", Lwing),- 
ford ietiGrd cl wN Cc Wing drag = %6.0f lbs\n',Dwing) ; 
i@epealgiane | jeakels ve Horizontal tail litt = 36.0f lJbs\n' ihoriz)- 
forint (fade Ve HorizZontalptail drag = -t6.0f lbs\n'7Rnorrz):, 
jd oPank panel allo Bau ce Vertical tail side force = %6.0f lbs\n',Lvert); 
EDLINGE Ela. \E Vertical tail drag = %6.0f lbs\n',Dvert); 
BoLiner (faa, \t Auxiliary thrust = %6.0f lbs\n',Taux) ; 
iptaint £ (Lid, + \et Tip path angle = %6.2f 

degs\n', alphaT*57.3) ; 

Pere r (hid, Vt Rotor coning angle = %6.2f 


deqs\n'*;betac*57.3); 

fprinth(fid, “\elocation of mean thrust (r/R)@="8360:.2f \n",xrTZ); 

pred (ec. \ t Collective pitch at .7 r/R = %6.2f 
degs\n*,;thetac*o7.3); 

forint i(rid,; *\t Ist lat cyclic term-Al.{deg) = $6.2f \n",thetale~s7<3), 
Perinti(fid, *\tist long cyclic. term-Bl (deg) = %6.2f \m",thetals*57-4), 


boment tL (fia, \t Solidity = @6.3f@ \n*,selidiry) 
BOrtntt (fid, * \t Disk loading = $6.2f Ibs/fle2\n PL) 
Printer (Lid, “\t Figure of Merit = %6.2£ \n',FM); 

Berintet (tid, Vt CT/sigma = %6.3f \n',CT sig); 
Bpreantf (fid, '\t CQ/sigma = %6.4f \n',CQ sig); 
mormemcr( fra.) Vt CH/sigma = %6.4f \n',CH_ sig); 
EoeIntei (iid, *\t Tip mach of the adv. blade = %6.3f \n',Machtip) ; 
Sorintitfid, *\t Advance ratio = %6.3f \n',mu); 
moeinti{fid, "\t Rotor thrust required (TPP) = %6.0£ lbs\n",T);; 
morinti( fid, *\t Rotor power required = %6.0f h.p.\n',Protor); 
Borantt (Lid, '\t ROLOE torques — sG.,0E te _Ibs\f  ,Orotrer)- 


me Lose (fid) ; 
% *** Configuring variables for output *** 


theta=theta*57.3; 

betat=[betat twist* (0.7-(Reff+(R-Reff) /2)/R)]*57.3; 
alpha=alpham*57.3;,alpha=[alpha zeros(size(psi))]; 
Mpsi=Mpsi(:,length(Mpsi(1,:))-1); 

adM=[dM ddM]; 

psi—psi*5/.3; 

r=(r (R-(R-Reff)/2)]; 

vi=[vi 0]; 


Sec SAIS, String, SfAND BY®POR CUITEUT...*) 
pause (3) 


structurel 

structure2 

Seu (HeSsTATUS; string ,°~) 
Perlormance JOuEpuUE 
closet hel ToMe in 
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APPENDIX AG. TRIM.M 


This script M-file 1s a subroutine to trim the rotor system. It is called in Perf.m. 


= Trimm 


$ Trim routine for collective/cyclic. 
$ JANRAD 98 VERSION 4.0 


set(H STATUS, 'String', 'EXECUTING ROTOR TRIM ROUTINE") 

set(H STATUS1, 'String',['Run Elapsed Time is ' numZstr(fix({toc)) | 
seconds']) 

pause (3) 


% *** calculation of required parameters *** 
mio—. OUZS77" (=. VO0OS1* PAt(—.002"tempti.. lo) ); 
% 7 * te TIeSt Guess at rotor pProfllevsdrage( nh toOrece m= -.: 


vr Vinkteaw 1659 
Drotor=0; 
else 
DrOtor—Vint* (rho/. 002377) ; 
end 


e=U. oO rnos Vint Z; 

Adisk=pi*R%*2; 

Vtip=omega*R; 

Dfuse=q*Afh; 
CDwing=CDowing+ (CLwing%2/ (ewing*pi* (bwing*2/Swing) )); 
SPnomd z—CDohori z+ (CLhNOriz 2/ {. 84pi*({bhoriz-27Shoriz))})> 
CDvert=CDovert+(CLvert*2/(.8*pi* (bvert*2/Svert) )); 
Dwing=q* CDwing*Swing; 

Dhoriz=q* CDhoriz*Shoriz; 

Dvert=q*CDvert*Svert; 
Dftotal=(Dfuse+Dwing+Dhoriz+Dvert) -Taux; 
Lwing=q*CLwing*Swing; 

ine ri z—q* CLhoriz* SnoriZ; 

Lvert=q*CLvert*Svert; 

Lftotal=Lwing+Lhoriz; 

alphaT=atan2 ((Dftotal+Drotor), (GW-Lftotal) ); 
mu=Vinf*cos (alphaT) /Vtip; 


% *** thrust, calculation *** 
af Vint < 6.9 
T=(1+(0.4*Afv/Adisk) )*GW; 


CT=T/ (Adusk*rho*Vtip 2) ; 
else 
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T= (GW-Lftotal) /cos(alphaT) ; 
T=T/ (Adisk* rho*Vtip%’2) ; 


end 

% *** setup blade radius elements, azimuth elements, 

% induced velocity distributions, and determination 
% of coning angle and tip loss parameter *** 


B=1-(sqrt(2*CT) /b) ; 

Reff=B*R; 

Rbar=Reff-e; 

dr=(Reff-grip) /nbe; 

r=grip:dr:Reff-dr;, r=rt+dr/2; 

Erl-0. 77% *** “first quess at err (=~ 
RbarT=rT1*Rbar; 

mblade=wblade/32.17; 

betao=asin ( (T/b*RbarT- (.5* (R-e) +e) *wblade) /((.5* (R- 
e)te)*2*omega*2*mblade) ); 

betat=twist* (0.7-(r/R)); 

psi=0: 360/naz: 360-360/naz;,psi=psi'/57.3; 


$$ set up vector of blade element chords and then varies them as 
$$ requested with the blade taper and blade taper start position 
$% rchord=root chord 

$% cbhlade=vector of blade element chord lengths 

$% tr=taper ratio (tip/root) 

$$ trst=taper ratio start position (r/R) 


cblade=rchord*ones(size(r)); % gives all elements same chord length 
ina tla lay 


if tr==0 % prevents division by zero later in code 


tr=1; % in case 0 is enter for taper ratio instead 
end % ‘of 1 £or no taper 
if trst== 
slope=(rchord-rchord*tr)/(Reff-grip); % Modifies each element 
cblade=cblade-slope*(r-grip); % chord length wrt input 
tchord=cbhlade (nbe); $ taper ratio which has been 
mchord=sum(cblade) /nbe; %$ been converted into a slope 
$ top portion takes into 
else % account the possibility that 
slope=(rchord-rchord*tr)/(R*(1l-trst)); %$ a O start position is 
really at 
Z=fix(nbe*trst) ; $ the start of the aero portion 
if z<=1 $ prevents beginning index fm being zero 
Zz=1; 
end 


cbhlade(z:nbe)=cblade(z:nbe)-(r(z:nbe)-xr(z))*slope; 
tchord=cblade (nbe); 
mchord=sum(cblade) /nbe; 

end 
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3 e#x* induced veli@ertey determinatwzen, “=~ 


if Vint <2.6- 9" 
A=4*pi; 
Bv=(b/2) *omega*a. *cblade; 
Tv=0; 
delT=T-Tv; 
while abs(delT) > .O1*T $ Prouty Eqns for Hover 
thetav=twist* (0.7-(r/R))+thetao; 
C=(-b/2) .*cblade*omega*2.*r*a.*thetav; 
Vi=(sBvVtsqrt (BV. 2-(4*A*C) ga a 
AT v= (b/2) * rho* ( (omega*r) .*2) *a.* (thetav- 
(vij3/ (omeqa*r) ))'.*eblade*dx; 
Tv=sum(dTv) ; 
delT=T-Tv; 
ne dell <70, 
thetao=thetao-0.5*thetao*abs(delT/T) ; 


else 
thetao=thetaot+0.5*thetao*abs (delT/T) ; 
end 
end 
else $ Wheatley Eqn for Fwd flt 
lamdaT=0; 
lamda=1; 


while abs (lamdaT-lamda)>le-4 
lamda=lamdaT; 
lamdaT=mu* sin (alphaT) +0.5*CT/sqrt (lamdaT*2+mu%2) ; 
end 
vi=lamdaT*Vtip-Vinf*sin(alphaT) ; 
Vi=vi*ones (size(r)); 
end 


% 77* first guess at theta *** 


thetalc=0.035* ( (0. 0006e-3*Vinf*2+0.244e-3*Vinf) /0.105) ; 
thetals=-0.087* ((0.0006e-3*Vinf%*%2+0.244e-3*Vinf)/0.105); 
theta=thetaotthetalc.*cos(psi)+thetals.*sin(psi) ; 


3 ~*~ EOLOL tCLEImmInNgGm@EOUutine *** 


Set(H STATUS, "String', ’TRIMMING COLLECTIVE' ) 

set (H STATUS], ’String’,[*Run Elapsed Time is * num2str(£fix(toc)) @ 
seconds']) 

pause (3) 


k=1; 
error0=(T*.02)+1; 


While abs(error0) > T*.02 
Tpsi=zeros (size(psi)); 
thireale 
Pikes, % Eccles change: These three lines were added. 


errorl; 


15] 


end 


error0=T-— (mean (Tpsi) *b); 
Pure rrorOsce—1* 02. 

thetao=thetao-0.35*thetao*abs (1.5*error0/T) * (1-mu) ; 
elseif error0 > T*.02Z, 

thetao=thetaot0.35*thetao*abs (1.5*error0/T) * (1-mu) ; 


end 
theta=thetaotthetalc.*cos (psi)+tthetals.*sin(psi) ; 
feek > 1, 
if abs(error0O) > abs(errorl), 
ele 
trim warning 
Seti GO, Enable*, off"); 
Sec(H RES, Enable”, off"); 
Scr(HeRUrl, Bnasle Off); 
set(H BK, Enable , on); 
Crror( 9+. = pigeNnoe Tram ***") 
end 
end 
errori=error0O; 
k=k+1; 


end 


set(H STATUS, *String', ‘TREMMING CYCLIC") 
Sser(heslAlUst string @aahuneelapsed Time 1s " numZser(fix(toc)) 
seconds']}) 

pause (3) 


tO0=clock; 
k=1; 
errord—( (1/b)*rlis {R-grip) )*,04) +1; 


Witte erroroe> ((lfb) rt la(R-grip) )*.04 
time=etime (clock, tO) ; 
if time > 15, 
set (H_STATUS, 'String', 'STILL TRIMMING ...") 
See(iectnluol, string , |" Run Blapsed Time 1s | 
NumZStret fix (coc) j=: seconds" | } 
pause (3) 
t0=clock; 
end 


Mpsi(:,k)=zeros(size(psi)); 
tmcalc 

theta=[theta theta(:,k)]; 
Mpsi=[Mpsi Mpsi(:,k)]; 


% *** calcullation of initial dthetadM *** 
if k <2, 
theta (:,k+1)=theta(:,k)+0.25/57.3; 


Mpsi(:,kt1)=zeros(size(psi)); 
k=k+1; 
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tmeale 

k=k-1; 

dthetadM= (theta (:,kt1})—theta{:,k)) -/ (Mpsi(:,k+1)—Mpsi(:,k)); 
end 


% *** calculation of M first harmonic parameters *** 


Mlc=2*sum(Mpsi(:,k).*cos (psi) )/naz; 
Mls=2*sum(Mpsi(:,k).*sin(psi) )/naz; 


% *** removal of first harmonic terms from Mpsi *** 


Mpsi(:,k+1)=Mpsi(:,k)-Mlc.*cos (psi)-Mls.*sin(psi) ; 
delM=Mpsi (:,k+1)-Mpsi(:,k); 
errorQ=max (delM) -min(delM) ; 
Lek Sele 
if error0 > error! 
eLe 
trim warning 
Secico,. Enable ;. off") 
set (RES, 'Enable','off'); 
set (RUPT, Enable’ ,;"off'); 
set (BK, 'Enable','on'); 
error('*** END OF PROGRAM ***') 
end 
end 
errorl=errorQ; 


% *** calculation of new theta *** 


de1lM=0.5* (1-mu) *delM; 
theta(:,k+1)=theta (:,k)+(dthetadM. *delM) ; 


if errorO <= ((T/b)*rT1*(R-grip))*.04 
thetalc=2*sum(theta(:,k).*cos(psi))/naz; 
thetals=2*sum(theta(:,k).*sin(psi))/naz; 
else 
thetalc=2*sum(theta(:,k+1).*cos(psi))/naz; 
thetals=2*sum(theta(:,k+1).*sin(psi))/naz; 
end 


theta(:,k+1)=thetao+tthetalc.*cos(psi)+thetals.*sin(psi); 
% *** calculation of new dthetadM *** 


theta=(theta theta(:,k+1)]}; 

Mpsi=(Mpsi Mpsi(:,k+1)]}; 

theta(:,k+2)=theta(:,k)+0.25/57.3; 
Mpsi(:,k+2)=zeros(size(Mpsi(:,k+1))); 

k=k+2; 

tmcalc 

k=k=-2; 
dthetadM=(theta(:,k+2)-theta(:,k))./(Mpsi(:,k+2)-Mpsi(:,k)); 
k=k+1; 
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end 


seu(heelAlUs, potting’), AVIUSTING COLLECTIVE’) 


SeGheor GUod, oerlng. 7,1 Run Elapsed Time is 
seconds']) 
pause (3) 


theta=theta(:,k); 
k=1; 
error0Q=(T* .01)+1; 


While abs(error0O) > T*.01 
Tpsi=zeros(size(psi)); 
thrcealc 
error0=T- (mean (Tpsi) *b) ; 
if errorO < -T*.0Ol1, 


Mums str (fix (coc) ) 


thetac—tnetaco-0.25*thetac*abs (1.25*erroro/T) * (1—mu); 


elseif error0O > T*.0O1, 


thetao=thetaot0.25*thetao*abs (1.25*error0/T) * (1—mw) ; 


end 
theta=thetao+tthetalc.*cos (psi)+thetals.*sin(psi) ; 
Pevke> 1, 
if abs(error0) > abs(errorl), 
eke 


trim warning 
error('*** END OF PROGRAM ***') 
end 
end 
errorl=error0; 
k=k+1; 
end 


= ***Scaleulating »arag moments *** 


DMpsi=zeros(size(psi)); 
dmcalc 


% =2* calculating focor H force *** 


Te Vink <a lono, 
Hrotor=0; 
aT=[(dT ddT]; 
dD=[dD ddD]; 
else 
aT=({daT ddT]; 
dD=[{dD ddD]; 
for i=l:length(r)+1, 
leo * sumi(anh |: , le ~Cosipsimy naz 
Hie (i) =24 umd +e ssn psas))) 7 naz 


end 
Hrotor=(((b*cos(alphaT) /2) * (sum(H1s) - 
sin(betao) *sum(Hlc) ))+Drotor)/2; 


end 
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% *** calculating=- new eae 

rT2=(( (mean (Mpsi(:,length(Mpsi(1,:))-1))/mean(Tpsi) )/R)+rT1)/2; 
% *** check Lotor drag and rT, retrim roton 1f teoquiredun.. 
while abs (DEOtOE—-HrEOotor) >= US2*Hrotoru| abs (rill —-r lz) ee le ere 


if abs (Drotor—Hrotor).> 0.2*hietor, 
set (H_STATUS, 'String', ‘ADJUSTING ROTOR DRAG') 
set(H STATUS], ‘String', | Run Elapsed Game is 
Mnum2sStr (fix (Oc) es seconds |) 
pause (3) 
end 
Drotor=h rotor, 
Lf absi(erl—riZ ee 0. OSs rT, 
set (H_STATUS, 'String', ‘ADJUSTING MEAN THRUST LOCATION ') 
set (H STATUS], “string 7 il Run Elapsed Time is‘ 
num2str(fix(toc)) ' seconds'])) 
pause (2) 
end 


set (H_ STATUS, 'String', 'RETRIMMING ROTOR') 

set(H STATUSI, 'String',[‘Run Elapsed Time is ' num2str(fix(toc)) ' 
seconds']) 

pause (3) 


aT=daT(:,1l:nbe); 
aD=dD(:,l1l:nbe); 


% *** recalculating parameters *** 


alphaT=atan ((Dftotal+Drotor) / (GW-Lftotal) ); 
mu=Vinf*cos (alphaT) /Vtip; 


Le Vinive—sle. 9, % Wheatley Eqn for Fwd flt 
T= (GW-Lftotal) /cos(alphaT) ; 
CT=T/ (Adisk* rho*Vtip%2) ; 
lamdaT=0; 
lamda=1; 
while abs (lamdaT-lamda) >le-4 
lamda=lamdaT; 
lamdaT=mu*sin(alphaT)+0.5*CT/sgqrt (lamdaT*2+mu%2) ; 
end 
vi=lamdaT*Vtip-Vinf*sin(alphaT) ; 
vi=vitones (size(r)); 
end 


B=1- (sqrt (2*CT)/b); 
Reff=B*R; 

Rbar=Reff-e; 
dr=(Reff-grip) /nbe; 
r=grip:dr:Reff-dr;, r=rtdr/2; 
RbarT=rT2* Rbar; 
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betao=asin( (T/b*RbarT-(.5* (R-e) te) *wblade) /((.5* (R- 
e) +e) *2*omega*2*mblade) ); 


% - a otoimeng collective *** 
tO0=clock; 
k=1; 


error0=(T*.02)+1; 


while abs(errorQ) > T*.02 
Tpsi=zeros (size(psi)); 
cthreaic 
error0=T- (mean (Tpsi) *b); 
Mier rord << .—T*=.02, 
thetao=thetao-0.35*thetao*abs (1.5*error0O/T) * (i-mu); 
elseif errorO > T*.02, 
thetao=thetaot+0.35*thetao*abs (1.5*error0O/T) * (1-mu) ; 
end 
theta=thetaotthetalc.*cos(psi)+thetals.*sin(psi); 
at * eee ihe 
if abs(error0) > abs(errorl), 
ele 
trim _warning 
error('*** END OF PROGRAM ***') 
end 
end 
errori=error0; 
k=k+1; 
end 


% Peete LU CaenyiG iG: = a 


k=1; 
errorQ=( ((T/b) *rT2* (R-grip) )*.04)+1; 


Wile error”d S(T /biert2*(R-grip)) > .04 
time=etime (clock,t0); 
at time > 15, 


Scie olaiwwo ering, , Still TRIMMING 227.) 

set (H STATUS], 'String', ['Run Elapsed Time is ' 
num2str(fix(toc)) ' seconds']}) 

pause (3) 

tO0=clock; 


end 
Mpsi(:,k)=zeros(size(psi)); 
tmcalc 

theta=[theta theta(:,k)]; 
Mpsi—[Mpsa Mpsi(:,k)1; 


% [re awCaleUlattOn OL Initial dthetadM *** 


fee 2, 
theta(:,k+1)=theta(:,k)+0.25/57.3; 
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Mpsi (27a) —zeros(size (psi) 
k=k+1; 
tmcalc 
k=k-1; 
dthetadM= (theta (:,k+1)-theta(:,k))./(Mpsi(:,k+1)- 
Mp sa iC gk) es, 
end 


% *** calculation of M first harmonic parameters *** 


Mlc=2*sum(Mpsi(:,k).*cos(psi))/naz; 
Mls=2*sum(Mpsi(:,k).*sin(psi))/naz; 


% *** removal of first harmonic terms from Mpsi *** 


Mpsi (:,k+1)=Mpsi(:,k)-Mlc.*cos (psi) -Mls.*sin(psi) ; 
delM=Mpsi (:,k+1)-Mpsi(:,k); 
error0=max (delM) -min(delM) ; 


i) age ee ae 
if CrEror0 > Greer, 
ellie 
trim warning 
error('*** END OF PROGRAM ***') 
end 
end 


errorl=error0; 
% *** calculation of new theta *** 


delM=0.5* (1-mu) *delM; 

theta(:,k+1)=theta(:,k)+(dthetadM.*delM) ; 

if error0 <= ((T/b) *rT2* (R-grip)}*.04 
thetalc=2*sum(theta(:,k) .*cos(psi) )/naz; 
thetals=2*sum(theta(:,k).*sin(psi) )/naz; 

else 
thetalc=2*sum(theta(:,k+1).*cos (psi) )/naz; 
thetals=2*sum(theta(:,k+1).*sin(psi) )/naz; 

end 

theta(:,k+1)=thetaotthetalc.*cos (psi) tthetals.*sin(psi) ; 


% *** calculation of new dthetadM *** 


theta=[theta theta(:,k+l)]; 

Mpsi=[(Mpsi Mpsi(:,k+1)]; 

theta(:,k+2)=theta(:,k)+0.25/57.3; 

Mpsi (:,k+2)=zeros(size(Mpsi(:,k+l1))); 

k=k+2; 

tmcalec 

k=k-2; 
dthetadM=(theta(:,k+2)-theta(:,k))./(Mpsi(:,k+2)-Mpsi(:,k)); 
k=k+1; 


ey, 


end 
% aL Cmetninind  CcoLlective *** 


theta=theta(:,k); 
k=1; 
error0=(T*.01)+1; 


While abs(errorQ) > T*.Q1 

Tpsi=zeros (size(psi)); 

threale 

error0=T- (mean(Tpsi) *b); 

af G€rror0. <= —17.0l, 
thetao=thetao-0.25*thetao*abs (1.25*error0/T) * (1-mu) ; 

elseif error0O > T*.0Ol, 
thetao=thetaot0.25*thetao*abs(1.25*error0/T) * (1-mu) ; 


end 
theta=thetao+thetalc.*cos(psi)+thetals.*sin(psi) ; 
rod, p< meee 
if abs(errorQ) > abs(errorl), 
(al We. 
trim warning 
error(****~ END OF PROGRAM ***"*) 
end 
end 
errorl=error0; 
k=k+1; 
end 
% “At “recalculating szotor H force *** 


Peavant« << 16.9. 
Hrotor=0; 
aT=[(dT ddT]; 
aD=[dD ddD}; 
else 
adT=[dT ddT]; 
aD=[dD ddD]; 
for i=l:length(r)+l, 
Hlc(i)=2*sum(dT(:,i).*cos(psi) )/naz; 
Hls(i)=2*sum(dD(:,i).*sin(psi))/naz; 
end 
Hrotor=(((b*cos (alphaT) /2) * (sum(H1s) - 
sin (betao) *sum(H1c) ))+Drotor) /2; 


end 

% amen caleulating rl *** 

ri l=riZ; 
rT2=(((mean(Mpsi(:,length(Mpsi(1,:))-1))/mean(Tpsi) )/R)+rT1)/2; 


end 


Se) ** recalenulacing drag moments *** 
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al=di(: -lenbeie 
dD=dD(:,1l:nbe); 
DMpsi=zeros (size (psi) ); 
dmcalec 

di=(ar dari: 

aD=[{dD ddD]; 


set(H STATUS, 'String', 'ROTOR TRIMMED') 

set(H STAs sering,, . hun slapsed: time 45° uumzst Yr (fix(tece) 
seconds ']) 

pause (3) 


save perftemp mchord DMpsi % Eccles addition - perf.m was 
$ not recognizing mchord and DMpsi. 
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APPENDIX AH. THRCALC.M 


This script M-file is a subroutine of Trim.m to calculate the rotor thrust. 


Threale.™ 

calculates the total thrust along a blade at 
each azimuth (psi) location 

JANRAD 98 VERSION 4.0 


oP oP oP oP 


Up=zeros (size(psi*r) ); 
Ut=zeros(size(Up)); 
daT=zeros(size(Up)); 
ddT=zeros(size(psi)); 


for i=l:length(psi), 


Up, :) =v. -cos (betac) +Vinti*sin(alphal) *cos(betao) -vint*cos(alphal) *caim 
(betao) *cos (psi(i)); 

Ut (1,;)=r. *omegatVinf*cos (alphaT) *sin(psi(i)); 
Ehi=atanZ(Upii,:), Uiteli, =) ); 

alpha=theta (i)+betat-phi; 


% Eccles added the following line for use with Mach dependent afoil 
pales. 
Mach = 
(Vtip.*cos(alphaT).*r./RtVinf.*sin(psi(i)))/(49.05*sqrt (temp+460) ); 


if afoil==1, 
{[CL,CD]=ool2clcd(alpha, Mach); 
elseif afoil==2, 
{CL,CD]=hh0O2clced (alpha) ; 
elseif afoil==3, 
(CL, CD] =vrl2clcd(alpha); 
elseif afoil==4, 
[CL,CD]=scl094r8clcd(alpha,Mach) ; 
elseif afoil==5, 
{[CL, CD] =scl095r8clcd(alpha,Mach) ; 
end 


lla, ) —Uneo? Pho. cColadetdr.=(Upiae).  2HUE(I, =). 2) (Cia COs (pings 
eb > sin(pha) ) >; 
Tpsi (i)=sum(dT(i,:)); 


=) ~** cCalculataioms for tip Loss varea ~~ 
Uptip=Vinf*sin (alphaT) *cos (betao) +Vinf*cos (alphaT) *sin (betao) *cos (psi (i 
)); 


Uttip=(R- (R-Reff)/2) *omega+Vinf*cos(alphaT) *sin(psi(i)); 
phitip=atan2 (Uptip, Uttip) ; 


16] 


oii Oo Liozem lade (noe) (0,510.5*cos (2*psi (1)))* (R- 
Refeieeuiot ip 2tUtea 2) > {—.009*sin(phitip) ) ; 

Tpsi (1) =Tpsi (i)+ddT (i); 

end 
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APPENDIX AI. TMCALC.M 


This script M-file is a subroutine of Trim.m to calculate the rotor thrust moment. 


fucaLesm 

calculates the total thrust moment along a blade 
at each azimuth (psi) location 

JANRAD 98 VERSION 4.0 


oP om OP oP 


Up=zeros (size(psi*r) ); 
Ut=zeros (size (Up) ); 
dM=zeros(size(Up)); 
ddM=zeros(size(psi) ); 


for i=l:length(psi), 


Up(i,:)=vi.*cos (betao) +Vinf*sin(alphaT) *cos (betao) +Vinf*cos (alphaT) *sin 
(petac)* cos (psi. (7))); 
Ut (i,:)=r.*omegat+Vinf*cos (alphaT) *sin(psi(i)); 


Bha-atanc(Upia, +) Ut (a, :) )s 
alpha=theta(i,k)+betat-phi; 


% Eccles added the following line for use with Mach dependent afoil 
files. 
Mach = 

(Vtip.*cos (alphaT) .*r./R+Vinf.*sin(psi(i)))/(49.05*sqrt (tempt+460) ); 


if afoil==1, 
[CL,CD]=ool2clced(alpha, Mach); 
elseif afoil==2, 
(CL, CD}]=hh02clcd (alpha) ; 
elseif afoil==3, 
[CL,CDJ=vrl2clcd(alpha) ; 
elseif afoil==4, 
[CL, CD]=sc1094r8clcd(alpha,Mach) ; 
elseif afoil==5, 
(CL, CDJ=scl1095r8clcd(alpha, Mach) ; 
end 


aM: )—O.0 Eno, colade.* 1 dr. (Up (i, :).°Z2tUGi(i,:). 2) -7 (el. “eos (phi — 
eben sin (phi); 
Mpsi(i,k)=sum(dM(i,:)); 


te. *** calculations fer tip loss areas **=* 
Uptip=Vinf*sin (alphaT) *cos (betao)+Vinf*cos (alphaT) *sin(betao) *cos (psi(i 
IG 


Uttip=(R- (R-Reff) /2) *omega+Vinf*cos (alphaT) *sin(psi(i)); 
phitip=atan2 (Uptip, Uttip) ; 
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adM(i)=0.5*rho*cblade (nbe) * (0.5+0.5*cos (2*psi (i) ) ) * (R- (R-Reff) /2) * (R- 
Romie Uo 2rUtcip 2) ~(=.009*sin (phitip) ); 

Mpsi (i,k) =Mpsi (i, k) +ddM (i); 

end 
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APPENDIX AJ. DMCALC.M 


This script M-file is a subroutine of Trim.m to calculate the rotor drag moment. 


DMCALC .M 

dmcalc calculates the total drag along a blade at 
each azimuth (psi) location. Called from Trim.m 
JANRAD 98 VERSION 4.0 


oP oP oP of 


Up=zeros (size(psi*r)); 
Ut=zeros(size(Up)); 
alpham=zeros (size(Up) ); 
dD=zeros(size(Up)); 
ddD=zeros (size(psi) ); 
ddDM=zeros(size(psi) ); 


for i=l:length(psi), 


Up(i,:)=vi.*cos (betao) +Vinf*sin (alphaT) *cos (betao)+Vinf*cos (alphaT) *sin 
(betao) *cos (psi(i)); 
Ut (i, :)=r.*omegat+tVinf*cos (alphaT) *sin(psi(i)); 


pol—dtanz(Up(1,:),Utl(i,:)); 
alpha=theta(i)+betat-phi; 
alpham(i,:)=alpha; 


% Eccles added the following line for use with Mach dependent afoil 
files. 
Mach = 


(Vtip.*cos(alphaT) .*r./R+Vinf. *sin(psi(i))}7 (49. 05*sqrt (temp+460) ) ; 


if afoil==1, 

[CL,CD]=ool2clcd(alpha, Mach); 
elseif afoil==2, 

[CL, CD] =hhO2clcd(alpha) ; 
elseif afoil==3, 

[CL, CD) =vrl2clcd(alpha) ; 
elseif afoil==4, 

[CL, CD) =scl094r8clcd(alpha, Mach) ; 
elseif afoil==5, 

[CL, CD] =sc1095r8clcd(alpha, Mach) ; 
end 


apis: )=0.5* rho. seblade*dr.* (Upti,. jn 2 UE (,: ). 2) (Cb. sin (phi) +epe= 
cos (phi) ); 

GDM=dD(1i,:)."r> 

DMpsi (1)=sum(dDM) ; 


$ *** calculations for tip loss area *** 
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Uptip=Vinf*sin(alphaT) *cos (betao) +Vinf*cos (alphaT) *sin(betao) *cos(psi(i 
))3 

Uttip=(R- (R-Reff) /2) *omegat+Vinf*cos (alphaT) *sin(psi(i)); 

phitip=atan2 (Uptip,Uttip) ; 


adi) =O" >> rho~ecblade (nbe)* (0.54+0.5*cos (2*psi(i)))*(R- 
Rete sOptip 2+ULtIp 2) “(2009*%cos (phi tip).); 

ddDM (i) =ddD (i) * (R-(R-Reff) /2); 

DMpsi (i) =DMpsi (1) +ddDM(i); 
end 
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APPENDIX AK. 0012CLCD.M 


This MATLAB function calculates a c, and cg for a NACA 0012 airfoil given 


angle of attack and Mach number. It is called in Thrcalc.m, Tmcalc.m and Dmcalc.m 


script M-files. 


Both 


3 
3 
3 
3 
3 
3 
% ‘alpha' 
3 
3 
3 
3 


function [CL,CD]=oolZclcd(alpha, 


[CL, CD]=ool2clcd(alpha, Mach) 


'"Mach' 


CL=zeros (size(alpha) ); 
CD=zeros (size(alpha) ); 


a=alpha*180/pi; 
aL = 
aD = 
Kl 
K2 


‘J MOF 


chk= (Mach<0.725 


SE=CL+chk.* ((0. 


chk=(Mach<0.725 


GL-—CL+chk.*((0. 


chk=(Mach<0.725 


eiL—-CL-cenk.*( (0. 


aoe. Ke) 


chk=(Mach<0.725 


CL=CL-chk.*((0. 


15> — 16. “Mach; 

ie e2 5.4. Mach; 

O@20233 + -0.342.* (Mach.*7.15); 
2205 -— 0.95." Mach; 


a>=0 & a<=aL); 


./sqrt(1l-Mach.%2) 


a>aL & a<=20); 


./sqrt(1-Mach.%2) 


a>=-20 & a<-aL); 


ey sqre(l-Mach. 2) 


a>=-aL & a<Q); 


ww sqrtetl-Mach. 2) 


Mach) 


“01. 


for Mach numbers < 0.725 and AOA inside +/- 


*Mach) 


.*Mach) 


. *Mach) 


. *Mach) 


%$ CL for Mach numbers > 0.725 and AOA inside +/- 


chk=(Mach>=0.725 & a>=0 & a<=aL); 
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.*abs (a) 


ool2clcd calculates CL and CD for the NACA 0012 
airfoil given angle of attack in radians and the 
local Mach number: 


are intended to be vectors 

the elements of which correspond to the rotor blade 
radial stations of interest in a blade element analysis. 
All equations are based on Ray Prouty's treatment of 
the 0012 in his text. 


20 deg: 


YS culics 


-*a - K1.* (a-aL) .*K2); 


- K1l.* (abs (a)- 


-*abs(a)); 


20 deg: 


Cr—-@Clmchik. = (1.0.07 7 —§02744.*Mach).*a) ; 


chk=(Mach>=0.725 & a>aL & a<=20); 
Cl-eracnk. W000? ——~60.)44.*Mach).*a — (0.0575-0.144.* (Mach- 
On Zone 4452" (a-—abL) 27 (K2) ) ¢ 


chk=(Mach>=0.725 & a<O & a>=-aL); 
CL=CL-chk.*((0.677 - 0.744.*Mach) .*abs(a)); 


chk=(Mach>=0.725 & a<-aL & a>=-20); 
@b-CGlo—chik., = ((0.677 107744. "Mach) .~abs (a) — (0.0575-@.144.* (Mach- 
0725). 0.44) .* (abs (ajieab) .~ (K2)au; 


$ CL for all Mach numbers and AOA outside +/- 20deg: 


chk=(a>20 & a<=161); 
CL=CL+ehk. * (1.15. *sami2iealpha) ); 


chk=(a>161 & a<=173); 
Ci—Clachnk. ~ (0. 7 )e; 


chk=(a>173 & a<=180); 
Cl=—Cilaecnrk.* (0.1.* (a—-130) ); 


chk=(a>=-180 & a<=-173); 
Ch=ChL+echk.*(0.1.* (att 180) ); 


chk=(a>-173 & a<=-161); 
CL=Chtenk.~ (0.7): 


chk=(a>-161 & a<-20); 
ClL=Citehnk.* (1. 15. *san(2.*alpha) }; 


%$ CD for Mach numbers < 0.725 and AOA inside +/- 20 deg: 


chk=(Mach<0.725 & a>=0 & a<=aD); 
CD=CD+chk.* (0.0081 + (-350.*a + 396.*a.%2 - 63.3.*a.%3 + 
Be OOm As) 42 LO. (= 6):).3 


chk=(Mach<0.725 & a>aD & a<=20); 
CD=CD+chk.*((0.008]1 + (-350.*a + 396.*a.*%2 - 63.3.%*a.%3 + 
3.66.*a.*%4).*10.%(-6)) + 0.00066.* (a-aD) .*2.54); 


chk=(Mach<0.725 & a<0O & a>=-aD); 
CP—-CH7chky (Oe008 1 + (—-350.*abs(a) + 396.*a.°2 - 63:3.*abs(a).°3 + 
3.06.7 a. 4) .* 10.>(=—6) )-: 


chk=(Mach<0.725 & a<-aD & a>=-20); 


GD=E€Drechke 1 (00-0001) + (-350.*abs(a) + 396.%a.°2 — 6323. *abs (a) .73 + 
BOcn awed 10. (—-o)) +.0.0006G, ~(abs(a)—aD).°2.54) ; 
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% CD for Mach numbers > 0.725 and AOA inside +/— 20 deg: 


chk=(Mach>=0.725 & a>=0 & a<=20); 
CD=CD+chkh eer m0 00s  (—-350.*aues96.*ai22e— Come asc, 
3.66.*a.°4).* 10." (=6)) pt 0.00035.*a. 02-54 .+ 2 (Mach-0. /25)eeo sae 


chk=(Mach>=0.725 & a<O & a>=-20); 

CD=CD+chk.*((0.0081 + (-350.*abs(a) + 396.*a.*2 - 63.3.*abs(a).*3 + 
3.66.*a.7*4). *10.*(-G)) + 0.00035. *absia).°2.54 + 21.* (Mach— 

Dr 2 ON eee 

% CD for all Mach numbers and AOA outside +/- 20deg: 


chk=(a>20 & a<=180); 
SU-Cbtenk. = (1 0Se-)1.02.*cos(2.*alpha) ); 


chk=(a>=-180 & a<-20); 
eP-CDtehnk.~ (1.03 = 91.02.*cos(2.*alpha) ); 
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APPENDIX AL. HH02CLCD.M 


This MATLAB function calculates a c, and cy for an HH-02 airfoil given angle of 


attack. It is called in Thrcalc.m, Tmcalc.m and Dmcalc.m script M-files. 


function [CL,CD]=hh02clcd (alpha) 


$ hhO2clcd calculates CL and CD for an HH-02 airfoil 
%$ given angle of attack (alpha) in radians 

$ [CL,CD]=hh0O2clcd (alpha) 
% JANRAD 98 VERSION 4.0 


CL=zeros (size (alpha) ); 
CD=zeros (size(alpha)); 
a=alpha*180/pi; 


chkl=(a>=20 & a<=180); 

CL=CL+chk1.* (0.425414+0.026863*a+5.5988e-4%*a.%2-2.1493e-5*a.*%341.5932e- 
T*a.*4-3.4659e-10*a.%5); 

CD=CD+chk1. * (-0.7179+0.061213*a-5.9861le-4*a.%2+7.3708e-6%*a.%*3-6.6605e- 
8*a.*44+1.913e-10*a.%5); 


chkl1=(a>=-180 & a<=-50); 

CL=CL+chkl.* (-4.6183-0.1923*a-3.5554e-3*a.*%2-3.3273e-5*a.*%3-1.4528e- 
T*a.*4-2.3003e-10*a.%5); 

CD=CD+chk1l.* (2.7093e-2-2.1309e-2*at2.0335e-4%a.*°%2+3.47e-T*a.*3- 
3.0586e-8*a.*4-1.2584e-10*a.%5); 


chkl=(a>-50 & a<-20); 

CL=CL+chkl.* (-2.5519-0.22847*a-9.5667e-3*a.*2-1.705le-4*a.*3-1.0909e- 
Ora. 4) > 

CD=CD+chk1.* (2.7093e-2-2.1309e-2*at+2.0335e-4*a.*%2+3.47e-7*a.*3- 
3.0586e-8*a.*4-1.2584e-10*a.%5); 


chk1=(a>=-20 & a<=-10); 

CL=CL+chk1l.* (-0.2+0.089*at+0.0034*a.%2); 

CD=CD+chk1.* (2.7093e-2-2.1309e-2*at2.0335e-4*a.%2+3.47e-7*a.%3- 
3.0586e-8%a.%4-1.2584e-10*a.%5) ; 


chkl=(a<20 & a>-10); 

CL=CLt+chkl.* (5.8766e-24+1.313le-1l*at2.4742e-3%*a.*%2-5.303e-4*a.%*3- 
1.5818e-5*a.*%4+1.28e-6*a.%5); 

chk2=a<-4; 

chk2=chk2.*chkl1; 

CD=CD+chk2 .* (1.3786+0.916*a+0.21396*a.*%2+2.0371le-2*a.*3+7.0076e- 
Aa A 

chk2=(a>=-4 & a<=7); 

chk2=chk2.*chkl; 
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CD=CD+chk2.* (9.732e-34+3.2326e-4* atl 
97a. 341.1826e-6*a.*44+1.527le-6*a.%*5); 

chk2=a>7; 

chk2=chk2.*chkl; 

CD=CD+chk2.* (1.842e-1-5.7532e-2*at5 


M104 


-4392e-4*a.*2-8.5073e- 


BOA Se [So cae Zot 2oOSse—4*a.*3); 


APPENDIX AM. VRI2CLCD.M 


This MATLAB function calculates a c) and cg for an VR-12 airfoil given angle of 


attack. It is called in Thrcalc.m, Tmcalc.m and Dmcalc.m script M-files. 


function [CL,CD]=vrl2clcd (alpha) 


vrl2clcd calculates CL and CD for the VR-12 airfoil 
given angle of attack (alpha) in radians 

[CL, CD] =vrl2clcd(alpha) 

JANRAD 98 VERSION 4.0 


of oP HH oP 


CL=zeros(size(alpha)); 
CD=zeros (size(alpha)); 
a=alpha*180/pi; 


chk=(a>=20 & a<=180); 
CL=CL+chk.* (1.1733-0.018879*a+1.5762e-3*a.*2-3.1925e-5*a.%*342.0949e- 
T*a.*%4-4.3807e-10*a.%5); 


chk=(a>=-180 & a<=-50); 
CL=CLt+chk.* (-4.6183-0.1923*a-3.5554e-3*a.*%2-3.3273e-5*a.%*3-1.4528e- 
7*a.*%4-2.3003e-10*a.%5); 


chk=(a>-50 & a<-30); 
CL=CL+chk. * (-0.22114+0.020857*a+2. 857le-4*a.%2) ; 


chk=(a>=-30 & a<=-10); 
CL=CLt+chk.* (-1.11-0.12383*a-0.01515*a.*2-6. 8667e-4*a.*3-le-5*a.%4) ; 


chk=(a<20 & a>-10); | 
CL=CL+chk.* (0.11976+0.12341*a+t5. 5841e- 4*a.*2-2.0652e-4*a.%3); 


chk=(a>=17 & a<=180); 
CD=CD+chk. * (-0.263764+0.017917*at6.9927e-4*a.*2-9.1137e-6%*a.%34+2.6277e- 
Sra. 4); 


chk=(a>=-180 & a<=-10); 
Rae *(-0.17486-0.034463*a-1. 0233e—-4*a.*2-2.8958e-6*a. > 3-4, 657 7e- 
4-1 oo e-10*a. 75) 


chk=(a>-10 & a<=0); 
CD=CD+chk.* (9.8678e-34+3.4934e-3*a+1.4844e-3*a.%2-1.3564e-4*a.%3- 
1.0936e-5*a.%4); 


chk=(a>0O0 & a<=15); 


CD=CD+chk.* (9.8e-3+7.0457e-4*at+5.6104e-5*a.%*%2-4.115le-5*a.%*3+3.8695e- 
St aac4a): 
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enk=(a>l5 & a<17); 
CDC) peni. (te ootl., ocoe-|*~a-Z.9e-3*a.*Z); 
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APPENDIX AN. SC1094R8CLCD.M 


This MATLAB function calculates a c and cg for an Sikorsky SC1094R8 airfoil 
given angle of attack and Mach number. It is called in Thrcalc.m, Tmcalc.m and 


Dmcalc.m script M-files. 


function [CL,CD])=Sc1094r8clcd (alpha, Mach) 


$ scl094r8 calculates CL and CD for a Sikorsky SC1094R8 airfoil 
% given angle of attack (alpha) in radians and Mach number (Mach) 
$ ([CL,CD}=sc1094r8clcd(alpha, Mach) 


CL=zeros (size(alpha) ); 
CD=zeros (size(alpha)); 
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for j = 1:size (alpha) 
Pe alpnal: Jj pl 
alpha(:,}) = alpha(:,j) + 2*pi; 
elseif alpha(:,j) > pi 
alpha(:,j) = alpha(:,j) - 2*pi; 
end 
end 


Mach = abs(Mach) ; 
a=alpha.*180/pi; 


CL 
CD 


diag(table2(cl_tab,a,Mach))'; 
diag (table2(cd_tab,a,Mach))'; 
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APPENDIX AO. SC1095R8CLCD.M 


This MATLAB function calculates a c and cy for an Sikorsky SC1095R8 airfoil 
given angle of attack and Mach number. It is called in Thrcalc.m, Tmcalc.m and 


Dmcalc.m script M-files. 


function [CL,CD]=Scl1095r8clcd (ALPHA, MACH) 


%$ scl095r8 calculates CL and CD for a Sikorsky SC1O95R8 airfoil 
$ given angle of attack (ALPHA) in radians and Mach number (mach) 
$ [CL,CD)=scl1095r8clcd (ALPHA, MACH) 


CL=zeros (size (ALPHA) ); 
CD=zeros (size (ALPHA) ); 
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for j = 1:size(alpha) 
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alpha (:,j))e— alphat:, )) = 2*pi; 
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end 


Mach = abs (Mach) ; 
a=alpha.*180/pi; 


CL 
O19) 


diag(table2(cl_tab,a,Mach))'; 
diag (table2(cd_ tab,a,Mach))'; 
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APPENDIX AP. VARIABLE _LIST.M 


This script M-file contains a list of variables used in JANRAD 98 version 4.0. 
The list is broken up into four parts. The first part lists variables used primarily for 
computation and analysis. Then global variables, structure variables and GUI graphic 


handles are listed. 


% This File lists the Variables used in JANRAD 98 Version 4.0 


%$ Computational Vars. 


% a lift curve slope of rotor system airfoil 

% Adisk area of rotor disk 

% Afh fuselage equivalent flat plate drag area 

% Afv vertical projected area (fuselage area under disk) 
% aod 1. rotor system airfoil type (HH02/VR12) 

% alpha angle of attack, rotor blade radial segment 
% alphaT rotor tip path plane angle 

% b number of rotor blades 

% B tip loss parameter 

% betao rotor coning angle 

% betat geometric angle, rotor blade radial segment 
% bhoriz span, horizontal tail 

% bvert span, vertical tail 

% bwing span, Wing 

% cblade chord, rotor blade 

% CD drag coefficient, rotor blade radial segment 
% CDohorsz profile drag coefficient, horizontal tail 

% CDovert profile drag coefficient, vertical tail 

% CDowing profile drag coefficient, wing 

% CDhoriz drag coefficient, horizontal tail 

% CDvert drag coefficient, vertical tail 

% CDwing drag coefficient, wing 

% CH rotor H-force coefficient 

% CHesig CH/solidity 

% CL lift coefficient, rotor blade radial segment 
% CiEnor LZ lift coefficient, horizontal tam 

% CLvert lift coefficient, vertical tail 

3 CLwing lift coefficient, wing 

% CQ rotor torque coefficient 

% CQ sig CQ/solidity 

% Ch rotor thrust coefficient 

% Ci saig CT/solidity 

% dD differential drag, rotor blade radial segment 
% dadD differential drag, rotor blade tip 

% ddDM differential drag moment, rotor blade tip 

% ddM differential thrust moment, rotor blade tip 
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© 0 CH © H HH © © ® HL © HC H © HH H WH HH HH HH HH HH HH HH HH HH CH CH CH CO CO OO CP CH CM CMH CH CH CO CH CH CH CH CH CO CH CH CH CH CH OH Of 


Gar 
delM 
Dftotal 
Dfuse 
DL 

dM 
DMpsi 
ay 
Drotor 
ollh 
Dhoriz 


dthetadM 


Dvert 
Dwing 
e 

ewing 


filename 


EM 
grip 
GW 
HEOrOE 
lamdaT 
Lftotal 
LHOLLZ 
Lvert 
Lwing 
Mic 
Mls 
Machtip 
mblade 
Mpsi 
mu 

naz 
nbe 
omega 
PA 

phi 
phaeip 
Proton 
psi 

g 
Orotor 
E 

R 

Rbar 
RbarT 
Reff 
rho 

ledk 


solidity 


Shoriz 
Svert 
Swing 
it 


differential thrust, rotor blade tip 

change in total thrust moment 

resultant of fuselage drag and aux thrust 

total drag generated by non-rotor bodies 

disk loading 

differential thrust moment, rotor blade radial seg 
total blade drag moment at specific azimuth angle 
rotor blade radial segment width 

rotor system drag 

differential thrust, rotor blade radial segment 
drag, horizontal tail 

change in cyclic pitch with change in thrust moment 
drag, vertical tail 

drag, wing 

effective hinge offset 

wing efficiency factor 

name of input file 

figure of merit 

length of inner non-aerodynamic portion of blade 
aircraft gross weight 

FOCOE H-fLorce 

forward flight induced velocity parameter 

total lift generated by non-rotor bodies 

Diet, horizontal tail 

lift, vertical tail 

lift, wing 

first harmonic (cosine) thrust moment coefficient 
first harmonic (sine) thrust moment coefficient 
Mach number at rotor blade tip 

mass of rotor blade 

total blade thrust moment at specific azimuth angle 
advance ratio 

number of azimuth sectors 

number of blade elements 

rotor rotational velocity 

pressure altitude 

inflow angle, rotor blade radial segment 

inflow angle, rotor blade tip 

power required by rotor 

azimuth angle 

dynamic pressure 

rotor Corque 

radius, rotor blade radial segment 

rotor blade radius 

Reff-e 

Phe Rbat 

effective rotor blade radius (tip loss) 

ambient air density 

location of resultant thrust vector 

solidity 

area, horizontal tail 

area, vertical tail 

area, wing 

Eobor thrust 


184 


$ Global Vars. 


$ COUNT 
called 

$ NAME 

%$ PICK 

$ REGIME 
$ MINUM 
$ MAXUM 
$ INTER 


Seeotructure Vars. 


$ S MATR_VEC 
$ S PERF INPUT 
$ S USER INPUT 


3 Taux auUxdala Gy Viel ieusit 

% temp ambient alr temperature 

% theta eveliegpr ech 

z thetalc first harmonic (cosine) of cyelic patch 
% thetals first harmonic (sine) of cyclic pitch 

% thetao collective pitch at .7 r/R 

$ Tpsi total blade thrust at specific azimuth angle 
$ eg rotor blade taper ratio 

$ twist geometric rotor blade twist 

% Up vertical component of velocity 

% Uptip vertical component of velocity at tip 

% Ut horizontal component of velocity 

$ Uttip horizontal component of velocity at tip 
$ vi induced velocity 

$ Vint forward airspeed 

S Vtip tip speed 

$  wblade weight of rotor blade 


Counter to determine where Performance Input was 


Input .mat file name 

Iteration Method Choice (1-7) 

Include HIGE Calculations Choice (l=yes, 0=no) 
Iteration Start Value 

Iteration End Value 

Iteration Interval 


Matrix/Vector structure 
Perf.m input structure 
user input structure 


%$ Graphics Handle Vars. 


H AL 
H_ANAL 
H_AS 

H ASPECT 
H ASPECT EDIT 
H_BK 

H_BT 
H_BTR 

H check] 
H_check2 
H_check3 
i e038) 

H CIM 
H_CNF 
H_datain 
H dataout 
H_ DISK 


oP 0 oP oP oP OP CP OP OP CP OP OP OCP OP OP OP Of 


Altitude Iteration Radio Button 

Analysis Figure Window 

Airspeed Iteration Radio Button 

Aspect Ratio Static Text Box 

Aspect Ratio Edit Text Box 

Iteration Method << Back Push Button 
Blade Twist Iteration Radio Button 

Blade Taper Ratio Iteration Radio Button 
Save Input Data Check box 

Save Output Data Check box 

Save Matrix & Vector Data Check box 
Change Input Data Radio Button 

Change Iteration Method Radio Button 
Create New Radio Button 

Save Input Data Edit Box 

Save Output Data Edit Box 

Horiz. Tail Under Main Rotor Disk Check Box 
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0 o®® o®® ® HO PH CH CH HP HH HH HP HP KH HP WH HH HK HW H HP WH oH HK HP CH HK CH KH HH HP HP HP HK CP CP 


H_EJANRAD 
H_ EREF 

H_ GO 

H_ GW 

H_ HIGE 
Heeneut tile 
ee 

H IT BOX 

H IT METH 
H_ JAN 

H LB 

H MEN 

H NI 

h_opt 

H OPTIONS 
Beemeputti Ne 
HP 

H PERF IN 
H PERF OUT 
H_PRDA 

H preintin 
Hepes Ent out 
BapErne vec 
H PSCA 

H_ RD 

H_ RES 

H_ RTB 

H_ RUPT 

H SAC 

H_ SOT 

H_ STATUS 

H STATUS1 
H STATUS2 
Hewecdata 
Hevect ie 
H WORK 

H WSA 


Exit JANRAD Radio Button 

Edit/Run Existing File Radio Button 
Analyze Push Button 

Gross Weight Iteration Radio Button 
Iteration Parameters HIGE Check box 

Input File Static Text box 

Iteration Parameters figure window 
Iteration Parameters Static Text Box 
Iteration Method figure window 

JANRAD 98 Figure window 

Input File List Box 

JANRAD 98 Options Menu handle 

No Iteration Radio Button 

Performance Output JANRAD Options Menu 
Options Figure Window 

Output File Static Text box 

Performace Radio Button 

Performance Input Figure Window 
Performance Output Figure Window 

Perform Rotor Dynamics Radio Button 

Print Input File Check Box 

Print Output File Check Box 

Print Matrix & Vector File Check Box 
Perform Stability and Control Radio Button 
Rotor Dynamics Radio Button 

Resume Push Button 

Return to Beginning Radio Button 
Interrupe Push Button 

Stability and Control Radio Button 
Start of Taper Iteration Radio Button 
Top Analysis Status Static Text Box 
Middle Analysis Status Static Text Box 
Bottom Analysis Status Static Text Box 
Save Matrix & Vector Data Edit Box 
Vector File Static Text box 

Working Directory Edit Box 

Wing Span Area Iteration Radio Button 


(Elapsed Time) 
(not used yet) 
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